GLOBAL SEMIGROUP OF CONSERVATIVE SOLUTIONS OF THE 
NONLINEAR VARIATIONAL WAVE EQUATION 



HELGE HOLDEN AND XAVIER RAYNAUD 



Abstract. We prove the existence of a global semigroup for conservative solutions of 
the nonlinear variational wave equation Utt — c(u)(c(u)u x ) x = 0. We allow for initial data 
u\t=o and u t \t=o that contain measures. We assume that < k _1 < c(u) < ft. Solutions 
of this equation may experience concentration of the energy density (uf + c(u) 2 u 2 ,)dx into 
sets of measure zero. The solution is constructed by introducing new variables related 
to the characteristics, whereby singularities in the energy density become manageable. 
Furthermore, we prove that the energy may only focus on a set of times of zero measure 
or at points where c'(u) vanishes. A new numerical method to construct conservative 
solutions is provided and illustrated on examples. 
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1. Introduction 

The nonlinear variational wave equation (NVW), which was first introduced by Saxton 
in |llj . is given by the following nonlinear partial differential equation on the line 

(1.1) utt - c(u)(c(u)u x ) x = 
with initial data 

(1.2) u\ t =o = u , ut\ t =Q = u\. 

The equation can be derived from the variational principle applied to the functional 

J J (u 2 — c 2 (u)u x ) dxdt. 

We are interested in the analysis of conservative solutions of this initial value problem for 
uo,ui £ L 2 (M). It is well known that solutions of this equation develop singularities in 
finite time, even for smooth initial data, see, e.g., [8]. The continuation past singularities 
is highly nontrivial, and allows for several distinct solutions. Thus additional information 
or requirements are needed to select a unique solution, and stability of solutions becomes a 
particularly delicate issue. We here study the conservative case where one in addition to the 
solution u itself, requires that the energy is conserved. For smooth solutions the energy is 
given by £ (i) = f R (u 2 + c 2 u 2 )(t,x) dx. However, as energy may focus in isolated points, one 
has to look at energy density in the sense of measures such that the absolutely continuous 
part of the measure corresponds to the usual energy density. The analysis resembles to 
a large extent recent work done on the Camassa-Holm equation and the Hunter-Saxton 
equation (see, e.g., [HI [13], [10] and references therein). Our main result is the proof of the 
existence of a global semigroup for conservative solutions of the NVW equation, allowing 
for concentration of the energy density on sets of zero measure. 

The NVW equation has been extensively studied by Zhang and Zheng p^l [T3l f!4l IT5l fT6l 
IT71ll8| . However, our approach is closely related to the approach by Bressan and Zheng [5], 
in that we introduce new variables based on the characteristics, thereby, loosely speaking, 
separating waves going in positive and negative direction. 

It is difficult to illustrate the ideas in this paper as there are no elementary and explicit 
solutions available, except for the trivial case where c is constant, which yields the classical 
linear wave equation. Thus one is forced to illustrate ideas numerically. Traditional finite 
difference schemes will not yield conservative solutions, but rather dissipative solutions due 
to the intrinsic numerical diffusion in these methods. Hence it is a challenge of separate 
interest to compute numerically conservative solutions of this equation to display some of 
the intricacies. This question is addressed and analyzed in Section [U 

Let us now turn to a more precise description of the content of this paper. We consider 
the variables R and S defined as 



(1.3) 



R = u t + c(u)u x , 
S = u t - c(u)u x . 
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By (11.1(1 . we have 
(1.4) 

or, on conservative form, 



Rt-cR x = ^-(R 2 -S 2 ), 

S t + cS x = Us 2 -R 2 ), 
4c 



(R^ + S 2 ) t -(c(R 2 -S 2 )) x = 0, 

< £ {R i-S 2 )) t -(R 2 + S 2 ) x = 0. 
c 

Let £(t) denote the total energy of the system at time t, i.e., 

(1.6) £{t)= I (u 2 + c 2 u 2 x ){t,x)dx = f(R 2 + S 2 )dx. 

Jr Jr 

We assume that the initial total energy that we denote £q is finite and that u is bounded 
in L°°. For smooth solutions of (II. ip we have % = 0. We also assume that c 6 C 1 (M) and 
c: R —> [k^ 1 , k] for some constant k > 0. 

From (II. 6|) we see that we need that the functions R and 5 belong to L 2 (M). It turns 
out that, as time evolves, the functions R 2 and S 2 can concentrate on sets of measure zero. 
The example presented in Figure Q], see Section 18.21 illustrates this phenomenon. In this 
case, we have a nontrivial solution u for t nonzero, which is however is identically equal to 
one at t = and ut(0,x) = 0. However, when we analyze this example closer, we see that 
the energy concentrates at the origin, indeed 

lim R 2 (t, x) dx = 5 and lim S 2 (t, x) dx = 25 
t->o t->o 

where 5 is Dirac's delta function. Clearly this complicates the existence and uniqueness 



Figure 1. Plot of u(t, x) for t = -3 (left), t = (center), t = 3 (right). 

question for this equation. As we want to construct a semigroup of solutions for this type 
of solutions, we have to know the location and amount of backward (R 2 ) and forward (S 2 ) 
energy that has concentrated on sets of zero measure, an information which is not given by 
the function u itself. (In our example, since, at t = 0, the function u is identically one and 
its time derivative ut identically zero, we cannot infer where the energy has concentrated.) 
Thus we introduce the set V whose elements, in addition to u, R, S, contain two measures, \i 
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and v, corresponding to forward and backward energy density. More precisely, the measures 
are nonnegative Radon measures that satisfy 

Mac = ^R 2 dx, u ac = ^S 2 dx. 

Our main contribution in this article is to present a rigorous construction of the semi- 
group of conservative solutions in V. Note that the set V is the natural set of solutions 
for conservative solutions, and the semigroup property can only be established in V, as 
illustrated by the example of Figure [TJ Furthermore, by incorporating the energy measures 
as independent variables the formation of singularities is natural, and it allows for more 
general initial data. The present approach also provides a natural numerical method for 
conservative solutions. 

As in [5], the construction of the solutions is achieved via a change of variables into a 
new coordinate system (X, Y) that straightens the characteristics. Even if we use different 
variables, the solutions we obtain are the same, but by extending the solutions to the set P, 
we are able to establish that the solutions we construct satisfy the semigroup property. We 
have to study in details the change of variables mapping — from the original variables to the 
new variables and vice versa — because, in order to prove the semigroup property, we have 
to establish that the two sets of variables match in an appropriate way. Compared to the 
variables used in [5], we prefer variables with a more direct physical interpretation. Namely, 
the variables we are considering are time, t(X, Y), space, x(X, Y), the solution function 
U{X, Y), which formally satisfies u(t(X,Y), x(X,Y)) = U(X,Y) and the energy potentials 
J and K. The definition of the energy potentials J and K follows from (|1.5p . which says 
that the forms \{R 2 + S 2 ) dx + \c{u){R 2 - S 2 ) dt and j^(R 2 - S 2 ) dx + \{R 2 + S 2 )dt 
are closed, so that, by Poincare's lemma, there exist functions, here denoted the energy 
potentials J and K, whose differentials are equal to the given forms. Thus the new set of 
variables we will be considering equals Z = (t, x, U, J, K) and, after rewriting the governing 
equations (11.3(1 and (11.4(1 in the new coordinate system (X, 1"), we get a system of equations 
of the form 

(1.7) Z XY = F(Z)(Z X ,Z Y ) 

where F{Z) : R 5 x R 5 — > R 5 is a bi-linear and symmetric operator, which depends only on 
U, cf. J233|). 

In the new coordinates, the initial data corresponds to the set Tq = {(X,Y) £ R 2 | 
t(X,Y) = 0}. In the smooth case, Tq will be a strictly monotone curve. However, in our 
setting, To may not even be a curve, and even if it is a curve, it may not be continuous nor 
strictly monotone. Indeed, it may contain horizontal and vertical segments, and further- 
more, rectangular boxes corresponding to the situation where both \x and v are singular at 
the same point. If Tq is a curve with no vertical or horizontal parts and the initial data 
is bounded in L°° (by initial data, we mean the values of Z, Zx and Zy on Tq), then 
the existence and uniqueness of solutions to (|1.7p is a classical result, see, for example, [3 
Ch. 4]. In the present paper we have to deal with unbounded data in V (u x and ut are 
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unbounded in L°°). The new coordinates (X, Y) are given by 

dx — c(u)dt = if and only if dY = 



and 



dx + c(u)dt 



if and only if dX = 



that is, the characteristics are mapped to horizontal and vertical lines. We denote by L the 
mapping from the possible initial data in V to the set T defined by To and the value of the 
initial data on T , thus L: V — > T , see Definition 13.81 From T we have to select one curve 
that can be used as initial data for the equation (II. 7|) . There is a certain nonuniqueness 
due to fact that T may not be curve. Let Go denote the set of all curves, including the 
information about the initial data. We let C denote the mapping that from a set Tq selects 
one possible curve, that is, C : T — > Go, see Definition ^. 51 The inverse map that from curve 
determines the corresponding set in T is denoted D, see Definition 13.71 Once we have a 
curve with the initial data, we can in principle compute the solution by solving (II. 7j) . To 
show the existence of a global solution we use the the bi-linearity of (|1.7p and an a priori 
bound on the energy potentials J and K, see Section HI We let the set of all possible 
solutions be denoted by TL, and let S : G — ► TL denote the map that computes the solution 
that passes through the curve in G, see Theorem 14.151 Here G is defined as Go without the 
constraint that t = 0, see Definition 13.21 Recall that as t now is a dependent variable, it 
does not make sense to compute the solution up to a specific time, but rather we determine 
the global solution for all times. Thus we need a mapping that extracts the solution Z for a 
given time T, that is, the intersection of the solution in TL with the set where t(X, Y) = T. 
Let E : TL — » Go denote the map that from any given solution in TL extracts the solution at 
t = 0, that is, in Go, see Definition 15.11 Next we define the operator tT'.H—*H that shifts 
time in a solution in TL by a given time T, see Definition 15.21 Now we can define the map 
St '■ J 7 — ► T by St = D o E o t T o S o C, see Definition 15.41 A key result is that St is a 
semigroup on T , see Theorem 15.51 Next we need to return to the original variables. Let 
M: T — > V denote that map, see Definition l6.il Thus the solution operator St - T> — ► V is 
defined by (Definition 16.4(1 



It remains to show that St is a semigroup. However, since M is not inverse of L, as 
LoM/ Idjr, the semigroup property of St still does not follow from (|1.8p . This fact is 
explained as follows. When changing variables, we have introduced a degree of freedom 
that we now want to eliminate. This degree of freedom can be identified precisely with 
the action of the group G 2 , where G denotes the group of diffeomorphisms of the real line. 
Indeed, by simply using the bi-linearity of (|1.7|) . one can check that if Z is a solution to 
(|1.7p . then Z(X,Y) = Z(f(X),g(Y)), where (/, g) G G 2 , is also a solution to the same 
equation. The transformation (X,Y) (f(X),g(Y)) corresponds to a stretching of the 
plane M 2 in the X and Y directions. Note that this transformation maps horizontal (resp. 
vertical) lines to horizontal (resp. vertical) lines and therefore preserves the directions of the 
characteristics. Moreover, this transformation does not affect the solution in the original 
coordinates. To illustrate this we ignore for the moment for the sake of simplicity, the 



(1.8) 



St = M o St ° L. 
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energies \i and v in the definition of V. The solution u(t, x) can be seen as the surface in M 3 
given by (t, x,u(t, x)) where (t,x) G K 2 are parameters. Through our change of variables, 
we obtain another parametrization of the same surface, namely, 

(1.9) (t(X,Y),x(X,Y),U(X,Y)) 

where (X, 1") e R 2 are the new parameters. Additional properties of the solution Z = 
(t, x, U, J, K) which are contained in the definition of H guarantee that the surface defined 
by (|1.9p does not fold over itself so that it is in fact a graph. It is then clear from (|1.9p 
that the transformation (X,Y) t— > (f(X),g(Y)) is simply a re-parametrization of the same 
surface, which defines u(t, x) uniquely. At the level of the set T , which corresponds to a 
parametrization of the initial data in the new coordinates, we can also define the action of 
the group G 2 that denote ip X (f,g) ■-► ip ■ (f,g) for any ip G T and (/, g) G G 2 . We prove 
that two elements which are equivalent correspond to the same element in T>, that is, 

(1.10) M(ip) = M(ip) 

where ip = ip ■ (f,g) for some (/,<?) G G 2 . From (jl.lOj) . it is now clear why LoM / Idjr 
as, in general, ip and ip are distinct. We introduce a subset T§ of T which corresponds 
to a section of T with respect to the action of the group G 2 , which means that the set 
Tq contains only one representative of each equivalence class so that T '/G 2 and T§ are in 
bijection. The system (II. 7p preserves the strict positivity of the quantities xx + Jx and 
xy + Jy and the set T somehow inherits this property which makes it possible to define 
the projection II: T — > Tq. The projection II associates to any element in T its unique 
representative in T§ which belongs to the same equivalence class. As expected, since we 
have now eliminated the degree of freedom we introduced by changing variables, we obtain 
that To and V are in bijection. We are then able to prove that St is a semigroup. 

Our main result, Theorem 17.91 reads as follows: 
Theorem. Given (uq, Ro, So, /io, ^o) £ let us denote (u, R, S, n, v)(t) = St(uo, Ro, So, fJ-o, vq). 
Then u is a weak solution of the nonlinear variational wave equation (jl.ip . that is, 

(1.11) / {<t> t - (c(u)(j>) x )Rdxdt+ / (<k + (c(u)<f>) x )Sdxdt = 

for all smooth functions (j) with compact support and where 

(1.12) R = m + c(u)u x , S = ut — c(u)u x . 

Moreover, the measures fj,(t) and u(t) satisfy the following equations in the sense of distri- 
bution 

(1.13a) {n + v) t - (c(n - v)) x = 
and 

(1.13b) {- c { i i-v)) t -{ i i + v) x = Q. 

The mapping St '-T> —*T> is a semigroup, that is, 

St+t' = St° s t > 
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for all positive t and t . 



Furthermore, we note the following important result (Theorem 17.10 
Theorem. The solution satisfies the following properties: 



(i) For allteR 
(1.14) + "(*) ( R ) = Mo(K) + ^o(K)- 



(ii) For almost every i£l, the singular part of n(t) and v(t) are concentrated on the 
set where d(u) = 0. 

In this article, we do not study the stability of the solutions. However, since the solutions 
we obtain coincide with the ones obtained in |5] for initial data which do not contain any 
singular measure, the solutions in that case also satisfy the stability result stated in [U 
Theorem 2]. To obtain a continuous semigroup of solution in V, we would like to follow 
the approach developed in [9], [4] for the Camassa-Holm equation and Hunter-Saxton 
equations. In these two papers, the conservative solutions are also obtained via a change of 
variables which is invariant with respect to relabeling (i.e., with respect to the action of the 
group of diffeomorphisms G). We define a distance between equivalence classes in the new 
coordinates. This distance is then mapped back to the original set of coordinates so that we 
obtain a continuous semigroup for this metric. In the case of the nonlinear wave equation, 
in particular, because of the truly two dimensional nature of the problem, it is not so easy 
to formulate a stability result in the new coordinates which holds when mapping back to 
the original set of variables. In Lemma f4. 121 we present a result in that direction. 

There is a lack of explicit solutions to NVW. In this paper we consider two explicit 
examples. The first example, see Section 18.11 is the simplest possible, namely the linear 
wave equation (with c constant), but with general initial data. We recover as expected the 
familiar d'Alembert solution. The energy measures are transported with velocity ±c. The 
second example, see Section f8T2l is a truly nonlinear case with velocity given by (|1.19|) . 
However, here we choose the simplest nontrivial initial data with energy concentration 
initially for both measures. The corresponding equation fj 1 . T[) is solved numerically, and 
the result is illustrated on Figs. [TJ I9HT21 

The numerical method that yields conservative solutions is described in Section [U 

1.1. Physical motivation for the nonlinear variational wave equation. The NVW 

equation was first derived in the context of nematic liquid crystals, see pTj HO]- More 
precisely, a nematic crystal can be described, when we ignore the motion of the fluid, by 
the dynamics of the so-called director field n = n(x, y, z, t) G M 3 describing the orientation 
of rod-like molecules. Thus |n| = 1. The Oseen-Franck strain-energy potential is given by 



where a, (5, 7 are constitutive constants. Consider next the highly simplified case of director 
fields of the type 



(1.15) 



VF(n, Vn) = a |n x (V x n)| 2 +/?(V-n) 2 + 7 (n- V x n , 



(1.16) 



n = n(x, t) = cos(u(i, x))e x + sin(u(i, x))e. 



■y 
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where e x and e y are unit vectors in the x and y direction, respectively. In this case the 
functional W(n, Vn) vastly simplifies to 

(1.17) W(n,Vn) = (/3cos 2 u + a sin 2 u)u 2 x , 

l 1 2 9 

and |nt| = uf. The dynamics is described by the variational principle 

(1.18) ^ J J (uf - c 2 {u)u 2 x )dxdt = 0, 
where 

(1.19) c 2 ( u ) = /?cos 2 u + a sin 2 u, 
which results in the nonlinear variational wave equation 

(1.20) u u - c(u)(c(u)u x ) x = 0. 

2. Equivalent system for the NVW equation 

In this section, we assume the existence of a smooth solution u = u(t,x) to (jl.ll) . We 
introduce the change of variables (t, x) i— ► (X, Y) which straightens out the characteristics: 
The forward characteristics, which are given by the solutions of ^ = c(u(t, x(t))), are 
mapped to the horizontal lines while the backward characteristics, which are given by the 
solutions of -fa = —c(u(t,x(t))), are mapped to the vertical lines. Formally, we can rewrite 
these conditions as 

(2.1) dx - c{u)dt = if and only if dY = 
and 

(2.2) dx + c{u)dt = if and only if dX = 0. 

Our goal now is to rewrite the governing equation (|1.1|) in terms of the new variables (X, Y). 
The variables (t, x) become functions of (X, 1") that we denote t(X,Y) and x(X,Y). We 
set 

(2.3) U(X,Y) =u(t,x). 
Since 

dx = xxdX + xydY and dt = txdX + tydY, 
we obtain from (|2.1I) and (|2.2I) that 

(2.4) x x = c(U)t x and xy = -c(U)ty. 

From IjOll . we infer that the forms \(R 2 + S 2 )dx + f (R 2 - S 2 )dt and ±(R 2 - S 2 )dx + 
\(R 2 + S 2 )dt are closed. Therefore, by Poincare's lemma, we infer the existence of two 
functions J and K for which these forms are the differentials, that is, 

(2.5) dJ =^{R 2 + S 2 )dx + ^{R 2 -S 2 )dt 
and 

(2.6) dK = ^-(R 2 -S 2 )dx + ^(R 2 + S 2 )dt. 
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We have, after using (|2.4 



dJ = ^(R 2 + S 2 )dx + ^(R 2 - S 2 )dt 

= ^(R 2 + S 2 )(x x dX + xydY) + ^{R 2 - S 2 )(t x dX + t Y dY) 

= ^R 2 x x dX + -S 2 x Y dY, 
and, similarly, we get 

R 2 , S 2 

dK = — x x dX xydY 

2c 2c 

so that 

(2.7) J x = c{U)K x and J Y = -c(U)K Y 

hold. Note the similarity between the relations (12.71) for the pair (J, K) and the relations 
(|2.4p for the pair (t, x). We want to compute the mixed second derivatives of our new 
variables, namely, t, x, U, J and K. By (12.41) . we obtain 

dt = t x dX + t Y dY = -x x dX - —x Y dY. 

c c 

By expressing the fact that the form dt is closed (since it is exact), we get 

d (\ \ d (\ 

t 

which implies 
Similarly, since the form 



dY \c dX 



XX = -77^7 ~ x y 



c 



x X y = 7T (uyxx + u x x Y ) . 
2c 



dx = xxdX + x Y dY = ctxdX — ct Y dY 

is closed, we obtain 



which implies 



d . . d 

W (ctx) = -dx {ctY) 



t X Y = ~7T i u xt Y + U Y t X ) ■ 
2c 



By using the relations (12.71) . the form dK can be rewritten as 

dK = -J x dX - —J Y dY 

c c 

and, expressing the fact that dK is a closed, we obtain 

W (~c Jx ) = JX \c JY 



which yields 



Jxy = y c ( j xUy + JyUx) ■ 
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Similarly, we can rewrite the form dJ as 

dK = cK x dX - cKydY 
and, expressing the fact dK is closed, we get 

Kxy = -r- (K X Uy + K Y Ux). 
2c 

Let us consider the forms 

R 1 

(2.8) ui = —dx + -Rdt 
v ' 2c 2 

and 

S 1 

(2.9) uj 2 = —dx - -Sdt. 
In the new variables, these forms rewrite 

wi = -5— — -(x x dX + xydY) + + cu x )(t x dX + tydF) 
2c 2 



= (litix + u x xx)dX (after using (|2.4 

(2.10) = l/ydX, 
and, similarly, we find 

(2.11) u 2 = -U Y dY. 
From (12. 8|) . by using (|1.4|) . we obtain 



R d R 1 

gL>i = ( oUt)dt A dx H — A (it 

2c 2c z 2 

-R+ — ci? r , , c'-R 1 , „ , 
= ^ dt Adx- + S)dt A 

c'(i? 2 - S 2 ) _ dRl 
8c 2 " 2c 2 2 

1 ^ 2 



(R + S) )dt A 



2~2 11^ + dx Adt = Y2 U t dX Adt ' 

and, furthermore, we obtain 

c' / 1 1 
duji = — 77 [ -(R 2 + S 2 )dx Adt + -RSdx A dt 
2c z V4 2 

d d 
(2.12) = -^d J A dt - —uJi A lo 2 



because 



RS , 

o>i A CJ2 = dx A dt. 

2c 
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We rewrite (|2.12|) in the new set of variables 

d d 
dui = —^{J x dX + JydY) A (t x dX + tydY) - — w x A uj 2 
2c L 2c 

d d 

= AJxXy + Jyx x )dX A dY H U x UydX A dY. 

2c 6 2c 

At the same time, by (I2.10|) . we have duj\ = —UxydX A dY , and therefore it follows that 

d d 
U X y = ^3 {Jxxx + Jyxx) - Y c UxUy ' 

Finally, we obtain following system of equations 



(2.13a) 


txY — 


-^(Uxty + U Y t x ), 
2c 


(2.13b) 


XXY = 


c < 

— {Uyx x + U x xy) , 


(2.13c) 


UxY = 


d 

7^3 (xyJx + Jyxx) - 


(2.13d) 


JXY = 


|- (JxUy + JyU X ) , 
2c 


(2.13e) 


KxY = 


-^(KxUy+KyUx) 
2c 



d 



2c 



Let Z denote the vector (t, x, U, J, K). The system (|2. 13|) then rewrites as 
(2.14) Z X y = F(Z)(Z x ,Zy) 

where F(Z) is a bi-linear and symmetric tensor from M 5 x M 5 to M 5 . Due to the relations 
(|2.4p . either one of the equations ()2. 13a|) and (I2.13bp is redundant: one could remove 
one of them, and the system would remain well-posed, and one retrieves t or x by using 
(|2.4p . Similarly, either one of the equations (|2.13dj) and (I2.13ep becomes redundant by 
(|2.7p . However, we find it convenient to work with the complete set of variables, that is, 
Z = (i, x, U, J, K). We will see later that the solutions of the system (I2.13P preserve these 
conditions. 

To prove the existence of solutions to (|2.13l) . we use a fixed point argument. The argu- 
ment is similar to the one that can be found for example in and in [5]. However, in order 
to take into account the non-regularity of the data {uq x and «oi are in I? and the energy 
can concentrate on sets of zero measure), we have to consider, in the new set of coordinates, 
data given on curves which have parts which are parallel to the characteristic directions. In 
particular, the curves are not given as graphs of a function. We are looking for a solution 
that satisfies a given initial condition at time t = 0. In the (X, Y) plane, the set of points 
which correspond to initial time, that is, t(X, Y) = 0, may be a curve, {X(s),y{s)) G M 2 , 
parametrized by s G R, but it may also be a more complicated set, see Figure [3] that we 
will comment on later. We consider curves of the following type. 
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Definition 2.1. We denote byC the set of curves in the plane M. 2 parametrized by (X(s),y(s)) 
with s£R, such that 

(2.15a) X — Id, y - Id G W 1 ' 00 ^), 

(2.15b) ^ > 0, 3> > 

ane? £/ie normalization 

(2.15c) ^ + = s, /or all s G R. 

We se£ 

(2.16) ll(^,^)|| c = ||Af - Id\\ LOO + ||^ - Id|| £00 • 

From the initial data (uo, Rq, So), we want to define the curve To = (X(s),y(s)) in 
C which corresponds to the initial time and the value of Z on this curve. To solve the 
governing equations (|2.13p . we need to know the values of Z, Zx and Zy on the curve IV 
In total, we have to determine 17 unknown functions. Given the initial data (uq, Ro, So), 
there is no unique way to define the curve (X(s),y(s)) and the values of Z on this curve 
in order to obtain to the desired solution. This fact is due to the relabeling symmetry, 
a degree of freedom which is embedded in the set of equations (|2.13l) that we precisely 
identify in Section [71 For now, the goal is to use this degree of freedom to construct 
an initial data which is bounded in L°°(M) on the curve. Let us now explain how we 
proceed for an initial data (uo, Rq, So) G [L 2 (M)] 3 for which energy has not concentrated 
and we will see later how to extend this construction to initial data containing singular 
measures. In this case, the function X and y are invertible and, slightly abusing notation, 
we denote by Z(s), Z X (X) and Z Y (Y) the values of Z(X(s),y(s)), Z X (X, y(X~ 1 {X))) 
and Zy (X (y^ 1 (Y)), Y), respectively. By definition, we have 

(2.17) t(s) = 0, 
and, naturally, we set 

(2.18) U(s)=uo(x(s)). 

From the formal derivation of the previous section, we have the following relations 

(2.19) J X (X) = c(v)K x (X) = ^(x)xx(X), JyO>) = -c(u)K Y (y) = \sHx)x Y (y), 

We have the compatibility condition 

(2.21) Z(s) = Z x (X(s))X(s) + z Y (y( s ))y(s). 

We have 17 unknowns (X,y,Z,Z x ,Z Y ) and 15 equations, namely (l2~T7V(l2~2Tl) . (l2~T5cll 
and (12. 4p . We use the two degrees of freedom that remain in order to obtain Z x and Zy 
bounded. We set 

(2.22) 2x x (X) + J X (X) = 1 and 2x Y (Y) + J Y (Y) = 1. 



THE NONLINEAR VARIATIONAL WAVE EQUATION 13 

Since xx and Jx are positive, it follows from (I2.22p that these two quantities are bounded. 
From the fact that 2xxJx = (c(U)Ux) 2 , it also follows that Ux is bounded so that Zx is 
bounded. The same conclusion holds for Zy. The normalisation (I2.22p is convenient but 
arbitrary, see Section [7l In particular, the coefficient 2 in front of xx and xy in (12, 22ft 
does not have any importance; it is used here to make the definition compatible with the 
normalization we will introduce in Section [3] for the general case. From (|2.21l) . (I2.20P and 
(gUID , we get 

(2.23) XX ( X ) = —?-p(x) and x Y (y) = — ^(x), 

(2.24) Jx(X) = -K X {X) = ^(x), Jy(y) = --Ky(y) = ^>(x), 

(2.25) u x {x)= , A 2R ° n2 A x), u Y (y) 2 ' s " ' 



c(4 + lZ§) v " c(4 + Sg) w 

Equation (|2.17|) implies 

o = t x (x)x + t Y (y)y = x x (x)x - x Y {y)y 

and, at the same time, we have by the chain rule 

x( s ) = x x (x)x + x Y (y)y 

and therefore 

(2.26) x x (x)x = x Y (y)y = ^. 

Hence, by (I2.15cp . (I2~23ll and (12361) . we get 
and we define x(s) implicitly as 

(2.27) 2x(s)+ / + S^)cfa = 2s. 

J — oo 4 

We have 

2x + j = 2x x (*)* + Jv(^)^ + 2xy(^)^ + Jy(W = 2 
because of (|2.21l) and (j2 . 15c() so that 2x + J = 2s. Hence, 

(2.28) J{S) = 8J_ ( R o + S %) dx 
and 

(2.29) ^g^«fc, 
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which are also defined as the integrals of the forms dJ and dK given by (12, 5p and (12,61) on 
the line (t,x) = {0} x (-oo,x(s)). From ([2361) and (j2~23l) . it follows that 

(2.30) *( a ) = x(s)(l + ij2g)(s( S )) and = x(s)(l + J^X^)), 
and we set 

(2.31) Af (a) = + - / R% dx and = x(s) + - / S" 2 dx. 

3. The initial data 

In order to construct a semigroup of conservative solutions, we have to take into account 
the part of the energy which has concentrated in sets of measure zero and we need to 
consider initial data in the set V that we now define. 

Definition 3.1. The set V consists of the elements (u,R,S,n,v) such that 

(u,R,S) € [L 2 (R)] 3 , 
u x = ^.{R — S) and [i and v are finite positive Radon measures with 

1 1 

(3.1) /x ac = -R 2 dx, i/ ac = -S" 2 dx. 

The measures [i and v correspond to the left and right traveling energy densities, respec- 
tively. Given the initial data (uq, Rq, Sq, Ho, fo), we have defined an element in Qq where 
the set Qq is defined below and which correspond to a parametrization of the initial data in 
the new system of coordinates. Elements of Q consists of a curve (X(s),y(s)) (for Qq, this 
curve corresponds to time equal to zero) and three variables, Z, V and W, that we now 
introduce. These functions correspond to the data that matches the solution Z to (12 . 13[) 
on the curve [X ,y) in the sense that 

(3.2a) Z{s) = Z(X(s),y(s)) 

and 

(3.2b) V{X{s)) = Z x {X{s),y{s)) and W(X(s)) = Z Y {X{s),y{s)) 

It is then convenient to introduce the following notation: To any triplet (Z, V, W) of five 
dimensional vector functions (we write Z = {Z\, Z2, Z%, Z^, Z§), etc), we associate the 
triplet (Z a ,V a ,W a ) given by 

(3.3a) Z$ = Z 2 - Id, V 2 a = V 2 -i W 2 a = W 2 -^ 

and 

(3.3b) Z? = Zi, Vf = V t , wt = m 

for t G {1,3,4,5}. 
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Definition 3.2. The set Q is the set of all elements which consist of a curve (X(s),y(s)) 
and three vector valued functions from R to R 5 denoted Z(s), V(X), W(Y"). We denote 
G = (X,y,Z,V,W) and set 

(3-4) \\&\\ g = \\U\\ LHR) + \\V a \\ L2 + \\W a \\ L , 

where we denote U = Z3 and 

(3-5) Hieilb = ||(^y)|| C + Wy^J-h^R) + \\ m+m \\L°°iR) 

+ \\z a \\ Lac + \\V a \\ Laa + \\W a \\ Lao . 

The element 6 Q if 
(i) 

||S||g < 00 and |||@|||g < 00; 

(") 

(3.6) V 2 ,W 2 ,V 4 ,>V4 > 0; 

(iii) for almost every s, we have 

(3.7) Z{s) = V(X(s))X(s) + W(y(s))y(s); 

(iv) for almost every X and Y , we have 

(3.8a) 2V 4 (X)V 2 {X) = [c{U)V z {X)f , 2W 4 {y)W 2 {y) = (c(^)W 3 (^)) 2 , 

(3.8b) v 2 (x) = c(u)V!(x), w 2 (y) = -c{u)m{y), 

(3.8c) Va{X) = c(U)V 5 (X), Wi{y) = -c(U)W 5 (y). 

(v) We require 

(3.9) lim JO) = 

where we denote J(s) = Z±(s). 
We denote by Qq the subset of Q which ■parametrizes data at time t = 0, that is, 

Go = {&£G\Z 1 = 0}. 



The requirement (13. 9p corresponds to a normalization of the energy potential (or cumu- 
lative energy) to zero at minus infinity. The variables Z, V and W are not independent of 
one another as it can be seen from (13. 7j) . (|3.8b|) . (I3.8cl) but selecting a set of independent 
variables will require an arbitrary choice that we prefer to avoid and that is why we consider 
all the variables at the same level. For £ 5o, we get by using (13.71) and (I3.8bl) . that 

(3.10) v 2 (x(s))x( s ) = m(y(s))y{s). 

By using the normalization (j2. 15c|) . we obtain that 

,3.11, * = nn , y- mx) 



v 2 (x) + w 2 (y)' v 2 (x) + w 2 (y) 
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{x ,Y ) 



r = {(X,Y) | t{x^ry= 0} 

T<X 2 ,Y 2 ) 
(X 2 ,Y 2 ) 



Figure 2. The domain t(X, Y) = in the X, y plane consists of the union 
of a graph of a strictly increasing function, vertical and horizontal segments 
and rectangular boxes. 



and, in principle, by integrating (|3.11|) . we recover X and y. However, there are two 
obstacles to that: The function V2 and W2 are in general not Lipschitz so that we cannot 
use the standard existence theorems for the solutions to (|3.11|) and, in addition, both V2 
and W2 may vanish (it is what happens in the case of a box) and (13. lip does not make 
sense any more. Given (uo, Ro, So, no, vq), in the case where fio = (Mo)ac an d vo = ( z/ o)aci 
we have defined 9 = (X,y,Z,V,W) G Go by pTTl) . (TTTHH . (12371) . f238l) . (I2~2%1) . (12331) . 
(I23TD and 

V 4 (*(«)) = c(£/)V 5 (*( S )) = j^L(x(s)), 

= -c(U)W 5 (X(s)) = j^s(x(8)), 

We do not prove here that, for this definition, we indeed have (X, y, Z, V, W) G Go because 
it will be done later in more generality, see Definitions 13.81 and 13.51 In the next section 
we consider {X,y, Z, V, W) £ ^0 and construct solutions of (I2.13P which satisfy (13.21) . 
However, the set Go is n °t adequate when it comes to parametrize initial data. In the case 
where there is no concentration of the measures, that is, no = (^o)ac and uq = (Vo)ao we 
can see from (12 .3 1 1) and (|2.27l) that X > and y > almost everywhere so that the curve 
does not contain strictly vertical or horizontal regions. This property is not preserved by 
the equation. In particular it means that at a later time, say T > 0, we can find a curve 
(X,y) G C such that t(X(s),y(s)) = T and X(s) = or y(s) = on an interval [sj,s r ], 
with si < s r . In general, the set of points 

T T = {(X,Y) G M 2 I t{X,Y) = T} 

is not a curve but a domain which consists of the union of a graph of a strictly increasing 
function, vertical and horizontal segments and rectangular boxes, see Figure [H If Ft 
contains regions with boxes or vertical or horizontal lines, it means that part of the energy 
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of the solution is concentrated at time T in sets of zero measure, see Section [6l We want to 
parametrize domains Tq (or IV) depicted in Figure [3j which give the solution at time zero 
(or a given time T) and which may contain boxes. The set Go defined above is inappropriate. 
When considering an element in Go, we choose a curve and in the case of a box, the choice 
of the curve which joins the two diagonal corners of the box while remaining inside the box 
is arbitrary. Thus we introduce an unwanted degree of freedom in the parametrization of 
the initial data. The domain Tq depicted in Figure [3] can be parametrized by using two 
nondecreasing functions x\(X) and x\{Y) and by considering the set {(X, Y) G R | x\(X) = 
x 2 (Y)}. Such sets consist exactly of the union of the graph of a strictly increasing function 
(when x[ > and x' 2 > 0), a horizontal segment (when x'^X) = for X G [X ,X ] and 
x' 2 (Yq) > 0), a vertical segment (when x'^Xi) > and x' 2 (Y) = for Y G [Yi,Yi]) and a 
rectangular box (when x[(X) = x' 2 (Y) = for X G [X 2 ,X 2 ] and Y G [I2, $2]), see Figure H 
This observation (partially) justifies the definition of the set T which is introduced below. 
The set T can be considered as a consistent way to parametrize initial data. However, to 
construct the solutions, we need to choose a curve and we use the description of the initial 
data given by Go so that, finally, both sets are needed. To define T ', we have to introduce 
the group G of diffeomorphisms with some regularity conditions. 

Definition 3.3. The group G is given by all invertible functions f such that 



We can now define the set T . 
Definition 3.4. We define the set T consisting of all function ip = (^1,^2) such that 



(3.12) 



/ — Id and f 



-1 



Id both belong to W 1 ' co (R), 



and 



(/ - Id)' G L 2 (R). 




(3.15) 



j[ = c{JJx)K[, J' 2 



c{U 2 )K' 2 



(3.16) 



x[j[ = (c(^)Vi) 2 , x' 2 J' 2 = {c{U 2 )V 2 )\ 
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(3.17) Xl + Jl, x 2 + J 2 eG, 

(3.18) lim Ji(X)= lim J 2 (Y) = 

X^— oo Y— +— oo 

ana 1 , /or any curve (X, y) G C shc/i i/iai 

xi(Af(s)) = x 2 (^(s)) /or a// sGl, 

we /lawe 

(3.19a) = 
/or nil s £ M and 

(3.19b) C/((^(s))^( S ) = E^(y(a))y(a) = Vi(Af(a))*(a) + F 2 (^(s))^(s) 

for almost all s G R. 

We show in Section [5] that, given a solution Z of (I2.13p . there exists a unique element 
ifi G ^ which describes in a unique way the set To = {(X,Y) G M 2 | f(X, y) = 0} 
and the values of Z, Zx and Zy on this set. The functions x\ and x 2 define the set To by 
To = {(X, Y) G M 2 | xi(X) = x 2 (Y)}. It means in particular that, for any curve (X,y) G C 
such that x\{X(s)) = x 2 (y(s)), we have t(X(s),y(s)) = 0. The functions U\ and U 2 give 
the value of U(X,Y) on the set To, as a function of X and Y. To be more concrete, 
let us consider the example where x\ and x 2 are smooth, invertible and the inverses are 
also smooth. In that case, which in fact corresponds to the case where uq, Rq and So are 
smooth and there is no concentration of energy, i.e., = (/ i o)ac and vo = (^o)ac, the set To 
is the graph of a strictly increasing function (there is no rectangular box and no vertical or 
horizontal segments). The curve To is given by either Y = x^ 1 o x 2 {X) or X = x^ 1 o X\{Y) 
and, just for this paragraph, for the sake of simplicity, we denote Y(X) = x^ 1 o x 2 (X) and 
X(Y) = x^ 1 o xi{Y). Then, we have 

(3.20) Ui{X) = U(X,Y(X)) and U 2 (Y) = U{X(Y),Y). 
The functions of V\ and V 2 give the partial derivative of U . We have 

(3.21) Vi{X) = Ux(X,Y(X)) and V 2 {Y) = U Y (X(Y),Y). 

As we can see in this example, the functions Ui, U 2 , V\ and V 2 are not independent from 
one another, and the way they depend one another is given by (|3. 19a|) and (|3.19bj) . The 
function J\{X) gives the amount of forward energy contained on the curve Y = Y(X) 
between —00 and X, that is, 

rX 



JipQ = f J x (X,Y(X))dX. 

J —00 



In the original set of coordinates, it gives Ji(X) = \ J^i^ R^dx. Similarly, the function 
J 2 (Y) gives the amount of backward energy which is contained on the same curve between 
—00 and y, that is, 

r-Y 



MY) = I J Y (X(Y),Y)dY. 
J —00 
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In the original set of coordinates, it gives J 2 (Y) = \jZ 2 ^ ] Sldx. We recall that these 
expressions hold only for smooth initial data with no concentration of energy. Still in this 
case, the functions X\ and x 2 are strictly increasing so that x' x > and x' 2 > and the 
conditions (|3.16h entirely determine the energy densities, which are given J[ and J 2 in the 
new sets of coordinates. We have 

1 J' 
= (^o)ac = ~;Ro(x) dx = -±ox^ l (x)dx 
4 x\ 

and the corresponding expression for vq. In the case where there is concentration of energy, 
the functions x[ or x' 2 vanish. The set where x\ (respectively x 2 ) vanishes corresponds to 
the region where the energy density fio (respectively vq) has a singular part. On those sets, 
the energy densities J[ and J 2 cannot be retrieved from (|3.16p . It is consistent with the 
fact that the singular parts of the energy fi and v cannot be recovered by the knowledge 
of the function u, R and 5, as illustrated in the example presented in the introduction. 
As we will see in Section [HJ the relations (|3.16l) correspond to a reformulation in the new 
coordinate system of (|3.1I) . 

We define a mapping C which to any given initial data ip £ T associate the corresponding 

data 9 = (x,y,z,v,w) eg - 

Definition 3.5. For any ip = (ipi,ip2) £ F , we define 

(3.22) X(s) = sup{X G R | xi(X') < x 2 {2s - X') for all X' < X} 
and set y(s) = 2s — X(s). We have 

(3.23) x 1 {X{s)) = x 2 {y{s)). 
We define 

(3.24a) t(s) = 0, 

(3.24b) x(s) = xx(X(s)) = x 2 (y(s)), 

(3.24c) U(s) = U^Xis)) = U 2 (y(s)), 

(3.24d) J(s) = Ji(X(s)) + J 2 (y(s)), 

(3.24e) K(s) = K^X^)) + K 2 (y(s)) 

and 

my) = 

m(Y) = V 2 (Y), 

mix) = j 2 (y), 

W 5 (Y) = K' 2 (Y). 

Let C be the mapping from T to Qq which to any ip S T associates the element (X, y, Z, V, W) 
defined above. 



Vx(X) 


1 


2c(U 1 (X)) 


V 2 (X) 


= ^(1), 


V 3 (X) 


= V 1 (X), 


V 4 (X) 


= Ji(X), 


V 5 (X) 


= K[(X), 
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Proof of the well-posedness of Definition 13.51 Let us prove that X is increasing. Given s > 
s, we consider a sequence X* which converges to X(s) with Xi < X(s). We have xi(Xi) < 
X2(2s — Xi) which implies x\(Xi) < X2{2s — X-i) because X2 is increasing. Hence Xi < X(s). 
By letting i tend to infinity, we get that X(s) < X(s). By continuity of x\ and X2, we have 
xi(X(s)) = X2(y(s)). We claim that X is Lipschitz with a Lipshitz constant no bigger than 
2, i.e., 

(3.25) \X(s) - X(s)\ <2\s - s\. 

Let us assume without loss of generality that s > s. If (I3.25P does not hold, we have 

(3.26) X(s)-X(s)>2(s-s) 

for some s and s in R. It implies y(s) < y(s). Then, by monotonicity of X2, 

Xl (x{ s )) = x 2 (y(s)) > x 2 (y(s)) = x 1 (x(s)), 

and therefore x\(X(s)) = x\(X(sj) because x\ is an increasing function and X(s) < X(s). 
It follows that x\ is constant on [X(s),X(s)]. Similarly, one proves that X2 is constant on 
[y(s),y(s)\. Let us consider the point (X,Y) given by Y = y(s) and X = 2s - y(s). We 
have 

X(s) = 2s- y(s) < X < 2s - y(s) = X(s) 
so that x\(X) = x\{X{s)) = X2(y(s)) = X2(2s — Y) and X < X(s), which contradicts the 
definition of X(s). Hence, (|3.26l) cannot hold and we have proved (|3.25|) . Let us prove that 
X — Id £ L°°. We have 

x( s ) - s = l -{x(s) - y{ s )) = l -(x(s) - x 1 (x(s)) + x 2 (y(s)) - y(s)) 

which is bounded as x\ — Id and X2 — Id belong to L°°. Let 

B = {s G R | X{s) > 1}. 
Since X + y = 2, we have y > 1 on B c . Hence, 

/ U 2 (s)ds= [ U 2 (s)ds+ [ U 2 (s)ds 
Jr Jb Jb c 

< [ uf(x( s ))x(s)ds+ [ u 2 2 {y(s))y(s)d S 

Jb Jb c 
^ ll^illis + ll^lll 2 • 

It is then straightforward to check that the remaining properties that enter in the definition 
of Go are fulfilled by (Z,V,W). To check that (|3.17|) is fulfilled, we use Lemma [3761 which 
is stated below. □ 

Lemma 3.6. Let a > 0. // / satisfies (I3.12p . then 1/(1 + a) < f^ < 1 + a almost 
everywhere. Conversely, if f is absolutely continuous, f — Id G L°°(R) and there exists 
c > 1 such that 1/c < < c almost everywhere, then f satisfies (13 . 12[) and 

11/ _ IdHjy-i.ooQuf) + ||/ 1 - Id | | w i.oo(k) - a 
for some a depending only on c and \\f — ld\\ Lao m\ . 
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The proof of this short lemma is given in [9]. In the opposite direction, to any element 
(X,y,Z,V,W) G Go, there corresponds an element (i/^i , 1^2) G F given by the mapping D 
that we define next. 

Definition 3.7. Given (X,y,Z,V,W) G Go, let V>i = Oi, U Xl J x , K u V x ) and ip 2 = 

(x 2 , U2, J2,K2, V 2 ) be defined as 

(3.27) Xl (X(s)) = x 2 (y(s))=x( S ) 
where we denote x(s) = Z 2 {s) and 

(3.28) = WOO) = U(s) 
where we denote U(s) = Z 3 (s) and 

(3.29) Jx{X[s))= [ V 4 (X(s))X(s)ds, WOO) = / W 4 (y(s))y(s) ds, 



-00 </ —00 

s /-s 



(3.30) K 1 {X{s))= V 5 (X(s))X(s)ds, K 2 (y{s)) = W 5 (y(s))y(s)ds, 

J— 00 J —00 

and 

(3.31) Vi = V 3 , V2 = W 3 . 

We denote by D i/ie mapping from Go t° J~ which to any (X, y, Z, V, W) G £?o associates 
the element ip as defined above. 

Well-posedness of Definition \3. 7l Since t(s) = we have 

= i( s ) = Vi(x( s ))x(s) + Wi(y(s))y(s) 

which implies that 

(3.32) v 2 (x(s))x( s ) = w 2 (y(s))y(s) 

by (|3.8bp . We check the well-posedness of (I3.27P and (|3.28|) . Let us consider s and s such 
that X(s) = X{s). Since X is increasing, it implies X(s) = and y{s) = 2 for all s G [s, s]. 
From (l3~52l . it follows that W 2 (3 ; (s)) = for all s G [s,s]. Hence, 

x(s) = v 2 (^(s))^(s) + w 2 (y{$))y(s) = 

and x(s) = x(s) so that the definition (I3.27P is well-posed. For s G [s,s], we have 
W 3 (y(s)) = 0, by (l3T8aD and the fact that W 2 (y(s)) = 0. Hence, 

LT(S) = V 3 (X(s))X(s) + W 3 (y(S))y(S) = 

and U(s) = U (s) so that the definition (I3.28P is well-posed. Let us prove that x\ is Lipschitz. 
We have 

x\{X{s)) — x\{X{s)) = x(s) — x(s) 

x(s) ds 

v 2 {x{s))x{ s ) + m(y(s))y(s) ds 
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= 2 I V 2 {X{s))X{s)ds (by $$M ) 

J s 

< ||V 2 || £ oo \X{s)-X(s)\. 
Hence, X\ is Lipschitz. One proves in the same way that x 2 is Lipschitz. Since 

o < Vi(x(s))x( s ) < Vi(x(s))x( s ) +m{y(s))y( s ) = j( s ) 

the function Vi(X(s))X(s) belongs to L (R). Assume that there exists an s < s such that 
X(s) = X(s). Since X is increasing, it implies that X(s) = for all s £ [s,s] and therefore 
J^ oo V4(X(s))X(s)ds = $l x Vt{X(B))X{a)da and the definition (ET291) of Ji is well-posed. 
The same results hold for J 2 . Let us prove that U\ is absolutely continuous on any compact 
set. We consider X\ < ■ ■ ■ < Xjy and Sj such X(si) = Xj. We have 



N N 

Ei 

i=l i=l 



\Ui(X i+ i) - U!(Xi)\ = \Ui(s i+1 ) - I7i( ai )| 



Ui(si,s i+1 ) 



ds 



U x {s) 

< [ (v 3 (x)x + w 3 (y)y)ds 



<\\V 3 \\ Lao Xds 

JLIi(si,Si +1 ) 



+ meas(U i ( Si , Sm )) 1 /2( / W 3 (y) 2 y 2 dsf' 2 . 

JUi(si,s i+1 ) 

By we get W| < 2k ||W 4 || L oc (r) W 2 , and therefore W%(y)y 2 < CW 2 (^)3> 2 

CV 2 {X)Xy, by (pQ2j) . for some constant C. Hence, 



/ m(y) 2 y 2 ds<c f v 2 {x)xds 

JUi(s i ,s i+1 ) JUi(si,Si +1 ) 

<C Xds = Cmeas(Ui(Xi,X i+1 )) 



for some constant C . Finally, 

JV 

J2\Ui(X i+ i) - U^XJl < CCmeasCUiC^^Xi+i)) +meas(l> i (X i ,X i+1 )) 1 / 2 ) 

8=1 

and f/i is absolutely continuous. After differentiating (|3.28l) . we get 

u{(x)x = v 3 (x)x + w 3 (y)y 

and, after taking the square of this expression, we obtain 

(3.33) U[(X) 2 X 2 < 2(V 3 {X) 2 X 2 + W 3 {y) 2 y 2 ). 
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cW 2 W 4 j 2 



c v 2 xw A y (by 



Since 

W3(^) 2 ^ 2 

we have that (I3.33P implies 

u[{x) 2 x < c(v 3 (x) 2 x + Wi(y)y) 

< C(V 3 (X) 2 X + j) 

and, after a change of variables, we obtain 

\\U{\\ 2 L2 <C(||V 3 || 2 : 2 + || J\\ Lao )<oo. 

Hence, U[ belongs to 1? . Similarly one proves that U 2 is absolutely continuous on any 
compact set and U' 2 G L 2 (R). To prove that the property (I3.17P is fulfilled, we use Lemma 
EJand the fact that 1/(V 2 + V 4 ), 1/(W 2 + W 4 ) G L°°(R). The other properties of .F that 
^ has to fulfill can be checked more or less directly from the definition of Q$. □ 

The sets T and Go are not in bijection; otherwise we would not have introduced and 
indeed one can show that CoD / Idg . However, we have D o C = Idj?r, as we will see in 
Lemma 15.31 

Now we define how, from any initial data in V, that is, in the set of original coordinates, 
we define the corresponding element in J- '. 

Definition 3.8. We define the mapping L: V — > T where, for any (u, R, S, fi^) £ V, 
ip = (ipijipz) = L(n, R, S, fi, v) is defined as follows. We set 



(3.34a) 

(3.34b) 

and 

(3.34c) 

and 

(3.34d) 

and 

(3.34e) 

and 

(3.34f) 



x\{X) = sup{x G 
X2(Y) = sup{x G 



x' + cxd, x') < X for all x' < x}, 
x' + u(y — oo, x') < Y for all x' < x} 



J 1 (X)=X - Xl (X), J 2 (Y)=Y -x 2 (Y) 
Ut(X) = u( Xl (X)), U 2 {Y) = u(x 2 {Y)) 



Vi(X) 



Ki(X) 



R 



2c([/i 



(x x (X)K(X), V 2 (Y) 



s 



x 



J[(X) 



dX, K 2 (Y) 



2c(U 2 
J' 2 (Y) 



(x 2 (Y))x' 2 (Y) 



dY. 



,c(Ui(X)) ~ ' zy ' J-ooc{U 2 (Y)) 

Before proving the well-posedness of this definition, we check that we end up with the 
same initial data that was obtained at the end of Section O where fi and v were assumed 
to be absolutely continuous with respect to the Lebesgue measure. Then, the functions 



/i( — oo, x) 



1 



R 2 dx and v{— oo, x 



1 



S 2 dx 
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are continuous, and furthermore, (I3,34al) and (13,34b|) rewrite as 

i-x 1 (X) i fX2(Y) i 

xi{X) + -R 2 dx = X and x 2 (Y) + / -S 2 dx = Y. 

J —OO J —OO ^ 

We sum these two equalities, and, since x\{X(s)) = X2(y(s)) = x(s) and X + y = 2s, we 
get 

f-x(s) 1 

2x{s)+ / -{R 2 + S 2 )dx = 2s 

J —oo * 

and recover (|2.27j) . In the definitions (I3.34p , we use the degree of freedom we have in the 
new set of coordinates to set the values of x\ and x 2 in such a way that their derivatives, 
x[ and x' 2 , are bounded. 

Proof of well-posedness of Definition Clearly, the definition of x\ yields an increasing 
function and limx^±oo x\ (X) = ±oo. For any z > x±(X), we have X < z + ^((—00,2:)). 
Hence, X — z < and, since we can choose z arbitrarily close to x\{X), we get X — 

x\{X) < /i(IR). It is not hard to check that x\{X) < X. Hence, 

(3.35) \xi(X) - X\ < n(R) 

and \\xi — Id|| ioo < /x(R). Let us prove that x\ is Lipschitz with Lipschitz constant at 
most one. We consider X, X' in R such that X < X' and x\{X) < x\{X'), It follows 
from the definition that there exists an increasing sequence, z^, and a decreasing one, 
Zi, such that lim^oo = x\(X), lim^oo z\ = x±(X') with [/,((— 00, z[)) + z[ < X' and 
fj>((— 00, Zi)) + Zi > X. Combining the these two inequalities, we obtain 

(3.36) m((-°°> z 'i)) ~ m((-o°> z i)) + A - Zi < X' - X. 

For j large enough, since by assumption x\{X) < x\{X'), we have z; t < z\ and therefore 
fj,((— 00, z'j)) — fi{(— 00, Zi)) = n([zi, z^)) > 0. Hence, z[ — z% < X' — X. Letting i tend to 
infinity, we get x\{X') — x±(X) < X' — X . Hence, x\ is Lipschitz with Lipschitz constant 
bounded by one and, by Rademacher's theorem, differentiable almost everywhere. Following 
[6], we decompose [i into its absolute continuous, singular continuous and singular part, 
denoted fi ac , [i sc and fi s , respectively. We have ^ ac = \R 2 dx. The support of /i s consists 
of a countable set of points. Let H(x) = fx((—oo,x)), then H is lower semi-continuous and 
its points of discontinuity exactly coincide with the support of /i s (see [6]). Let A denote 
the complement of x^ 1 (supp(/i s )). We claim that for any X £ A, we have 

(3.37) n((-oo, Xl (X))) + xi(X)=X. 

From the definition of X\(X) follows the existence of an increasing sequence Zi which 
converges to x\{X) and such that H(zi) + z% < X. Since H is lower semi-continuous, 
rmij_ >00 H(zi) = H(x\{X)) and therefore 

(3.38) H(xi(X)) + xi{X) < X. 

Let us assume that H(x\(X)) + x\(X) < X. Since x±(X) is a point of continuity of H, we 
can then find an x such that x > x\{X) and H(x) + x < X. This contradicts the definition 
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of x\{X) and proves our claim (|3.37|) . In order to check that (I3.16P is satisfied, we have to 
compute the derivative of x\. We define the set B\ as 

B x = \x e R I lim — u((x - p, x + p)) = \r 2 (x) 1 . 
[ pio 2p yy 4 J 

Since |i? 2 (x) dx is the absolutely continuous part of p, we have, from Besicovitch's deriva- 
tion theorem (see P), that meas(fif) = 0. Given X £ srH^i). 

we denote x = x\(X). We 

claim that for all i G N, there exists < p < A such that x — p and x + p both belong to 
supp(^ s ) c . Assume namely the opposite. Then for any z G (x — j,x + A) \ supp(// s ), we 
have that z' = 2x — z belongs to supp(/i B ). Thus we can construct an injection between the 
uncountable set (x— A, x + A) \supp(/i s ) and the countable set supp(^ s ). This is impossible, 
and our claim is proved. Hence, since x\ is surjective, we can find two sequences Xi and 
X[ in A such that \{x\{Xi) + xx(X-)) = x\{X) and x x {X'^ - x x {Xi) < A. We have, by 
(|3.37|) . since x\{Xi) and xi{X[) belong to A, 

(3.39) p([x 1 (X i ),x 1 (X' i ))) + Xl (Xl) - x 1 (X i ) = X'i-Xt. 

Since xi( Xj) £ supp(^ s ), we infer that //({xi(Xj)}) = and p,([xi(Xi), xi(X-))) = p({x\{X i ) 1 x l {X! j ))). 
Dividing (13, 39p by X[ — X,i and letting i tend to oo, we obtain 

(3.40) x' x {X)^R 2 {xi{X)) + x[(X) = 1 

where x\ is differentiable in xj~ 1 (i?i), that is, almost everywhere in x^ l {B\). We will use 
several times this short lemma whose proof can be found in [9]. 

Lemma 3.9. Given an increasing Lipschitz function f : M. —> R, for any set B of measure 
zero, we have f' = almost everywhere in f^ 1 (B). 

We apply Lemma I3U1 to B\ and get, since meas(-BJ) = 0, that x[ = almost everywhere 
on x^ 1 (Bl). On x^ 1 (B\), we proved that x[ satisfies (|3.40p . It follows that < x[ < 1 
almost everywhere, which implies, since J[ = 1 — x[, that J[ > 0. Prom (|3.40l) . we get 

xx(Xyj[(X) = x'^Xf-^xriX)) = (ciU^X^X)) 2 . 

Let us prove that U\ is absolutely continuous on any bounded interval. We consider a 
partition X\ < ■ ■ ■ < Xjy. We have 

N 

J^pxiXi+^-U^Xi)] < / \u x \ dx. 

Given M > 0, for any e > 0, there exists 5 such that for any set A C [—M,M], we have 
that meas(j4) < 5 implies f A \u x \ dx < e, because u x S L\ oc . We have 

meas(U? =1 (xi(Xi),xi(X i+ i))) < ||a;'i|| Loo meas{\jf =1 {Xi, X i+1 )) 

and since x' x G L°°, it follows that for any partition such that \X-i — Xi + \\ < 5, we 
have YaLi I U\(Xi + i) — Ui(Xi)\ < e and U\ is absolutely continuous in any compact. Let 
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us consider a curve (X,y) G C such that a;i(A'(s)) = X2(y(s)) (such curves exist, see 
Definition 13.5(1 , We differentiate Z7i(<Y) and obtain 

U[(X)X = u x ( Xl {X))x' x {X)X = 
Here we have use the fact that 

x[(x)x = x ' 2 (y)y 

which follows from x\{X) = xzty). From (|3.35|) . we obtain || Ji|| ioo < /u(R). We have, since 
J[ > 0, 

I Kill £,2 < IrilLoo H^ilLi ^ IKiIIl 00 - 

After a change of variables, we get 

.2 



/ yi(X) 2 dX < ^- / R 2 (x)dx < 

JR 4 J R 



and 

7V2/v\JV_ / 2/ _/2 



C7f(X)dX= / u%{x 1 {X))x'{{X)dX 
Jr 

< [ ul(x 1 {X))xf l {X)dX = I u 2 x (x)dx <oo, 



so that both V\ and U[ belong to L 2 . Similarly, one proves that V 2 and U 2 belong to L 2 . 
Let B 3 = {£ G E | x'j < \}. Since x' x - 1 > 0, B 3 = G R \ \x[ - 1| > |}, and, after using 
the Chebychev inequality, as — 1 = — J[ E I 2 , we obtain meas(i?3) < oo. Hence, 

/ Uf(X)dX = [ Ul{X)dX+ [ Uf(X)dX 
Jr Jb 3 Jbc 

< meas(i?3) Unll^oo + 2 / (uoxi) 2 x'idX 

Jb§ 

< meas(S 3 ) [|it||ioc + 2 \\u\\ 2 L2 , 

after a change of variables, and U\ G L 2 . Similarly, one proves that U 2 £ L 2 . □ 

In this section we have shown how to construct, from a given initial data in V, an 
element in J- (via the mapping L) and then, from an element in and element in Qq 
(via the mapping C). From an element in Qq, we can finally construct the corresponding 
solution of (12331 . 

Now, we turn to the existence of solution to (I2.13P for given data in Q. 
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4. Existence of solution for the equivalent system 

4.1. Short-range existence. We first establish the short-range existence of solutions to 
(|2.13l) . The difficulty here consists of taking into account initial data defined on a curve 
which may be parallel to the characteristic curves X = constant or Y = constant. In the 
following, we will denote by any rectangular domain of the type 

n=[X l ,X r ] x [Y t ,Y r ], 

and we denote si = \{Xi + Y{) and s r = \{X r + Y r ). We define curves in £1 as follows. 

Definition 4.1. Given = [Xi, X r ] x [Y], Y r ], we denote by C(O) the set of curves in ft given 
by (X(s),y(s)) for s G [sj,s r ] which match the diagonal points ofVL, that is, X(si) = X\, 
X(s r ) = X r , y(si) = Yi, y(s r ) = Y T , and such that 

(4.1a) X - Id, y - Id € W 1 ' 00 ^, s r }), 

(4.ib) x > o, y > 0, 

(4.1c) \^ X ^ + y ( s ^ = s ' f° r aU seM ~ 

We set 

y)\\ C (n) = u* - ^\\L^([ Sl ,s r ]) + \\y - l d\\ L oo {[ShSr]) . 

In this subsection we will construct solutions on small rectangular domains Q. We intro- 
duce the set £/(Q) which is the counterpart of Q on bounded intervals. Elements of G(Q) 
correspond to a curve in C(Q) and data on this curve. 

Definition 4.2. Given a rectangular domain O = [Xi, X r ] x [Yi, Y r ), let = (X, y, Z, V, W) 
where (X,y) £ C(Q) and Z(s),V(X),W(Y) are three five- dimensional vector-valued mea- 
surable functions. Using the same notation as in (13. 3p . we set 

\\ @ \\g(n) = ll^llx 2 ^,^]) + \\ va \\^{[x u x r ]) + H^ll^a^y,.]) 
where we denote U = Z3 and 

1 



l®llb(c) - ll(*>^)|| c( n) + || _2 \\L°°([x„x r ]) + W w „1 Wa ll £00 (K. y r]) 



+ \\ za \\L^([ Sl<Sr }) + ll^ a ll£°°(pQ,x r ]) + W wa \\L^([Y t ,Yr]) • 
The element Q belongs to Q(fl) if the following four conditions hold: 
(i) 

lll®lllg(Q) < 00 

and therefore \\@\\g(Q\ < 00 because we here consider a bounded domain. 

(") 

(4.2) V 2 ,W2,^4,V4,W4 > 0. 

(iii) For almost every s, we have 

(4.3) Z{s) = V{X{s))X{s) + w{y{s))y{s). 



28 HOLDEN AND RAYNAUD 

(iv) For almost every X and Y , we have 
(4.4a) 2V 4 {X)V 2 {X) = (c(U)V 3 (X)) 2 , 2w 4 {y)w 2 {y) = (c(U)W 3 (y)) 2 , 

(4.4b) v 2 {x) = c(u)Vi(x), w 2 {y) = -c(u)m(y), 

(4.4c) V 4 {X) = c(U)V 5 (X), Wt(y) = -c(U)W 5 (y). 



We introduce the Banach spaces W^°°(Q) and Wy 00 ^) defined as 
(4.5) 

W%°°(n) = L 0O (\Y l ,Y r ],W l > 0O ([X l ,X r ])) 1 W Y >°°(n) = L~([X l ,X r ],W 1 >°°(\Y l ,Y r ])), 

and the Banach spaces L^(Q) and Ly ($7) defined as 

L%(Q) =L°°([Y l ,Y r ],C([X l ,X r })), L?(Q) = L^{[X h X r lC{[Y h Y r ])). 

Let us consider (X,y,Z,V,W) E G{&)- By definition, the functions X and y are increas- 
ing. To any increasing function, one can associate its generalized inverse, a concept which 
is exposed for example in Brenier [2J. More generally, an increasing function (not neces- 
sarily continuous) can be identified as the subdifferential of a convex function (which is a 
multivalued function). The generalized inverse is then the subdifferential of the conjugate 
of this convex function. We do not use this framework here and prove directly the results 
we need. 

Definition 4.3. Given Q = [Xi,X r ] x and a curve (X,y) G C(Q), we define the 

generalized inverse of X and y, respectively, as 

(4.6) a(X) = sup{s E [s h s r ] | X(s) < X} for X E {X h X r ], 

(4.7) f3(Y) = sup{s E [si, s r ] | y(s) < Y} for Y E {Y\, Y r \. 

We denote X~ x = a and 3 7 " 1 = f3. 

The generalized inverse functions X~ l and y~ l enjoy the following properties. 

Lemma 4.4. The functions X~ l and y~ l are lower semicontinuous nondecreasing func- 
tions. We have 

(4.8) X o X^ 1 = Id and y o y- 1 = Id, 
and 

(4.9) X' 1 o X(s) = s for any s such that X(s) > 0, 

(4.10) y^ 1 o y( s ) = s for any s such that y(s) > 0. 
Definition 14.31 extends naturally to curves in C and Lemma l4~4l still holds. 

Proof. We prove the lemma only for X -1 , as the results for y~ l can be proved in the same 
way. Let us prove that a is nondecreasing. For any X < X, there exists a sequence Sj 
such that linij^oo s» = a(X) and X(si) < X. Hence, X(si) < X which implies Sj < a(X), 
which after letting i tend to infinity, gives a(X) < a(X). Let us prove that a is lower 
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semicontinuous. Given a sequence Xi such that lim^oo Xi = X, for any e > 0, there exists 
s G [si, s r ] such that 

(4.11) a(X) > s > a(X) - e 

because a(X) > s/ for all X G (Xi, X r \. It implies X(s) < X as, otherwise, X < X(s) would 
yield a(X) < s, which contradicts (|4.11|) . Thus, for large enough i, we have X(s) < Xi so 
that s < a(Xi). Combined with (14. lip , it implies 

a(X) — e < s < liminf a(Xi) 

and, as e is arbitrary, we get that a is lower semicontinuous. Let us prove (|4.8p . Given 
X G (Xi,X r ], we consider an increasing sequence Si such that lim^oo Si = a(X) and 
X(s.j) < X. Letting i tend to infinity, since X is continuous, we get X(a{X)) < X. Assume 
that X(a(X)) < X, since X is continuous, there exists s such that X(a(X)) < X(s) and 
X(s) < X. The latter inequality implies that s < a(X) which, by the monotonicity of 
X, yields X(s) < X(a(X)) and we obtain a contradiction. Let us prove (|4.9p . We denote 
N = {s £ [si, s r ] | X(s) > 0}. We consider a fixed element so £ M. We have 

(4.12) aoX(s )<s . 

Indeed, by the monotonicity of X, for any s £ {s £ [si,s r ] \ X{s) < X(sq)} 7 we have 
s < sq and therefore, after taking the supremum, we obtain (I4.12p . Let us assume that 
a o X(sq) < so- We denote si = a o X(so). By (14. 8p . X(s±) = X(so), and from the 
monotonicity of X , it follows that X(s) = X(sq) = X(si) for all s G [so, si]- It implies that 
X(so) = 0, which contradicts the fact that so G M. □ 

In the following and when there is no ambiguity, we will slightly abuse the notation 
and denote 3^ o X~ 1 {X) and X o y~ l (Y) by y(X) and X(Y), respectively. The curve 
(X(s),y(s)) is almost a graph as it consists of the union of the graphs of the functions 
X i ► y(X) and (after rotating the axes by ^) Y i— > ^(l^). We prove the existence of 
solutions to (|2.13|) on rectangular boxes. First we give the definition of solutions. 

Definition 4.5. We say that Z is solution to (I2.13P in ft = [Xi,X r ] x if 

(i) we have 

z g w 1 ' 00 ^), z x g w^in), z Y g iyi ,00 (^); 

(ii) and for almost every X G [X;,X r ], 

(4.13) (Z X (X, Y)) Y = F(Z)(Z X , Z Y )(X, Y); 
and, for almost every Y G [Y^,!^], 

(4.14) (Z Y (X,Y)) X = F(Z)(Z X ,Z Y )(X,Y). 

We say that Z is a global solution to (|2.13p if Z is a solution to (I2.13P as defined above, for 
any rectangular domain 17. 

The regularity that we impose is also necessary to extract the relevant data on a curve 
from a function defined in the plane, as it is explained in the following lemma. 
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Lemma 4.6 (Extraction of data from a curve). We consider a five- dimensional vector 
function Z in such that 

Z g W 1 ' 00 ^), Z x g Wp°°(Q), Z Y g W]i°°(p) 

for any rectangular domain VL. Then, given a curve (X,y) G C, let (Z,V,W) be defined as 

(4.15) Z(s) = Z(X(s),y(s)) for alls£R 

and 

(4.16a) V(X) = Z x (X,y(X)) for a.e. lEl, 

(4.16b) W{Y) = Z Y {X(Y),Y) for a.e. Fel, 

or, equivalently, 

V(X(s)) = Z x (X(s),y(s)) for a.e. s G M such that X(s) > 0, 
W(y(s)) = Z Y (X(s),y(s)) for a.e. s£R such thaty(s) > 0. 

We have (Z,V,W) G i~(R) and we denote 9 = (X,y,Z,V,W) by 

Zm{X,y). 



Proof. We consider a domain Q = [Xi,X r ] x [Yj,y r ]. We claim that for any / G Wy°°(fi) 
then f(X) = f(X,y(X)) is measurable and f(X,y(X)) G L°°{[X h X r }). It suffices to show 
that the linear mapping f t-> f from Wp°°(n) to L°°([X h X r ]) is well-defined on simple 
functions and continuous. We assume that / is a simple function, that is, 

N 

f(X,Y) = J29j(y)XA 3 (X) 

3=1 

where xa denotes the indicator function of the set A, Aj are disjoint measurable sets and 
9j G W^ilYuYr]). Then, f(X) = Ef=i 9j(y(X)) XAj (X) is measurable (as X ^ y(X) 
is lower semicontinuous) and 



f(X) < max esssup \gj(y(X))\ 
je{i,-,N} xe[x lr x r ] 



esssup 

X£[Xi,Xr 



so that / G L°°([Xi,X r ]). Note that we need gj G W^ 00 [\Y u Y r ]) as, if gj only belongs to 
L°°([Yi,Y r ]), we do not have in general esssup Xe[Xi Xr ] \gj(y(X))\ < \\9j\\ LaC ([ Yl Xr]) as the 
function X y(X) may send sets of strictly positive measure to a set of measure zero (for 
example if y is constant on an interval). Therefore the continuity in the Y direction which 
is necessary to make meaning of (I4.16p . Using the same type of estimate, one gets that 



'([x,,x r ]) ^ M\\w*-°°(n) 
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which concludes the proof of the claim. Similarly one proves that, for any / £ W^°°(Q), 
the mapping Y i— ► f(X(Y),Y) is measurable and belongs to L°°([Xi,X r ]). Hence, we get 
that (Z,V,W) £ L™ C (R), ' □ 

The decay of Z at infinity in the diagonal direction is more conveniently expressed in 
term of the function Z a which we now define as 

(4.17) Z 2 a = Z 2 - -{X + Y) and Zf = Z { for % £ {1, 3, 4, 5}. 

Even if we are not concerned yet with the behavior at infinity, it is convenient to introduce 
Z a already here to write the estimate in a convenient way. We now introduce the set 7i(Q) 
of all solutions to (I2.13P on rectangular domains, which satisfy additional properties. 

Definition 4.7. Given a rectangular domain ft = [Xi,X r ] x [YJ,Y r ], let be the set 

of all functions Z which are solutions to (I2.13|) in the sense of Definition \4-5\ and which 
satisfy the following properties 



(4.18a) x x = c(U)t x , x Y = -c(U)t Y , 

(4.18b) J x = c(U)K x , Jy = -c(U)K Y , 

(4.18c) 2J x x x = (c(U)Ux) 2 , 2J Y x Y = (c(U)U Y f , 

(4.18d) x x > 0, J x > 0, 

(4.18e) x Y > 0, J Y > 0, 

(4.18f) x x + Jx > 0, x Y + J Y > 0. 



We have the following short-range existence theorem. 

Theorem 4.8. There exists an increasing function C such that, for any ft = [X[,X r ] x 
[Y h Y r ] and 9 = (X,y,Z,V,W) in Q{9), if s r - si < 1/C(|||0|||g(n)), then there exists a 
unique solution Z £ TL{Q) such that 

(4.19) B = Z*(X,y). 

Proof. We use a Picard fixed-point argument. Define B as the set of elements (Zh, Z v , V, W) 
such that 

Z h e[L^f, Z v e[m\ V£[L?] 5 , We[L%] B 

and 

5 

(4.20) + II^IL- + H^Hif + H^ll^) < 2|||e||| e( n) 
i=i v Y 

where we use for Z^ and Z v the same notation given in (|4.17l) for Z. For the fixed point, the 
functions Z^ and Z v coincide and are equal to the solution Z, see below, but it is convenient 
to define both quantities in this proof and keep the symmetry of the problem with respect 
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to the X and F variables. We introduce the mapping V given, for any (Z^, Z v , V, W) G £>, 
by V(Z h , Z v , V, W) = (Z h , Z v , V, W) where 

r-x 



(4.21a) Z h (X,Y) = Z{y-\Y))+ [ V(X,Y)dX 

JX(Y) 

for a.e. Y G [Y u Y r ] and all X G pT/,X r ], 

(4.21b) Z„(X,F) = + / W(X,Y)dY 

Jy(x) 

for a.e. X G [Xj,X r ] and all Y e [Y h Y r ], 

(4.21c) K(X,F)=V(X) + / Y F(Z h )(V,W)(X,Y)dY 

Jy{x) 

for a.e. X G [X h X r ] and all F G [Y h Y r ], 

(4.21d) W(X,F) = W(F)+ / F(Z h )(V;W)(X,y)dX 

JX{Y) 

for a.e. F G [F/,F r ] and all X G pQ,X r ]. Let us consider a solution Z to (I2.13P which 
satisfies (I4.19p . For any (X, F) G $7, we have 

Z(X,y(8))=Z( 8 )+ / Z x (X,^(s))dX 

which, after taking s = y~ l (Y), yields 

r x 

Z(X,Y)=Z(y- 1 (Y))+ Z x (X,Y)dX, 

Jx(Y) 

by (14. 8p . Similarly, one proves that 

Z(X,F) =Z( < %- 1 (X))+ / Zy(X,F)dF. 

Jy(x) 

For every almost every X G [^,X r ] and all F G [F;,F], we have 

Z X (X,Y) =Z x (X,y(X))+ [ Y F(Z)(Z x ,Z Y )(X,Y)dY, 

Jy(x) 

and therefore, by (14.19p . we get 

(4.22) Z X (X,Y) = V(X)+ [ Y F(Z)(Z x ,Z Y )(X,Y)dY. 

Jy{x) 

Similarly, we have 

(4.23) Z Y (X,Y) = W(Y)+ [ X F(Z)(Z x ,Z Y )(X,Y)dX 
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for all X G [Xi,X r ] and a.e. Y G [Yi,Y r ]. Thus, if Z is a solution to (12 . 13|) (in the sense 
of Definition I4.5P which satisfies (|4.19l) then (Z,Z,Zx,Zy) is a fixed point of V. Since 
< X < 2 and < Y < 2, we have 

(4.24) \X - X(Y)\ = \X(a{X)) - X(a(Y))\ < X(s r ) - *(*,) < 2( Sl - s r ) < 25 
and, similarly, 

(4.25) \Y -y(X)\ <2{si- s r ) <25. 

We can choose 5 small enough, depending only on |||0|||g(n), such that the mapping V maps 
B into B. Let us check this in more details only for the second component of Zh. We have 

Z a K2 = Z ha -\(X + Y) = Z 2 {y-\Y)) -\(X + Y) + [ X V(X,Y)dX 
2 2 Jx(Y) 

and, denoting s = y^ 1 (Y), 

Z 2 (s) -l(X + Y) = Z 2 (s) ~\{X + y{s)) = Z 2 (s) - l -{X{s) + y{s)) + ~(X - X{s)) 

= Z^s) + l -{X-X(Y)). 

Hence, 

IfelUtn) ^ W\ z 2\hm + 5(1 + 2|||e||| s(n) ) 
by (I4.24p and (|4.20|) . After doing the same for the other components, we get 

5 

DK*IIl ? + \\z°A\ Lf + \\Ml S + II^IM < wmkm + SC 

1=1 

for a constant C which depends only on |||0|||g(n). Hence, by taking 5 small enough, the 
mapping V maps B into B. Using the fact that F is locally Lipschitz (because it is bi- linear 
with respect to the two last variables and depends smoothly on U = Z3), we prove that V 
is contractive. Hence, V admits a unique fixed point that we denote (Zh, Z v , V, W). Let us 
prove that Zh = Z v . It basically follows from the fact that Wx = Vy. Let us now denote 
by J\f x the set of points X G [X h X r ] for which (|4.21bh and (|4.21cD hold (by definition, 
the set Mx has full measure). Similarly we denote by My the set of points Y G for 
which (I4.21ah and (I4.21dll hold. We have meas([Xi, X r ]\J\f x ) = meas([Y/, Y r ] \Afy) = so 
that meas(fi \flfx x Ny) = 0. For any (X, Y) G Mx x My, we have 

Z h (X,Y) - Z V (X,Y) = Z(y-\Y)) + / V(X,Y)dX 

Jx(Y) 

r Y . . 

(4.26) -Z(X~ 1 (X))- W(X,Y)dY. 



iy{x) 

Since the terms involving F cancel, we obtain by t|4.21c|) and (|4.21dj) that 

r-X pY pX i-Y 

/ V(X,Y)dX- W(X,Y)dY= V(X)dX- W(Y)dY. 

Jx(y) Jy(x) Jx(Y) Jy{x) 
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Returning to the rigorous notation, we get 
(4.27) 

/ v(x,Y)dx- / w(x,Y)dY= / v(x)dx- / w(Y)dY 

Jx(Y) Jy(x) Jx{y-\Y)) Jy(x-^{x)) 



where we have also used (|4.8p . We proceed with a change of variables in the two integrals 
on the right-hand side of (I4.27|) and get 

,x{x-\x)) ryiy-HY)) 



/ V(X)dX- / W(Y)dY 

Jx(y-i(Y)) Jy(x-Hx)) 

ry-HY) , 

W(y{s))y{s) +V(X(a))X{s) \ ds 



'X- 1 (X) 

t-y-HY) . 

= - / Z{s)ds by (1431) 

Jx-^ix) 

(4.28) =Z(X-\X))-Z(y- 1 {Y)) 

and combining (l4~26i (14371) and (I43H1) . we get that Z ft (X, Y) = Z„(X, Y) for all (X, Y)_G 
A/"x X A/y, that is, almost everywhere. We denote Z = Zh = Z v . For any (X, Y") and (X, Y) 
belonging to Mx x A/y, we get, by using (|4.21a|l and (|4.21bj) . that 

Z(X,Y) - Z(X,Y) = [ V(X,Y)dX+ f W(X,Y)dY. 
Jx Jy 

Hence, by using the bound (|4.20l) . 

\Z(X,Y) - Z(X,Y)\ < 2|||e||| e(n) (|X - X\ + |Y - Y|) 

and Z is Lipschitz in Mx x My- It implies that Z is uniformly continuous in Mx x My 
and there exists a unique continuous extension of Z to the closure of Mx x My, that is, Q. 
From (j4.21aj) and rt4.21bj) , we get that 

(4.29) Z X (X,Y)=V(X,Y) 
for all Y G [Y/, Y r ] and a.e. X G [X,X r ] and 

(4.30) Z Y {X,Y) = W(X,Y) 

for all X G [X/,X r ] and a.e. Y G [Y,Y r ]. By using the fact that (Z, Z, Zx > Zy) is a fixed 
point in B and (143911 and (14301) . we can check that Z x G lYy'°°(0) and Zy G w£°°(fi). 
By density, we can prove that 

(4.31) Z(X,Y)- Z(X,Y)= [ Z x (X,Y)dX, 

Jx 

not only for almost every Y G [Y/,Y.] as (I4.21al) yields, but for all Y G [Y,Y r ]. Indeed, 
for any Y G [Y,Y r ], there exists a sequence Y n G My such that lin^^oo Y n = Y as 
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meas([l/,y r ] \A/y) = and we have 

(4.32) Z(X,Y n )-Z(X,Y n ) = f Z X {X, Y n ) dX. 

Jx 

Since Z x G W Y '°°(Ci), we have that \\Z X ( • , ^)||l°°([Xi,X,.]) - ll^xllv^.^n) - 2 lll llb(O) 
and, for a.e. X G pQ,X r ], lim n _>oo Zx(X, Y n ) = Zx(X,Y) for almost every X. Hence, by 
Lebesgue dominated convergence theorem and the continuity of Z, (I4.32P implies (|4.31|) . 
It remains to check that Z satisfies (|4.19l) . Since (Z, Z, Zx , Zy) is a fixed point of V, we 
have, by (l4T21cl) and (I4.21djl . that 

Z x (X,y(X)) = V{X) and Z Y (X(Y),Y) = W{Y) 

for a.e. X G pQ,X r ] and F G [l;,y r ], respectively. It remains to check that Z satisfies 
(|4.15ll . On one hand, we have that 

(4.33) Z(X(s),y(s)) = Z{X-\X{s))) 

by (|4.21hj) and, by (l43]l . it implies (l4~15ll for all s G [s z ,s r ] such that X{s) > 0. On the 
other hand we have 

(4-34) z{x(s),y{s)) = z{y-\y( s ))) 

by (I4.21b[) and, by ([4TTuD . it implies (05jl for all s G [s h s r ] such that > 0. Since 

X + y = 2, the set of all s G [s;,s r ] such that X(s) > or y(s) > has full measure 
and therefore, for almost every s G <|4.15jl holds. By continuity, we infer that (|4.15p 

holds for all s G [aj, s r ]. Hence, we have proved that Z is a solution to (|2.13|) which satisfies 
(|4.19l) if and only if it is a fixed point of V. Since the fixed point exists and is unique, 
we have proved the existence and uniqueness of the solution. Let us define the functions 
v G W Y '°°(n) and w G in as 

v = xx — c(U)tx and w = xy + c(U)ty. 

We want to prove that t; and tu are both zero. After some computations using the governing 
equations (I2.13p . we obtain 

V y = X X Y ~ c'(U)Uyt X ~ c(U)t X Y 

(4.35) = {Uyv + U x w) 



and 



(4.36) = py^y- (Uyv + U x w 



2c(U) 

Wx = X X Y + c'(U)Uxty - c(U)t X Y 
2c(U) 



It follows that 



v(X,Y)= f Y f£l(Uyv + U x w) 
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r x c'(U) 
Jx{x) 2c {U) 

which implies, after using (|4~24l) . ([4251) and (1420) . that 

c'(U) 



< 



r, (TT\ (IMIl?? \\ W \\L 



r - 2c(U) ^ " ,,JJ * 

< 52|||e||| g(n) max|^(||u|| i? + ||u>|| L? ). 

Hence, by taking 5 smaller if necessary, we get IMI^oo = ||io|| i00 = 0. Let us now introduce 
z G Wy'°°(n) asz = 2J x x x - (c{U)U x ) 2 . We have 

z Y = 2J XY x x + 2J x x XY - 2c(U) 2 U x U XY - 2c(U)c'(U)U Y {U x f 

and z(X,y(X)) = for X G [Xi,X r ], the unique solution to (|4.37j) is z = 0. One proves 
in the same way that J Y x Y = 2(c(U)U Y ) . Let us now prove (|4.18f|) . Since the initial 
data belongs to G(Sl), we have \\l/(x x + J X )(X, y(X))\\ Lao ^ Xl ,x r ]) < lll@lllg(n) as l/(x x + 
J x )(X,y(X)) = 1/(V 2 + W 2 )(X) for a.e. X G [X h X r \. For all fixed X such that l/(x x + 
Jx)(-X", lypO) < that is f° r almost every X £ pQ,X r ], we define 

F* = inf{y G [^,Y r ] I y < and (x x + Jx)(X,Y') > for all Y' > y} 

and similarly 

Y* = sup{y G [Yi , y] | y > ypQ and (xx + > for all Y' < Y}. 

On (y,y*), we have (xx + J X )(X,Y) > and we define 

q{Y) = (x x + J x )(X,Y)- 
Let us assume that y < Y r and therefore, by continuity, 

(4.38) (x x + J X )(X,Y*) = 0. 

On (y,y*), we have <?(y) > and, since Jxxx > by (|4. 18c|) . it implies that 

(4.39) > and J x > 0. 
By using (|2.13l) . we obtain 

xxy + Jxy 

q Y = 

(4.40) 

From (I4.18cp . we infer that 



(xx + Jx) 2 

c!{U) U Y (x x + Jx) + (Jy + x y )Ux 
'2c{U) (xx + Jx) 2 



lUxl = ^^ V ^-2~^ {JX + XX) - 
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Hence, (|4.40l) yields 

\c'(U)\ 

qy < VtTttCI^I + l J ^l + l Xy D9 ^ C ^ 
2c(f7) 

for some constant C which depends only 

(4.41) esssup (|*7| (X,Y) + |Z y | (X, Y)), 

^G[Vi,V r ] 

that is |||@|||g(T2), by (I4.20p . By Gronwall's lemma, it follows that q cannot blow up in Y*, 
contradicting (I4.38P and 

(4.42) ^TT^ Y ) ^ V^V iX ^ ClY ~ yiX)l 

for a.e. X G pQ,X r ] and all Y G [Yi,Y r ]. In the same way, one proves that Y* = Y\, Hence, 
we have proved that, for almost every X G [Xi,X r ], 

xx {X, Y)>0 and J X (X, Y) > for all Y 
and ||l/(xx + Jx)(X, Y)\\ woo ^ Yi Yr \ is bounded, for almost every X £ [Xi,X r ], by a con- 
stant which is independent of X and therefore l/(xx + Jx) G Wy°°{&)- This concludes 
the proof of (|4.18d|) and the first identity (|4.18f|) while (14. 18e|) and the second identity in 
(|4.18fj) can be proven in the same way. □ 

4.2. A priori estimates. Given a positive constant L, we call domains of the type 

D = {(X, Y) G R 2 | \Y - X\ < 2L} 

for strip domain. Strip domains are correspond to domains where time is bounded. We have 
the following a priori estimates for the solution of (|2.13p . The energy J(X,Y) is bounded 
in the whole plane while Z a (that is, Z, up to a shift in the second component) and its 
derivatives are bounded in every strip domain. 

Lemma 4.9. Given Q = [X h X r ] x [Y h Y r ] and G = (X, y, Z, V, W) G G(n), let Z G H{9) 
be a solution to (12 . 13|) such that G = Z • (X,y), Let £q = ||^4||£,oo(r s s j) + ll^sllioo^g ar n. 
Then the following statements hold: 

(i) Global boundedness of the energy, more precisely, 

(4.43) < J(X, Y) < £ for all (X, Y) G fl and \\K\\ Loo(n) < (1 + k)£ 

where J = Z4 and K = Z5. 

(ii) The function Z and its derivatives remain uniformly bounded in strip domains. 
More precisely there exists a nondecreasing function C\ = C\(L, \\\@\\\g(n)), such 
that, for any L > and any X and Y such that \X — Y\ < 2L, we have 

(4.44) \Z a (X,Y)\ < Ci, \Z X (X,Y)\ < C u \Z Y (X,Y)\ < C x 
and 

(4.45) ] (X,Y)<C h — (X,Y)<d. 
XX + J X Xy + Jy 

The condition (ii) above is equivalent to the following condition (Hi): 
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(iii) For any curve (^,3^) S C(O), we have 
(4-46) %Z.{X,y)lg{Ci) <Ci 

where C\ = C\{\\{X , y)\\c(o.)-, Ill®lll<7(n)) is a given function which is increasing with respect 
to both its arguments. 

The inequalities in (Oil hold in fact in W Y '°°(n) and W x '°°(p,), respectively. 

The inequalities in (|4.45j) hold in Wy'°°(f2) and W x ,co (Q), respectively. 

Proof. Given P = (X,Y) G 0, let s = y- l {Y) and s l = X~ l (X). We denote P = 
(X(so), Y(so)) and P\ = (X(si), Y(si)). We assume that so < -si (the proof for the other 
case is very similar). Since X and y are positive, it implies that X = X(si) > X(so) and 
Y = Y(so) < 3^(-si). Then, because Jx > 0, Jy > and Z^ > 0, we have 

(4.47) < Z 4 (s ) = J(P ) < J(P) < J(Pi) = Z 4 (si) < £ 

which gives the first inequality in (14,43p . By (14,18b|) . we get 

\K(P)-K(P )\<k(J(P)-J(P )) 

which implies 

(4-48) \K(P)\ < \K(P )\ + k(J(P) - J(P )) 

and 

\K(P)\ < (1 + k)£ 

by (I4.43p . Since xx > 0, we have 

(4.49) x(P) > x(P ) = Z 2 (s ) > -|||e||b ( n) + s . 
Since \{X + Y) = Y + ±(X - Y) < Y(s ) + L, it follows that 

x(P) -l(X + Y)> -\\\@\\\g m +s - y(s ) — L> -2|||e||| g(n) - L. 
Similarly, using that xy > 0, we get 

(4.50) x(P) < x(Pt) = Z 2 ( Sl ) < \\\@\\\ g{ n) + si- 
and 

(4-51) x(P)-^X + Y)<2\\\@\\\ gm +L. 

Hence, \x(P) -\{X + Y)\< 2\\\Q\\\ g{n) + L. We have 



(4-52) \t(P) 



/ t Y (X,Y)dY < -^. dY <n{x{P 1 )-x{P)). 

Jy Jy c(U) 



Since 

(4.53) x(Pi) - x(P) < x(Pi) - x(P ) = Z 2 {s x ) - Z 2 (s ) < 2|||e||| e(n) + si - s 
and 

(4.54) 5l - So = Sl - - (s - X(s )) + Y-X< 2\\\G\\\g {n) + 2L, 
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it follows from (gSZjl that \t(P)\ < K(4|||e||| g(n ) + 2L). We have 

ry('i) 

\U(P)\ < \U(Pi)\ + UydY. 

By (|4.18c|) . we have that 

(4.55) |[/ y |<_^_( Jy + Xy ). 
Hence, 

\U(P)\ < \U(P 1 )\ + JL= |J(P X ) + x(P 1 ) - J(P) - x(P)\ 

(4.56) < |||e||| g( n) + ^=(5o + 4|||e||| g(n ) + 2L) by and flUZJ. 

To prove that and remain bounded, we use the bi-linearity embedded in the gov- 
erning equation (I2.14p . We use first the linearity of F(Z) with respect to the first variable 
and for almost every X G [A^,X r ], we get, after applying Gronwall's lemma, that 

f-y-Hx) 

\Z X (X,Y)\ < \Z x (X,y(X))\eMj \F(Z)(-,Z Y )\ dY) 

f-y-Hx) 

(4.57) = \V(X)\ exp( J \F(Z)( -,Z Y )\dY). 

Here F(Z)(-,W) denotes the matrix V i— > F(Z)(V,W) and we use any matrix norm as 
they are all equivalent. We also assume (as earlier) that Y < y~ 1 (X) (otherwise we have 
to interchange the bounds in the integral) and we denote Pi = (X, y^ 1 (X)). We have 
|V(X)| < |||6|||g;( n) . After using (|4.18all . (I4.18bll and (l4~55ll . we obtain that 

\F{Z){;Zy)\ < C(\t Y \ + \x Y \ + \U Y \ + \Jy\ + \K Y \) 
and we have used here the linearity of F(Z) with respect to its second variable. Hence, 

\F(Z)(; Z Y )\ = C (hxy + J Y ) + x y + J y + \U y \ 
< C(x Y + J Y ) 

for a constant C that depends only on c(U) and therefore only on |||0|||g(n) and L, by 
(gJaED . Hence, 

r y-Hx) j-y-Hx) 

J \F(Z)(-,Z Y )\dY <C I (x Y + J Y )dY 

= C{x{P l )-x{P) + J{P 1 )-J{P)) 
<C7(£ + 4|||e||| g(n) + 2L), 

by (14471) . (jUSJ) and (14541) . Combined with (I457D . it yields 

\Z X (X,Y)\<C 
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Figure 3. Construction of the global solution. 



for some constant C which depends only on L and |||@|||g(n)- Similarly, one proves the bound 
on Zy ■ The estimate (I4.45P follows from the estimate (|4.42l) in the proof of Theorem 14.81 
as the constant C in (|4.42l) depends only on L and |||©|||g(n), by (|4.44h . and 

\y - y(x)\ = \y-x + x(x~\x)) - y(x~\x))\ <l + |||e||| g(n) . 

□ 

4.3. Global existence. We obtain the following existence and uniqueness lemma. 

Lemma 4.10 (Existence and uniqueness on arbitrarily large rectangles). Given a rectan- 
gular domain £1 = [Xi,X r ] X [Yi,Y r ] and = (X, y, Z, V, W) in Gift), there exists a unique 
Z G Ti.(Q) such that 

@ = z.(x,y). 

Proof. Let N denote an integer that we will set later and 5 = Sr ^ Sl ■ For i = 0,...,N, 
let Si = id + si and we consider the sequence of points Pj = (Xi,Yi) = (X(si),Y(si)). 
For i,j = 0,...,N, we construct a grid which consists of the points P{j = (Xi t j,Yij) 
where X^j = Xi and Y^j = Yj, see Figure [31 We denote by by fijj the rectangle with 
diagonal points Pij and Pi + ij + ±. Let Q n denote the rectangle with diagonal points given 
by (Xo,Yq) and (X n ,Y n ). We prove by induction that there exists a unique Z G Tt(£l n ) 
such © = Z • (X, Y) (Here we use the same notation for G G G{Q) and (X,y) G C(f2) and 
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their restriction to J7 n which belong to Q(fl n ) and C(0 n ), respectively). On fij, we can 
choose N large enough and depending only on |||0|||g(ni) < sucn that 

si - s < 5 < cdiiem^))- 1 < c7(|||e||| g(Ql) )- 1 , 

and, by Theorem 14.81 there exists a unique solution Z £ TL{Q,i) such that Q = Z • (X, Y). 
We assume that there exists a unique solution Z £ Q n and prove that there exists a solution 
on O n+ i. On fl n<n , we get the existence of a unique solution by Theorem 14.81 as 

s n+ i -s n <5< CdllOIH^))- 1 < CdllOIH^))- 1 . 

For j = n—1, ... ,0, we construct iteratively the unique solution in £l n j and Qj, n as follows. 
We treat only the case of Fl n ,j- We assume the solution is known on O nj+1 , then we define 
© = (X,y, Z,V,W) G G(tt n ,j) as follows: The curve X(s),y(s) is given by 

X (s) = X n , y(s) = 2s-X n 

for \{Yj +X n )<s< |(Y i+ i + X n ), 

X{s) = 2s-Y j+1 , y(s)=Y j+1 

for \{Yj + \ + X n ) < s < \{Yj + i + X n+ i) and set 

Z(s) = Z{X(s),y(s)) for s G [^+X„), + 

V(X) = Z x (X,Y j+1 ) for a.e. X € [X n ,X n+1 ], 
W(Y) = Z Y (X n ,Y) for a.e. X E [i^+i]. 

Using Lemma [4791 we can check that |||0|||g(n n j) is bounded by a constant C2 that depends 
only on L and |||0|||g(r2)- We have 

i(Y- +1 + x n+1 - y,- - x n ) = + A'K+i) - y(sj) - X{s n )) < 25. 

Here we have used that X and y are Lipschitz with Lipschitz constant smaller than 2. By 
taking N large enough so that 25 is smaller that C(C2) _1 , we can apply Theorem 14.81 to 
U n j and obtain the existence of a unique solution in H(Q n ,j)~ Similarly we get the existence 
of a unique solution in H(Qj in ), Since 

^n+l = U (U" =o r2j i?1 ) U (Uj =Q Q n j), 

we have proved the existence of a unique solution in fin+i- □ 

In Lemma 14.91 we establish L°°-bounds on the derivatives on a strip domain. It turns 
out that we can also establish L 2 -bounds on the derivatives as stated in the next lemma. 
In this context, by L 2 -bounds, we mean that we can bound the integrals of the differential 
forms {Z^) 2 dX and {Z Y ) 2 dY along a curve in C. It is useful to have in mind that, for any 
given time T, we can find a curve (X,y) £ C which corresponds to this given time T, that 
is, t(X(s),y(s)) = T for all s € M.. Thus the L 2 -bound we now establish is fundamental to 
obtain L 2 -bounds in the initial set of coordinates. 
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Lemma 4.11 (A Gronwall lemma for curves). Given = [Xi,X r ] x [Y;,Y r ] ; Z £ 7i(Q) 
and (X,y) GC(fi). Then, for any (X,y) GC(fi), 

(4-58) \\z.(x,y)\\g m <c\\z.(x,y)\\ gm 

where C = C(\\(X ,y)\\^^, \\\Z • (X,y)\\\g(n)) ^ s a given increasing function with respect to 
both its arguments. 

Proof. Note that, for any function in / S Wy'°°(J2) (and respectively g £ W X °°(Q)), the 
forms f(X,Y)dX (respectively g(X, Y) dY) are well denned while the forms f(X,Y)dY 
(respectively g(X,Y) dX) are not. Given Z G H(O), we can consider the forms XJ 2 dX, 
U 2 dY, {Z a x ) 2 dX and (Z£) 2 dY. For any curve f = (X,y) £ C(fi), (X,y,Z,V,W) = 66 
£7(0), we have (by definition of the integral of a form along a curve and the definition of 

z»(x,y)) 

[_(U 2 dX + U 2 dY) = 2 [ " Zi(s)ds (as X + y = 2s), 
Jr is, 

and 

\_{Z a x ) 2 dX = [ Xr V a (X) 2 dX, [_(Z^) 2 dY = f W a (Y) 2 dY. 
Jr Jx, Jr Jy, 



We can rewrite 



\Z.(X, Y)\\l (n) = j_ Q[/ 2 (dX + dY) + (Z x ) 2 dX + (Z£) 2 dVj . 



Thus, we want to prove that 



(4.59) J_ {^U 2 (dX + dY) + {Z a x ) 2 dX + (Z^) 2 dY^j 

<C J Q[/ 2 (dX + dY) + (Z x ) 2 dX + (Z£) 2 dY^j . 



We decompose the proof into three steps. 



Step 1. We first prove that (I4.59P holds for small domains. We claim that there exist 
constants 5 and C, which depend uniquely on ||(<^>3^)|| C (q) and \\\Z • ('f ,J0lllc7(n) such 
that, for any rectangular domain Q = [Xi,X r ] x [l/,Y r ] with s r — si < 5, (|4.59l) holds. 
We denote by C a generic increasing function of ||(<^\50|| C (q) an d \\\Z • (X ,y)\\\g(ny By 
Lemma [4.91 we have 

H^1lL°°(n) + II-^xIIl 00 ^) + II-^yIIl 00 ^) ^ 

Let 

A = sup J Q[/ 2 (dX + dY) + (Z x ) 2 dX + (Z^) 2 dY^j 



THE NONLINEAR VARIATIONAL WAVE EQUATION 43 

where the supremum is taken over all T = (X,y) £ C(Q). Since Z is a solution of (|2.13|) . 
we have for a.e. X £ [Xi,X r ] and all Y £ \Yi,Y r ], that 

(x x - \)\X,Y) = (x x - \)\X,y{X)) + f (x x - \)x XY dY 
1 1 Jy(x) 1 



{x x - l -f(x,y(x)) 



+ r C ^((^-h 2 (U x + U Y ) + U x (x x - l) + U Y (xx- \))dY 
Jy(X) c(U) 2 2 2 



and 



(4.60) h Xx -hfdX< [ (x x -\) 2 dX + C ^ \{Z\f + (Z Y f)dXdY. 



For any Y £ [Yj,y r ], the integral f x r (Z X ) 2 (X,Y) dX can be seen as the integral of the 
form (Z x ) 2 dX on the piecewise linear path T going through the points (Xi,Y{), (Xi,Y), 
(X r ,Y), (X r ,Y r ), so that f x ;{Z x ) 2 {X,Y)dX < A. Similarly, f£{Z x ) 2 (X, Y) dY < A, 
for any X £ [X h X r ], Hence, (14611 yields 

j_{x x - i) 2 dX < J^x x - i) 2 dX + C(Y r - Y x + X r - XfiA 

< f {xx ~ \? dX + C5A 

as (Y r — Yi + X r — X[) = s r — si. By treating similarly the other components of Z x and 

Z Y , we get 

(4.61) f_(Z x ) 2 dX < f (Z x ) 2 dX + 6C5A and (_{Z Y ) 2 dY < I \z Y f dY + 6C5A. 
J r j r •/ r «/r 

For the component U, we have 

U 2 {X,Y) = U 2 {X,y{X)) + / 2UU Y dY 

Jy(x) 



iy(x) 

<-Y ,-Y 



<U 2 {X,y(X))+ [ U 2 dY+[ U^dY 
Jy(x) Jy(x) 



y(x) Jy(x) 

and it follows, as before, that 

(4.62) J U 2 dX < J U 2 dX + C5A. 

Similarly, we obtain 



(4.63) J u 2 dY< J 



U z dY< I U 2 dY + C5A. 
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After adding (14T6TL (14T621) . (14T6H and recalling that ds = \{dX + dY), we obtain 

/(^[/ 2 {dx + dY) + (z^) 2 dx + (zf-) 2 dy) 

Jr 2 

< [ (\u 2 (dX + dY) + (Z£) 2 + (Z^) 2 dY) + 13CM 
which yields, after taking the supremum over all curves F, 

(1 - 13C5)A < [ (-U 2 (dX + dY) + (Z x ) 2 dX + (Z^) 2 dY) 
and p~59D follows. 

Step 2. For an arbitrarily large rectangular domain £1 = [Xi,X r ] x [Yj,y r ], let us prove 
that p~59l) holds for the curves f = (X,y) G C(U) such that 

y(s) - X(s) > y(s) - X(s) for all s G (s h s r ), 

that is, the curve f is above T and intersects T only at the end points. Similarly one proves 
that (14391) holds for curves f = (X, y) G C(Q) such that y(s) - X(s) < y(s) - X(s) 
for all s G [s/,s r ]- For a constant K > that we will determine later, we have for a.e. 
X G [Xi,X r ], that 

e -K(y(x)-x ){xx _ ]_f { x,y {x)) _ e -Kmx)-x ){xx _ l )8(Jf)W) 

= / -Ke- K V~ x \ Xx - \) 2 dY+ e- K ^ x - Y \x x - \)x X y dY 

which implies, since Z is solution and by the estimates of Lemma l4~9l that 



(4.64) 



J e- KiY - X \xx ~ \? dX - j e~ K ^- x \x x - \) 2 dX 
r x r r y(x) 

< / / -Ke- K{ y- x \x x --) 2 dXdY 
Jx l Jy(x) 2 

r x r f y(x) 

+ C / e- Ki - Y - x \{Z a x ) 2 + {Z Y f) dXdY. 
Jx t Jy(x) 

Note that (|4.64l) corresponds to an application of Stokes's theorem to the domain bounded 
by the curves V and V. We treat in the same way each component of Z x and obtain that 

(4.65) j e~ K{Y - x \z x ) 2 dX - J e- Ki - Y - x \z a x ) 2 dX 

< / / -Ke' K{Y - x Hz a v ) 2 dXdY 



x> 

Xi Jy(x) 
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x r ,-y(x) 
x l Jy{x) 



+ C / e- K V- x X{Z a x ) 2 + (Z Y f) dXdY. 
Jx l Jy(x) 

As far as Z y is concerned, we get 

(4.66) J e^ Y - x \Zyf dY-j e~ K ^- x \z Y f dY 

r Y r r X(Y) 

<-/ / Ke- K{ y~ x \z\fdXdY 

JYi Jx(Y) 

f Y r r X(Y) 

+ C / e- K V- x X(Z a x f + (Z Y ) 2 ) dXdY. 

JYi Jx(Y) 



Let us prove that the sets 



J Xi < X < X r \Yi<Y <Y r 

Ml = \y(x)<Y<y(x) and M2 = \x(Y)<x<x(Y) 

are equal up to a set of zero measure. Let us consider (X,Y) G J\f\. We set s\ = X~ 1 (X), 
s 2 = y-\Y), s 3 = y- x (Y) and s 4 = X~ l (X). Since ^(X) < Y < y(X), we get 

y(x) = y(si) < y(s 2 ) = y(s 3 ) = y < y(s 4 ) = y(x). 

Hence, si < s 2 and S3 < S4, which implies 

X = X{ Sl ) < X{s 2 ) = X(Y), X(Y) = X(s 3 ) < X{s A ) = X 

and therefore X(Y) < X < X(Y). Thus we have prove that Mi C M2 up to a set of zero 
measure. Similarly, one proves the reverse inclusion. Hence, by adding (|4.65|) and (|4.66l) . 

we get 

(4.67) J_e- KiY ' x \(Z x ) 2 dX + (Z Y ) 2 dY) - J e~ K{y - x) ((Z a x ) 2 dX + (Z Y ) 2 dY) 

<-K f e- K(Y - x) ((Z x ) 2 + (Z Y ) 2 ) dXdY + C [ e- K ^ Y - x \(Z x ) 2 + (Z Y ) 2 )dXdY. 
J Mi J Mi 

As far as U is concerned, we proceed in the same way and get 

e -K(Y- X)U 2 dx _ f e -K(Y- X)u 2 dx 

"X r r-y(x) 
q Jy(x) 



r 



(4.68) 
and 



= / / e- K{Y ~ x) (-KU 2 + 2UU Y )dXdY 
Jxi Jy(x) 

< [ e- K{Y - x) (-KU 2 + U 2 + U Y )dXdY 
J Mi 



J e- K V~ x ^U 2 dY- J e- K( ^- x ^U 2 dY 
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(4.69) < f e~ K{Y ~ x) {-KU 2 + U 2 + U 2 x )dXdY. 

J Mi 

Combining (l4~65ll . (l4~66i (l4~68ll . (l4~68ll . we get 

(4.70) / e" x(y - x) (-?7 2 (dX + dY) + (Z£) 2 dX + (Z£) 2 dY) 
Jr 2 

- / e- Ki > Y - x \-U 2 {dX + dY) + {Z a x ) 2 dX + {Z$) 2 dY) 
Jr 2 

< / {C - K)e- K{Y - x) {U 2 + {Z a x f + {Z^f)dXdY. 
JMx 

We choose K sufficiently large so that the right-hand side in (|4.70l) is negative and we 
obtain that 

e -. K \\X-y\\ LOO I A V 2 ( dx + dY) + ( Z a )2 dx + (z a )2 rfy j 

Jr 2 

/ (ic/ 2 (dx + dy) + (z£) 2 dx + {z^f dY) 

Jr 2 



< e K\\x-y\\ 



and (I459D follows. 



Step 3. Given any rectangle = [Xj,X r ] x [y,y], we consider a sequence of rectangular 
domains Oj = [Xi, Xj+i] x [y , y+i] for % = 0, . . . , N— 1 where X,- L and y are increasing and 
X = X, y = Y h X N = X r , Y N = y and such that (X,y), (X,y) belong to g(Qi) for 
s G [sj, Si_|_i]. We construct the sequence of rectangles such that either Sj + i — Sj < 5 (and 
Step 1 applies) or y(s) - X(s) < y(s) -X(s) or y(s) -X(s)< y(s) - X{s) for a € [s fj s i+1 ] 
(and Step 2 applies). Hence, 

N-1 

1 2 II rr , ■r, 1 1 2 



^•(^)&n) = EII Z *(*>^ 



g(n) ii ^ w '^^iig(^) 

i=0 
AT-l 



< y, c \\ z ' {xmIm ( b y ste P s 1 and 2 ) 



i=0 

<C||Z.(^,^)|| 2 



10(H) • 

We can construct the sequence of rectangles as follows. Let N be an integer such that 
Sl N Sr < | and we set Sj = si + j| for j = 0,...,N, We take so = si and define Sj 
iteratively: Given Si and ji G {0, . . . , N — 1} such that > i, A'(sj) = <V(sj), 3^-Si) = 3^(s«) 
and Si G [s^ , s^+i] • If ji + 1 = X, we set X = i + 1, Sj+i = s r and we are done. Otherwise, 
there exists an index k > ji + 1 such that 3^(s) — (s) < y(s) — X(s) for all s G [sj+i, s^) 
or lP(s) — ^(s) > y(s) — X(s) for all s G [sj+i, s^) and there exists s G [s^, s^+i] such that 
^(s)-^(s) = y(s)-X(s) (which implies that Af(a) = X(s)smd y(s) = y(s)). We then set 
ji+l = k and choose Si+i G [sfc,s"fc+l] such that ^f(si + i) = X(si+i) and ^(si+i) = ^(si+i)- 
Since ji > i, the iteration stops in a finite number of steps. □ 
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Si = SQ 



Figure 4. The interval [si, s r ] is divided into large intervals (in this example 
[so,si] and [s3,S4]) where one curve is over the other and small intervals of 
length smaller than 5 (in this example [si,S2] and [412,83]) where the curves 
can cross an arbitrarily number of times. 

Given two solutions Z and Z, we want to compare along curves in C the forms Zx dX 
and ZxY dY with Zx dX and ZxYdY, respectively. By using the same argument as in 
the proof above, we obtain the following stability result in L 2 . This is a stronger result 
than the one that could be established from the fixed point argument in Lemma 131)1 as the 
latter would only hold in L°° . 

Lemma 4.12 (Stability in L 2 ). Given = [X h X r \ x \Y h Y r ], Z,Z G H(P) and (X,y) G 
C(O). Then, for any (X,y) G C(O), 

(4.71) \\(Z - Z) . (X,y)\\ g{ n) < C\\(Z - z) . (x,y)\\ g(a) 

where C = C(\\(X, ^)llc(n), 111^ • (X ,y)\\\g(n)), \\\Z • (X, y)\\\g(n)) is a given increasing 
function with respect to both its arguments. 

In the definition below of global solutions we include a condition about the decay of 
the solutions along the diagonal. This condition is necessary to guarantee that, given a 
solution, the curves which correspond to a given time T belong to C. 

Definition 4.13 (Global solutions). Let Tt be the set of all functions Z G W^^°(M?) such 
that 

(i) Z G Tl(£l) for all rectangular domains and 

(ii) there exists a curve (^,3^) G C such that Z • (^,3^) G Q . 

The condition (ii), which corresponds to a decay condition, does not depend on the 
particular curve for which it holds, as the next lemma shows. In particular, we can replace 
condition (iv) in Definition 14.131 by the requirement that Z • (Xd,yd) G Q for the diagonal 
(Y = X), which is given by Xd(s) = 3^rf(s) = s. We then denote 

\\Z\\ n = \\Z.(X d ,y d )\\ g and |||Z||| w = |||Z.(* d ,3^|||g. 
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Lemma 4.14. Given Z G TL, we have Z • (X,y) G Q for any curve (X,y) G C. Moreover, 
the limit lims^oo J (X (s) , y (s)) is independent of the curve (X,y) G C. 

In this lemma we denote as before Z 4 by J where G = (X ,y, Z,V,W) = Z • (X,y). 
Later, we will see that the limit of J at infinity corresponds to the total energy and the 
lemma would allow us to prove that the total energy is conserved. 

Proof of Lemma \4-14\ For an y curve (X,y) G C, we have to prove that 

(4.72) \\\Z»(X,y)\\\g <oo and \\Z • (X,y)\\ g < oo 
and 

(4.73) lim J(s) = 



where J = Z A with (X,y,Z,V,W) = Z • (X,y). For any real number s G R that will 
eventually tend to infinity and denote Us = [X(— s), X(s)] x [y(—s),y(s)]. Let 



(4.74) 

and 

(4.75) 



'y-\y(s)) a y(x(s))<y(s), 

X~ 1 (X(s)) otherwise, 



y-Hy(-s)) if y{x(-s)) < y(-s), 

X~ l (X{—s)) otherwise, 



see Figure[5]for an example. One can check that by construction s m i n < — s < s < s max and 
we denote Q s = [^(smin),^(s m ax)] x [y(smm),y(sm& x )]- We have C ttg- We construct 
the curve which consists of a (vertical or horizontal) straight line joining (X(s m i n ),y(s m i n )) 
and (X(s),y(s)), the curve (X(s),y(s)) for s G [— s, s] and another (vertical or horizontal) 
straight line joining (X(s),y(s)) and (X(s max ), y(s max )), see Figure We denote by 
(X,y) this curve and we have that (X,y) and (X,y) belong to Q(Cls). By Lemma [479l we 
get 

\\\z*(x,y)\\\ gm <\\\z.(xM\ gi n s ) 

< CiCIK^P,^)!^^), IIIQllbco^)) < Ci(ll(^,3>)llc |||€»|||^) 
and by letting s tend to infinity, we get \\\Z • {X ,y)%g < oo. By Lemma f4. Ill we get 

(4.76) \\z. (x,y)\\ gm <\\z. (x, y)\\ g{Cls) <c\\z.(x, y)\\ gm <c\\z.(x, y)\\ g 

where the constant C depends on \\(X, Y)^ g ,^ ^ and |||Z»(A', 3 ; )|||g(Q_), that is, on ||(X, Y")||g 
and \\\Z • {X,y)\\\g, which are independent on s. By letting s tend to infinity in (|4.76l) . we 
get || Z • (X, y)\\g < oo. It remains to prove (I4.73p . We know that J is positive. We denote 
J(s) = Z±{s) with (X ,y, Z,V,W) = Z • (X, y) and, slightly abusing the notation, we de- 
note also by J, J(X,Y) = Z 4 (X,Y). For any s£l, let ax = X~ 1 X(s) and s 2 = y _1 y(s). 
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S: 



Figure 5. Prolongation of the curve T. 



If si < s 2 , then X(s) = X(s\) < X{s 2 ) and y{s\) < y{s) = y(s 2 ). By the monotonicity of 
J(X, Y"), we get 

j(x( Sl ),y(si)) < j(x( s ),y(s)) < J(x( S2 ),y(s 2 )). 

If S2 < Si, we get a similar result so that, finally, 

(4.77a) mm{J(X(s 1 ),y( Sl )),J(X(s 2 ),y(s 2 ))} < J(X{s),y(s)) 

and 

(4.77b) J(X(s),y(s)) <m a x{J(X( Sl ),y( Sl )),J(X(s 2 ),y(s 2 ))}. 

Since 

|si -s\< \x( Sl ) - Sl | + |^( s ) - s \ < \\(x,y)\\ c + \\(x,y)\\ c , 

we have that lim^-too s± = ±oo similarly we obtain that lim^j-oo s 2 = ±oo. Hence, (I4.77P 
yields 

lim J(X(s),y(s)) = lim J(Af(a),y(s)). 

s— >±oo s— >±oo 

Thus, these limits are independent of the curve T which is chosen. The existence of the 
limits is guaranteed by the monotonicity and boudedness of J. The identity (|4.73|) follows 
from (pTQj) . □ 

From Lemma 14.101 we infer the following global existence theorem for the equivalent 
system. 

Theorem 4.15 (Existence and uniqueness of global solutions). For any initial data 

(X, y, Z, V, W) in Q, there exists a unique solution Z £ Tt such that = Z • (X,y). We 

denote this solution mapping by S: Q — > 7i. 
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5. Semigroup of solution St in T 

From a solution function Z G TL in the whole plane, we want to extract the data at 
a given time. It is enough to do it at t = 0, and the definition below describes how we 
proceed. 

Definition 5.1. Given Z G TL, we define 

(5.1) X(s) = sup{X G R | t(X', 2s - X') < for all X' < X} 

and y(s) = 2s- X(s). Then, we have (X(s),y(s)) G C and Z • (X,y) G Q . We denote 
by E the mapping from TL to Qq that associates to any Z £TL the element Z • (X, y) G Go 
as defined here. 

Proof of the well-posedness of Definition I5.il First we prove that X is increasing. Let X n 
be a sequence such that X n < X(s) and X n — > X(s). We have t(X n , 2s — X n ) < for any 
s > s. Since t is decreasing with respect to the second variable (as ty < 0), we have that 
t(X n ,2s — X n ) < and therefore X n < X(s). After letting n tend to infinity, we obtain 
X(s) < X(s) so that X is an increasing function. Let us prove that X is Lipschitz with 
Lipschitz coefficient smaller than 2. Let us assume the opposite, i.e., there exists s > s such 
that 

(5.2) X(s)-X(s)>2{s-s). 

It implies that y(s) > y(s). Since tx > and ty < 0, we have, for any (X, Y) G 

[x(s),x(s)]x[y(s),y(s)] 

o = t(x( s ),y(s)) < t(x,y( s )) < t(x,Y) < t(x,y(s)) < t(x(s),y(s)) = o 

and, therefore t(X,Y) = on $7 = [Af(s),^(s)] x [3^00, y(s)]. Let us consider the point 
(X,Y) G Q for Y = y(s) and X = 2s - y(s). We have t(X,Y) = 0, X + Y = 2s and 
X < X(s) : which contradict the definition of X at s. Thus, we have proved that X is 
Lipschitz. To show that (^,3^) G C, it remains to prove that — 3^lliL°°(]R) < 00 ■ We 
claim that there there exists L such that 

(5.3) liminft(X + L,X) > 1 

X— »oo 

for any L > L. Let us prove this claim. By using the fact that xx = c(U)tx and c(U) > ^ 
we get 

rX+L 

t(X + L,X) = t(X,X) + / t x (X,X)dX 

Jx 

(5-4) > -|||Z||| W + ±- + - [ + (x x -h(X, X) dX. 

Ik k Jx 2 

We look at the domain £lx,L = [X, X + L\ x [X, X + L\. We consider the curve Xa(s) = 
yd(s) = s (the diagonal) and the curve (X,y) which consists of a horizontal and a vertical 
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segment given by 

{x(s) = 2s-X, y{s) = X for s £ [0, X + -§] 

[X(s) = X + L, y(s) = 2s - (X + L) forse [X + %,X + L]. 

By Lemma f4. Ill we get 

rX+L i 

(5-5) (x x - -) 2 (x,x)dx < \\z.(X,y)\\ 2 g(nxL) <c\\z.(x d ,y d )f g{nxL) 

where the constant C depends on L and \\\Z • (X d , yd)\\\g(n x l) which 1S bounded by |||Z|||%. 
We have 

rX+L 

lim \\Z.(X d ,y d )f g{n ) = hm / (C/ 2 (X,X)+(Z^) 2 (X,X)+(Z?-) 2 (X,X))dX = 0, 

and therefore (15. 5|) and (15. 4p yield 

liminfi(X + L,X) > -IIZIL. + — 
X—oo " n 2k 

which, for L large enough, implies (|5.3p . Using the same type of argument, we prove that 



there exists L such that 

(5.6) liminft(X + L,X) > 1, limsupt(X - L, X) < -1, 

X-^oo ' X^oo 

(5.7) liminf t{X + L,X) > 1, liminf t(X — L,X) < — 1. 

V— ►— oo V— >— oo 

Let us prove that limsup s ^ 00 (Af(s) — s) < -|. We assume the opposite and then, there 

exists s <E K such that i(s + f,s - f) > 1, by (151)1) . and > s + f . It impli 



les 



that y(s) = 2s — X{s) < s — 2", and, using the monotonicity of t (that is tx > and 
iy < 0), we get 1 < t(s + f,s - §) < t(^(s),y(s)) = 0, which is a contradiction. 
Similarly one proves that limmf s ^ OQ (X(s) — s) > — limsup s ^_ OG (A'(s) — s) < -| and 
liminf s ^_ 00 (^'(s) — s) > — k and it follows that 

limsup \X(s) — s\ < — . 

s— >±oo 2 

Hence, the condition (I2.15aj) is satisfied and (X,y) £ C. By Lemma 14.141 we have 
(X,y,Z,V,W) = Z • (X,y) £ Q and by construction Z x {s) = t(X(s),y(s)) = so 
that Z • (X, y) £ Qq. □ 

Definition 5.2. Given any T, let us introduce the mapping ty: TL — > 7i defined as follows. 
For any Z £ TL, let t^(Z) = Z £TL be given by 

i(X,Y) = t(X,Y)-T 

and 

x(X,Y) =x(X,Y), U(X,Y) = U(X,Y), 

J(X,Y) = J(X,Y), K(X,Y) = K(X,Y). 
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We have 

(5.8) tx+T' = ° tx 1 ' • 

We have now denned the following mappings: 

c s 

D 



C S 

(5.9) T^^Go^^Hj 

I A I." 



The mapping is used to extract the solution at any given time T by only using the 
operator E, which is designed for time zero. Indeed, by taking Eoty, we recover an element 
in Go which corresponds to the solution at time T. In the following lemma, we prove that 
T and TL are in bijection, which also justify the introduction of T: It is a consistent way to 
parametrize initial data: To any element in there corresponds a unique solution in TL, 
and vice versa. The Go does not fit that role as Go and TL are not in bijection. 

Lemma 5.3. We have 

(5.10) CoDoE = E, DoC = Id 

and 

(5.11) EoSoC = C, SoE = Id. 

It follows that S o C = (D o E)~ x and the sets T and TL are in bijection. 

Proof. Step 1. We prove KWh . Given Z G TL, let (X,y,Z,V,W) = E(Z),_(^_i,^ 2 )_ = 
T)(X,y,Z,V,W) and {X,y,Z,V,W) = Ctyujk)- We want to prove that (X,y,Z,V,W) = 
(X,y,Z,V,W). We have to prove that X = X, the rest will easily follow. For any s£l, 
we claim that for any couple (X, Y) such that X < X(s) and X + Y = 2s then we have 
either 

(5.12) xx(X) < xi{X{s)) or x 2 (Y) > x 2 (y(s)). 

Let us assume the opposite, that is, there exist s, X and Y such that X < X(s), X + Y = 2s 
and 

x x (X) = Xl (X(s)) = x(s) = x 2 (y(s)) = x 2 (Y). 

Here,_a;(s) denotes Z 2 (s), see (1^271) . Let s = X~ l {X) and s x = y^iY). Since X < X{s) 
and Y > J^(s), we have sq < s < s\. We have 

x(s ) = x 1 (X(s )) = xi(X) = x(s) 

and, similarly, we obtain that x(si) = x(s). We consider the rectangular domain Q = 
[X(so), X(s\)] x [y(so),y(si)]. Since x(s) = x(so) = x(s\) for all s £ [sq,«i], we have 
x = on [so,si] because x is nondecreasing. We have x = = V 2 (X)X + W 2 (y)y on 
[s ,si], which implies that V 2 (X) = for a.e. X £ [X (s ) , X (si)] and W 2 (Y) = for a.e. 
Y G [y(s ),y( Sl )]. By (l^8bl and (l3^aTl it implies Vi = V 2 = V 3 = on [Af(s ),#(si)] and 
Wi = W 2 = W 3 = on [3 ; (so),^(si)]. Then, we can check that Z given by 

t(X,Y) = 0, x(X,Y)=x(s), U(X,Y) = U{s) 
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and 

J(X,Y) = Ji(X) + J 2 (Y), K(X,Y) = K X {X) + K 2 {Y) 

is a solution to (l2~T3l) in (that is, Z G and Z • (#,;)>) = (X,y,Z, V, VV). By 

uniqueness of the solution, we get Z = Z. In particular, we have t(X,Y) = such that 
X + Y = 2s and X < X(s), which contradicts the definition of X(s) given by (15,11) . This 
concludes the proof of the claim (15.121) . Since xi(X(s)) = x 2 {2s — X(s)) we get, by (|3.27l) . 
that 

X(s) < X(s). 
We have by the continuity of x\ and x 2 that 
(5.13) Xl (X(s))=x 2 (y(s)). 

Let us assume that X(s) < X(s), then, by the claim (|5.12l) we have proved, we have 
either si (#(*)) < xi{X(s)) or x 2 {y{s)) > x 2 (y(s)). If a?i (#(«)) < a* (#(*)), then, as 

y(s) > y( s ) 

Xl (x( s )) < Xl (x( s )) = x 2 (y(s)) < x 2 (y( s )) 

which contradicts (I5.13p . Similarly, we check that if x 2 (y(s)) > x 2 (y(s)) then we obtain 
a contradiction to (I5.13p . Hence, X = X and therefore 3^ = 3^- Then, x(s) = x 2 {X(s)) = 
x 2 (X(s)) = x(s) and similarly, we treat the other components of Z. From the definitions 
of C and D, we have that V = V and VV = W. Hence, we have proved that C o D o E = E. 
The fact that D o C = Id follows directly from the definitions of C and D. 

Step 2. We prove ijCTjl . Given (^1,^2) G let (X,y,Z,V,W) =_C(^i,^), Z = 
S(X,y,Z,V,W) and (X,y,Z,V,W) = E(Z). We have to prove that ^ = AT, the rest 
will easily follow. Since Z G TL is a solution with data (X,y,Z,V,W) G £/o, we have 
i(^f(s), y(s)) = 0. Hence, from the definition of E, we get X(s) < X(s), Assume that there 
exists sGK such that X(s) < X(s). By the definition of E, we have t(X(s),y(s)) = 0. Let 

s = X- 1 (X(s)) and *i = 3T 1 (?(*))■ 

We have 

X(s ) = X(s) < X(s) and y(s 1 )=y{ s )>y( s ), 

and therefore sq < s < s\. Due to the monotonicity of t(X, Y) (that is, tx > and < 0), 
since t(AT(a ),y(so)) = *(^(a),3?(a)) = t(AT(si),y(si)) = 0, we get that t(X,Y) = on 
the rectangle = [X(sq),X(si)] x [3^(so), It implies that xjf = c(U)tx = and 

= c(U)ty = on f!. Hence, the function x(X, Y) is constant on f2 and we have 

Xl (x( s )) = xx(x( So )) = x(x( So ),y(s )) = x (x( Sl ),y( Sl )) = x 2 (y( Sl )) = x 2 (y(s)). 

However, the fact that xi(X(s)) = x 2 (y(s)) and X(s) < X(s) contradicts the definition of 
X in (|3.22l) . and therefore we have proved that X = X. Then, y = y and 

z{ s ) = z(x( s ),y(s)) = z(x( s ),y( s )) = z( s ). 

Similarly, one proves that V = V and VV = W. Thus we have proved that E o S o C = C. 
The fact that S o E = Id follows from the uniqueness of the solution for a given data. □ 
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Definition 5.4. For any T > 0, we define the mapping St' T ' — ► T by 

(5.14) S T = DoEot T oSoC. 

Theorem 5.5. The mapping St is a semigroup. 

Proof. We have 

St ° St> = DoEotfoSoCoDoEo tjv o S o C 

= D o E o t t ° tjv o S o C (by Lemma 15. 3|) 

= DoEot T+T / oSoC (by (JSIHJ)) 

= St+t' ■ 

□ 

6. Returning to the original variables 
Definition 6.1. Given tp = (ipi,tp 2 ) £ 3~ ', we define (u, R, S, (jl,u) £ V as 
(6.1a) u(x) = Ux(X) ifxi{X) = x 

or, equivalently, 

(6.1b) u(x) = U 2 (X) ifx 2 (X) = x 

(6.1c) n=(x 1 )#{J{(X)dX), 

(6.1d) v=(x 2 ) # (4(Y)dY), 

(6.1e) R(x) dx = (xi) # {2c(U 1 {X))V 1 (X) dX) , 

(6.1f) S(x) dx = (x 2 )# {-2c{U 2 (Y))V 2 (Y) dY) . 

The relations (|6.1e|) and (|6.1f|) are equivalent to 

(6.2a) R{x 1 {X))x' l {X) = 2c(U 1 (X))V 1 {X) 

and 

(6.2b) S{x 2 (Y))x' 2 (Y) = 2c{U 2 (Y))V 2 (Y) 

for a.e. X and Y . We denote by M: T — » V the mapping that to any tp £ T associates 
(u, R, S, fj,,u) G T as defined above. 

We have to prove that the measures (aci) # (c(£7i(X))T<i(X) dX) and (x 2 ) # (-c(U 2 (Y))V 2 (Y) dY) 
are absolutely continuous with respect to the Lebesgue measure and that (u, R, S, /u, v) be- 
longs to V so that the definition is well-posed. It will be done in the proof of the following 
lemma where an equivalent definition of the mapping M is given. 



lr The push-forward of a measure A by a function / is the measure /#A defined by /#A(B) = A(/ 1 (B)) 
for Borel sets B. 
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Lemma 6.2. Given if> = (ipi,ip2) £ F > let (u,R,S,fi,u) = M^i,^) as defined in Def- 
inition Then, for any = (X ,y, Z,V,W) G Go such that (ipi,ipi) = D0, we have 



(6.3a) u(x) = U(s) if x = x(s) 

and 

(6.3b) n = x#(y A {X(s))X{s)ds), 

(6.3c) v = x # {W A {y{s))y{s)ds), 

(6.3d) R{x) dx = x # (2c(U{s))V 3 (X(s))X(s) ds^j , 

(6.3e) S(x) dx = x # (-2c(U \s))W 3 {y (s))y (s) ds} . 

The relations (|6.3djl and (|6.3eP are equivalent to 

(6.4a) R(x(s))V 2 (X(s)) = c(U(s))V 3 (X(s)) for any s such that X(s) > 
and 

(6.4b) S(x(s))W 2 (X(s)) = -c{U{s))W 3 (X(s)) for any s such that y(s) > 0, 

respectively. 

Proof. We decompose the proof into 5 steps. 

Step 1. We prove that ([6J]) imply ((631). If x = cci(X), let s = X~ l {X). Then, we have 
x = x 1 {X{s)) = x(s) and U X (X) = Ui(X{s)) = U{s). Hence, (RUa!) implies (lOaD . For any 
measurable set A, we have 

fi(A) = [ J[{X)dX 

= [ JiW.))*(.)d. (after a change of variables) 

J(x 1 oA')- 1 ( J 4) 

= / V 4 {X(s))X(s)ds (by the definition of D) 

Jx- 1 {A) 

and therefore /i = x#(Vi(X(s))X(s) ds). One proves in the same way the other identities 
in ([OH . 

Step 2. We prove that u is a well-defined function L 2 that is Holder continuous with 
exponent 1/2. Given x such that x = a;(so) = ^(si). Since x is nondecreasing, it implies 
that x = V 2 (^)<Y + W 2 {y)y = in [a ,si]. Hence, V 2 (X)X = W 2 (^)3> = as both 
quantities are positive. By (13.8ap . it implies that V 3 (X)X = W 3 (y)y = and therefore 

u = v 3 (x)x + w 3 (y)y = o 

in [so,si] an( l U(sq) = U(si). The definition of u is therefore well-posed. We have 



u 2 (x)dx= / u 2 (x(s))x(s) ds < \\U\ 
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and u e L 2 (M). We have 
(6.5) u{x(s)) - u{x(s)) = I U(s)ds = I (Va(X)X + W 3 (y)y)d8. 



Since £ 5o, we have t(s) = which implies that V±(X)X = W\{y)y and, therefore, 
V 2 {X)X = W 2 (y)y, by dObb . and 

(6.6) x = 2v 2 (x)x = 2w 2 (y)y. 

By using the Cauchy-Schwarz inequality and (I3.8ap . we get 

rs / rs \ 1/2 / ns \ 1/2 

y iv 3 (^)i^ds< (y v|(^)^dsj (j xdsj 

( r s \ 1/2 

(6.7) - C \J ±ds ) = c ( x ( s )- x (s)) 1/2 

where the constant C depends on |||0|||g and \X(s) — X(s)\. Similarly, one proves that 
// \W3(y)\yds < C{x{s)-x{s)) 1 / 2 . Hence, JEH) implies that u is locally Holder continuous 
with exponent 1/2. 

Step 3. We show that the measures x# (c(U(s))V 3 (X(s))X(s) ds^j and 
x # (-c(U(s))W 3 (y(s))y(s) ds^j are absolutely continuous and (|6.4aD and (|6.4b|) hold. The 

inequality (|6.7|) proves that the measure V 3 (X)X ds is absolutely continuous with re- 
spect to xds. For any set A of zero measure, we have f x -i^ x ds = and therefore 

fx- 1 (A) V3(X)X ds = 0. It follows that 

x # (2c(U(s))V 3 (X(s))X(s)ds) (A) = [ 2c(U)V 3 (X)X ds = 
V ' Jx- 1 {A) 



and the measure x# f 2c{U(s))V%{X(s))X(s) dsj is absolutely continuous. In the same way, 

one proves that x# (^—2c(U(s))W 3 (y(s))y(s)ds S j is absolutely continuous and R(x) and 
S(x) as given by (lOdl and (IQel) are well-defined. We have 



(6.8) / R(x(s))x(s) ds = / 2c(U(s))V 3 {X(s))X{s)ds 

for any measurable set A. For any measurable set B, we have the decomposition x~ l (x(B)) = 
BU(B c C\x~ 1 (x(B))) . Let prove that the set B c Hx~ 1 (x(B)) has measure zero with respect to 
x(s) ds. We consider a point s 6 B c C\ x~ 1 (x(B)), there exists s £ B such that x(s) = x(s). 
Since s is increasing, it implies that x(s) = and therefore fs^x^ixiB)) x ( s ^ ds = so that 
the set B c r\x~ 1 (x(B)) has zero measure with respect to x(s)ds. Since V 3 (X(s))X(s) ds is 
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absolutely continuous with respect to x(s) ds, it implies that JB c r\x~ 1 {x{B)) ^(.^( s ))^-( s ) ds = 
0. Hence, taking A = x(B) in (|6.8p . we obtain 



R(x(s))x(s)ds = / 2c(U(s))V 3 (X{s))X(s)ds 

B Jb 

for any Borel set B. Hence, 

(6.9) R(x(s))x(s) = 2c(U(s))V 3 (X(s))X(s) 
which yields 

(6.10) R(x(s))V 2 (X(s)) = c(U(s))V 3 (X(s)) 

after simplifying by X{s). Similarly, we obtain (I6.4b|) . 

Step 4. We show that R and 5 belong to L 2 and u x = We hi 

R 2 (x)dx = / R 2 (x(s))xds 



2 / R 2 (x(s))V 2 (X(s))Xds 

2 / mmimxjs 



'{seR|v 2 (A-(s))>o} V 2 (X(s)) 
< 4 / V 4 (X)Xds (by (IQaTl and dOajl ) 

<4||J|| LOO(R) <4|||9||b 

because J(s) = V4,(X)X + VV , 4(3 ; )3 ; and V4 and W4 are positive. Hence R £ L 2 and, 
similarly, one proves that SG L 2 , For any smooth function 4> with compact support, we 
have 



u(f) x dx= / u(x(s))4> x (x(s))x(s) ds = / U(s)4>(x(s)) s ds. 
1 Jr Jr 

After integrating by parts, it yields 

[ U( f> x dx = - [ U(s)<P(x(s))ds= [ (V 3 (X)X + W 3 (y)y)(p(x(s))ds 
Jr Jr Jr 

' (R(x(s))V 2 (X)X - S(x(s))W 2 (y)y)cP(x(s)) ds 



c(U) 

^—(R(x(s)) - S(x(s))<f>(x(s))x ds (by (|M 

— — — <pdx, 
2c(u) 

after a change of variables. Hence, u x = in the sense of distribution. 

Step 5. We show that /x ac = jR 2 dx and f ac = \S 2 dx. Let 
(6.11) A = {s I V 2 (^(s)) > 0} and 5 = {x(A c )) c . 
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We have 

meas(-B c ) = / x ds = 

because x = almost everywhere on A c , by (|6.6p . The set B has therefore full measure. 
We have x~ 1 (B) C A. Indeed, for any s G we have / x(s) for all s € A c . For 

any measurable subset B' C B , we have, by definition of that 



(6.12) KB')= Vi(X(s))X{s)ds. 

Hence, since x~ 1 {B l ) C A, 

i t>i\ [ Va{X)V 2 {X) ■ 



I est) 

A-i(B') ^V 2 (<<) 

- / i? 2 (x(s))x(s)ds fbv dEaj) and (lOll) 



R 2 dx 



B' 



and therefore ^ ac = ^R 2 dx. Similarly, one proves that ^ ac = |5 2 dx. □ 

Lemma 6.3. Given (uo, Rq, Sq, fio, uq) G 2?, /ei ms denote (u, i?, S, /i, v)(T) = M o S T o 
L(ito, -Roj Sb> [J-o, ^o) and Z = S o L(uo, i?0) S'o) ^o)- Then, we have 

(6.13) u(t(x, y), x(x, y)) = u(x, y) 

for all (X,Y) G R 2 , and 

(6.14a) R(t(X,Y),x(X,Y))x x (X,Y) = c(U(X,Y))U x (X,Y), 

(6.14b) S(t(X, Y),x(X, Y))x Y (X, Y) = -c(U(X, Y))U Y (X, Y), 

for almost all (X, Y) G M? such that xx(X,Y) > and xy(X,Y) > 0. We have 

(6.15) u t = l(R + S) and u x = -^—{R-S) 

2 2c{u) 

in the sense of distributions. 

Proof. We consider a solution Z G TC. Given (X, Y) G M 2 , let us denote i = t(X,Y) and 
x = x(X,Y). Let (X,y,Z,V,W) = Eot,-(2). By definition, we have t(X(s),y(s)) = i, 
and, slightly abusing the notation, x(s) = Zi(s) = x(X(s),y(s)) and U(s) = Z^{s) = 
U(X(s),y(s)) for all s G K. By Lemma [6T2l we have u(i,x) = U(s) for any s such that 
x(s) = x. It implies that, for any s such that 

(6.16a) t(X{s),y(s)) = t = t(X,Y) 

and 

(6.16b) x(X(s),y(s)) = x = x(X, Y), 
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we have 

u(t,x) = U(X(s),y(s)). 
Then, (I6.13P will be proved once we have proved that 

(6.17) U(X(s),y(s)) = U(X,Y). 

Let us prove that when (16.16P hold, then either (X,Y) = (X(s),y(s)) or 

(6.18) x x = x Y = U x = Uy = 0, 

in the rectangle with corners at (X,Y) and (X(s),y(s)), so that (I6.17P holds in both cases. 
We consider first the case where X(s) < X and y(s) < Y . Since the function x is increasing 
in the X and Y directions, we must have xx = xy = in the rectangle X] x [3^00, Y] 

and, by (|4.18cl) . Ux = Uy = in the same rectangle so that U is constant and we have 
proved (|6.17|) . In the case where X(s) < X and y(s) > Y, since the function t is increasing 
in the X direction and decreasing in the Y direction, it follows that tx = ty = in the 
rectangle [^(s),!] x [Y, ^y(s)]. Hence, xx = xy = and, as before, we prove (|6.17|) . The 
other cases can be treated in the same way and this concludes the proof of (I6.17P and 
therefore (I6.13P holds. Let us prove (|6.14ap . By (|6.4a|) and the definition of E, we get 

R(i,x(s))x x (X(s),y(s)) = c(U(X(s),y(s)))U x (X( S ),y(s)) 

so that 

(6.19) R(t(X,Y),x(X,Y))x x (X(s),y(s)) = c(U(X(s),y(s)))U x (X(s),y(s)) 

for any s such that (|6.16|) holds. We have proved that when (|6.16p is satisfied, then either 
(X,Y) = (X( s),y(s)) or (I6~T8D holds. Hence, (l6~T3aD follows from (I6J9D . Similarly, one 
proves (I6.14bp . For any smooth function (f>(t,x) with compact support, we have 



(6.20) / u{t,x)(t> t {t,x) = 

u{t(X, Y),x(X, Y))4> t {t{X, Y),x(X, Y)){t x xy - x x ty) dXdY 



where we have used (14. 18a|) and (I6.13p . By differentiating the function <f>(t(X, Y),x(X, Y)) 
with respect to X and Y, we get that 

c/>(t(X, Y), x(X, Y)) x xy - <f>{t{X, Y),x(X, Y)) Y x x = Mt(X, Y),x(X, Y)){t x x Y - t Y x x ). 
We insert this identity in (I6.20P and obtain, after integrating by parts, 



u(t, x)4> t {t, x) dtdx = - / ((Uxy)x - (Ux x )y)(X, Y)</>(t(X, Y),x(X, Y)) dXdY 

Jr 2 

= - [ {{U x x Y -U Y x x ){X,Y)(l){t{X,Y),x{X,Y))dXdY. 
Jr 2 

We use (|6.14D and get 

/ u(t, x) (fit (t,x) dtdx = — [ (( — (fi) o (t, x) x x xy dXdY 
Jr 2 Jr 2 v c[u) 
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= - / (-(R + 5)20) o (t, x) (t x x Y - txyxx) dXdY 

= -/ -(R + S)(t,x)(/)(t,x)dtdx. 

This proves the first identity in (16, 15ft : the second one is proven in the same way. □ 

We can now define the semigroup mapping St on V, the original set of variables. 
Definition 6.4. For any T > 0, let St '■ T> — > V be defined as 

St = M o St ° L 

Given (uo, Ro, So, /xo, uq) 6 V, let us denote (u, R, S, fi,v)(t) = St(uo, Ro, So, /io, vq). In 
the theorem that follows, we prove that u(t, x) is a weak solution of the nonlinear wave 
equation. However it is not clear if St is a semigroup. Indeed, we have 

St o St> = M o St o L o M o St> ° L- 

By the semigroup property of St, it would follow immediately that St is also a semigroup 
if we had L o M = Id, but this identity does not hold in general. It is the aim of the last 
section to show that St is a semigroup. 

7. Relabeling symmetry 

We consider the set of transformations of the M 2 -plane given by 

(X,Y)»(f(X),g(Y)) 

for any (/, g) S G 2 , where G is the group of diffeomorphisms on the line, see Definition 
13.31 It is a subgroup of the group of diffeomorphisms of M 2 . Such transformations let 
the characteristics lines invariant. Indeed, vertical and horizontal lines, which correspond 
to the characteristics in our new sets of coordinates, remain vertical and horizontal lines 
through this mapping. In this section, we show that the subgroup G 2 plays an essential role 
by exactly capturing the degree of freedom we have introduced when changing coordinates 
and introduced the equivalent system (12 . 13|) . Given / and g in G, the M 2 plane is stretched 
in the X and Y direction by the transformations X ^ f(X) and Y \— > g(Y). The solutions 
of (|2.13p are preserved and we can define an action of G 2 on the set of solutions H. 

Definition 7.1. For any Z £ TL, f and g in G, we define Z £ TL as 
(7.1) Z(X,Y) = Z{f{X),g{Y)). 

The mapping from TL x G 2 to TL given by Z x (/, g) ^ Z defines an action of the group G 2 
on TL and we denote Z = Z ■ (f,g). 

Proof of well-posedness of Definition 1 7. A For any ft, given Z £ TL(ft) and (/, g) £ G 2 , let 
X l = f(Xt), X r = f(X r ), Y x = g(Yi), Y r = g(Y r ) and 

n=[X l ,X r ] x [Y h Y r ]. 
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Let us prove that Z G H(£i). We have 

Z X (X,Y) = f'(X)Z x (f(X),g(Y)), Z Y (X,Y) = g'(Y)Z x (f(X),g(Y)). 

and 

Z XY (X,Y) = f'(X)g'(Y)Z XY (f(X),g(Y)). 
By using the linearity of the mapping F(Z) in (I2.14p . we get 

Z X y = f'g'Z XY (f,g) 

= f'g'F(Z(f,g))(Z x (f,g),Z Y (f,g)) 

= F(Z(f,g))(f'Z x (f,g),g'Z Y (f,g)) (by the linearity of F(Z)) 

= F(Z)(Z X ,Z Y ) 

and Z is a solution of (I2,13p . Since / and g belong to G, there exists 5 > such that 
f'(X) > 5 for a.e. X G M and g'(Y) > 5 for a.e. FgR, see Lemma [Ml We have to check 
that Z fulfills (14.181) . It is not difficult to do so once one has observed that the equalities 
and inequalities in (I4.18P enjoy the required homogeneity properties. For example, we have 

2J x x x = 2f' 2 J x (f,g)x x (f,g) = f' 2 (c(U(f,g))U x (f,g)) 2 = (c(U)U x ) 2 

and 

x X = f'x x (f,g) > 0. 

We will prove that Z fulfills the condition (ii) in Definition 14.131 after we have introduced 
the action of G 2 on Q. □ 

We can define an action on C as follows. This action corresponds to a stretching of the 
curve in the X and Y direction. 

Definition 7.2. Given (^,3^) 6 C, we define {X ,y) £ C such that 

(7.2) X = f- 1 oXoh y = g - 1 yoh 
where h is the re-normalizing function which yields X + y = 2 Id, that is, 

(7.3) (f~ 1 oX + g~ 1 oy)oh = 2Id. 
We denote {X,y) = {X,y) ■ (f,g). 

Proof of wellposedness of Definition \7.2l Let us denote v = f~ l o X + g~ x o y. We have 
u-Id G W 1,00 (M) because / _1 -Id, 5 _1 -Id, X-Id and ^-Id all belong to W 1,0O (R). There 
exists 5 > such that (/ _1 )' > 5 and {g^ 1 )' > 5 a.e. and therefore v > 5{X + y) = 25. 
Hence, by Lemma l3~6l we have that v is invertible so that h exists and h — Id G W l,ao {M). 
One proves then easily that (X ,y) £ C. □ 

We define the action on Q so that it commutes with the • operation and the actions on 
H and C, that is, 

(7.4) (z • r) • 4> = (z ■ </,) . (r • d>) 

for any Z G Tt, T = (X,y) G C and <f> = {f,g) G G 2 . We obtain the following definition. 
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Definition 7.3. For any G = (X, y, Z, V, W) G and f,g€G, we define G = (AT, J, Z, V, W) G 
as follows 

(x,y) = (x,y)-(f, g ) 

and 

(7.5) V(X) = f'(X)V(f(X)) W(Y) = g'(Y)W(g(Y)) 
and 

(7.6) Z = Zoh 

where h is given by (17. 3p . The mapping from Q x G 2 to Q given by Q x (/, 5) 1— ► Q defines 
an action of the group G 2 on Q that we denote 6 = G • (/, g) . 

To check that this definition is well-posed, we have to check that G £ Q. This can be 
done without any special difficulty and we omit the details here. Let us however prove (|7.4p 
in details as we will use it several times in the following. For any Z £ Ti., T = (X, y) £ C 
and </> = (f,g) £ G 2 , we denote G =_(X,y,Z,V,W) =Z»T,@ = (X,_y, Z, V, W) =@(f), 
f = (X,y) = r • 4> and Q =_ Z • f . We want to prove that Q = Q. We have Z(s) = 
Z(X{s),y(s)) and, by JESJ. = M*)> ^ ° &(*)). Hence, by fZ2), 

J = Z(#, 3>) = Z(/ o X, g o y) = Z(X oh,y oh) = z. 

We have V(X(s)) = Z x (X(s),y(s)) and, by 

V(*(a)) = /'(*(a))V(/ o #(*)) = /'(^( a ))V(Af o 

= /'(*(a))Zx(#o/i( s ),;yo/i( s )) 

and 

V(*(a)) = Zx(*(a), = f'(X(s))Z x (f o *( s ), 5 o y( s )) 

= f'(X(s))Z x (Xoh(s),yoh(s)). 

Hence, V = V. Similarly one proves that W = VV, which concludes the proof of (|7.4|l . 

i?nd of proof of well-posedness of Definition \7.1\ For any (p £ G 2 and Z G 7i, it remains to 
prove that Z, as defined by (|7.1j) . fulfills the condition (ii) in Definition 14.131 By this same 
condition, for any Z £ 7i, there exists a curve T = (X,y) G Q such that Z »Y £ Q. From 
(|7.4p . it follows that, for the curve T = V ■ (p, we have Z • T G <5 because (Z »T) ■ cp £ Q and 
therefore Z fulfills the condition (ii) in Definition 14.131 □ 

Definition 7.4. For any tp = (ipi,ip 2 ) G T and f,g £ G, we define ip = (^1,^2) £ J 7 as 
follows 

X!(X) = Xl (f(X)), U^X) = UxtfiX)), MX) = Ji(/(X)), 

x 2 (X) = x 2 (g(Y)), U 2 (Y) = U 2 (g(Y)), J 2 (Y) = J 2 (g(Y)), 

and 

Vt(X) = Vi(/(X))/'(X), V 2 (Y) = V 2 (g(Y))g'(Y). 



THE NONLINEAR VARIATIONAL WAVE EQUATION 63 

The mapping from T x G 2 to J- given by ifi x (/, g) i— > ih defines an action of the group G 2 
on J- , and we denote 

if) = if> ■ (f>. 



Proof of well-posedness of Definition \7.4\ We have to check that if> = (ipi,ip2) G J~ '■ We 
will only check that the identities (I3.19P in the definition 13.41 of T are fulfilled, as the 
other properties can be checked without difficulty. For any curve (X,y) G C such that 
xi{X) = x 2 (y), let (X,y) = (X,y) ■ (/,<?). We have 

= x l of- 1 oX( s ) = x 1 oXoh- 1 {s) = x 2 oyoh- 1 (s) =x 2 of- 1 oy( s ) = X2 (y(s)) 

and therefore, since if) G T , Ui(X(s)) = U 2 (y(s)) for all s£l, which implies 

u l (x) = UioXoh = u 2 oyoh = u 2 (y) 

and this proves (I3.19aj) for if). Similarly, one proves that (I3.19bll holds for V- □ 



In the following lemma, we show that all the mappings given in (|5.9|) are equivariant 
with respect to the action of the group G 2 . 

Lemma 7.5. The mappings E, tx, C, D and C are G 2 -equivariant, that is, for all <f> = 

(/,</)« c 1 . 

(7.7a) E(Z • 0) = E(Z) ■ (f>, 

(7.7b) t T (Z • <f>) = t T (Z) ■ (f> 

for all Z £TL and 

(7.7c) S(0 • 4>) = S(0) • <f> 

for all G and 

(7.7d) D(0 • 0) = D(0) • 

/or a// G Go and 

(7.7e) C(V> • <f>) = C(iP) ■ 

for if) G T. Therefore St is G 2 -equivariant, that is, 
(7.8) S T (if> • 0) = 5 T (^) • </> 

/or a// ip £ J 7 . 

Proof. Let us prove (1773 . We denote = (X,y,Z,V,W) = E(Z), Z = Z • 0, = 
(Af, 3>, Z, V, W) = E(Z) and = (X, y, Z, V, W) = ■ 0. We want to prove that = 0. 
First we prove that T = T where T = (X,y) and f = (X,y). By definition (15, lj) . we have 

X(s) = sup{X G K | t(/(X'), 5(1" )) < for all X' < X and F' such that X' + Y' = s}. 

We have 

t(f(X(s)),g(y(s))) = t(X o h(s),y o h(s)) = 0, 
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by the definition of {X,y) given by (|5.1I) which implies that t(X(s),y(s)) = for all sGK, 
Hence, 

(7.9) X < X. 

Assume that X(s) < X(s) for some point s. We have t(f(X(s)),g(y(s))) = and due to 
the monotonicity of t, it implies that t(X, Y) = for all (X, Y) G [/ o X(s), f o ^ (s)] x[go 
y( s )>9 ° 3^( s )]- If / ° -^( s ) < 2/i(s) - go X(s), we obtain a contradiction. Indeed, if we set 
X' = 2h{s) — go y{s) and Y' = g o 3^(s), then we have 

X' < 2/i(s) - 5 o = 2/i(s) -yo h(s) = X o 

so that i(X', Y') = and X' + Y' = 2h(s), which contradicts the definition (EH]) of (X,y) 
at /i(s). If / o X(s) > 2h(s) -go X{s), then we set X' = f o X(s) <Xo h(s) = f o #( s ) 
and Y' = 2h(s) - f o X{s). We have t(X', Y') = and X' + Y' = 27»(a), which also leads 
to a contradiction of (|5.ip . Hence, we have proved that X = X and therefore F = F. It 
means that 

e(z • 0) = (z • 0) • F = (z ■ <j>) • f = (z • 0) • (r • 0). 

Hence, by (17. 4p . it yields 

E(Z ■ </>) = (Z • r) • = E(Z) • 0, 

and we have proved (I7.7ap . The identity (I7.7bl) follows directly from the definition of ty. 
Let us prove (lT77cj) . For any 9 = (X,y,Z, V,W), we denote Z = 5(G), Z = 5(6 • </>), 
Z = Z ■ (f). We want to prove that Z = Z. By definition of the solution operator 5, we have 

z • (r • 4>) = e • </>, 

and 

z • (r • 0) = (z • 0) • (r • </>) = (z • r) • <f) = e • <j>, 

by (|7.4I) . Hence, Z and Z are solutions that match the same data on a curve. Since the 
solution is unique by Theorem 14.151 we get Z = Z. The property (|7.7d|) follows directly 
from the definitions. Let us prove (I7.7el) . For any ip = (ipi,tp2) G F ■, <P G G 2 , we denote 

$ = (^1,^2) = v • 0, e = (*,;y,f ,v,yv) = c(v), e = (*,y,z,v,w) =_c(^and 

8 = • ^. We want to prove that 6 = 0. First, we prove that F = T where F = (X,y) 
and F = (X,y). By definition (|5.ip . we have 

(7.10) X(s) = sup{X G R I n o /(X') < x 2 o 5 (Y'), 

for all X' < X and Y' such that X' + Y' = s}, 

and 

X{s) = sup{X G K I x x (X') < x 2 (Y'),for all X < X and Y"' such that X' + Y' = s}. 

By continuity, we have x\(X(s)) = X2(y{s)) so that xi o X o h = £2 o y o h. Hence, 
xi o f o X = X2 o g o y and it implies, by (|7.10|) . that 

X < X. 
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The proof then resembles to what was done above after (|7.9I) , Let us assume that X(s) < 
X(s) for some point s. We have f{X{s)) < f(X{s)) and g(Y{s)) < g{y{s)). By using the 
monotonicity of x\ and x 2 , we get 

x\ o / o X(s) < x\ o / o X(s) = X2 o g o y(s) < X2 o g o 3^(s) = X\ o / o 

Hence, x\o f a X(s) = x\ o f o A'(s) and x 2 o g o 3^(s) =i2°ff° 3^( s )- Since »i and X2 are 
decreasing, it follows that xi and X2 are constant on [foX(s), foX(s)] and [<7°3^(s), <7°CP(s)], 
respectively. If / o X(s) < 2h(s) —go X(s), then we obtain a contradiction. Indeed, let us 
set X' = 2h(s) — go y(s) and Y' = go y{s), we have 

X' < 2h(s) - go y( s ) = 2h-yo h(s) = Xo h(s) 

so that xi{X') = x 2 (Y') = and X' + Y' = 2h(s), which contradicts the definition (1X22]) 
of (X, y) at h(s). If / o X(s) > 2h(s) - go X(s), then we have 

Xl (f o X(s)) = x 2 (2h(s) - f o X(s)), 

which, in the same way, leads to a contradiction of (|3.22|) . Thus we have proved that T = f . 
We then prove that Z = Z, V = V and VV = VV. It is just a matter of applying directly 
the definitions. For example, we have 

U(s) = U x o X(s) = U 1 ofo X{s) = U 1 ofo X{s) = U x oXo h(s) = Uo h(s) = U(s) 

and 

v (X) - f(x)v (f(x)) - fWdfW) - - _ v (x) 

v 1 (x) - f pov l( /(x)) - 2c{Ui{f{x))) - - ^^yy - MX). 

The equivariance property (|7.8j) of St follows directly from the definition of St and the 
equivariance properties (|7.7p . □ 



Definition 7.6. We define by T ' /G 2 the quotient of T with respect to the action of the 
group G 2 on T , that is, 

if) ~ ?/> if there exists 4> S G 2 such that ip = ip ■ (p. 

Definition 7.7. Let 

Fo = = {ifruifo) 6 T I x\ + Ji = Id and x 2 + J 2 = Id} 

and II: T — ► T§ be the projection on given by ip = (tpi,ip 2 ) = Il(V') where ip E J-q is 
defined as follows. Let 

(7.11) f(X) = x 1 (X) + J 1 (X), g(Y)=x 2 (Y) + J 2 (Y), 

we set 

tp = if] ■ ifr 1 . 

Lemma 7.8. The following statements hold: 
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(i) We have 

(7.12) tp~tp if and only if U(tp) = U(tp) 
so i/iai the sets T j G 2 and J-q are in bijection. 

(ii) We have 

(7.13) M o n = M 
and 

(7.14) LoM|^ =Id|^ and M o L = Id 

so i/iaf i/ie seis 2?, .Fo anc ^ J~ I G 2 are in bijection. 

(iii) We have 

(7.15) noS T on = no5 T . 

Note that the first identity in (17.14p is equivalent to 

(7.16) LoMon = n. 

Proof. Step 1. We prove (|7.12j) . If tp ~ tp, then there exists </> E G 2 such that tp = tp ■ <p. 
Let <p = (/, g) and (p = (/, g) be given by (17. lip for tp and tp, respectively. One can check 
that <p = po p and therefore 

n($ = ^ • (0)- 1 = (tp • 0) • (0 o = ^ • o o = ^ • = n(v). 

Conversely, if H(tp) = n(^) then f/S • p~ l = tp ■ (ft -1 so that tp = (tp ■ <p~ l ) ■ <p = ip ■ (<p~ l o cp) 
and V> and -0 are equivalent. 

Step 2. We prove that Lo M = Idjr . Given tp = (tpx, tp2) E Fo, let us consider 
(u, R, S, fj,,v) = li(tpi,tp2) and $ = (tp\,tp2) = M(u, R, S, /J,,v). We want to prove that 
tp = tp. Let 

(7.17) g(x) = sup{X eR\ xi(X) < x}. 

It is not hard to prove, using the fact that X\ is increasing and continuous, that 

(7.18) xi{g(x))=x 

for all x E R and x|f 1 ((— oo, x)) = (— oo,g(x)). For any x E R, we have, by (I6.1cj) . that 

(7.19) K(~oo,x)) = I J[(X)dX= [ 9{X J[(X)dX = Ji(g{x)) 

J x-^ 1 {(—oo,x)) J —oo 

because J\(— oo) = 0. Since tp E Fo, x\ + J\ = Id and we get, by (I7.18P and (|7.19j) . that 

(7.20) n((-oo,x)) + x = g(x). 
From the definition of x\, we then obtain that 

(7.21) x\(X) = sup{x E R | g(x) < X}. 

For any given X E R, let us consider an increasing sequence z% tending to x\(X) such 
that g(zi) < X; such sequence exists by (|7.21|) . Since x\ is increasing and using (17.18(1 . it 
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follows that Zi < x\(X). Letting i tend to oo, we obtain x\(X) < x\(X). Assume that 
x\{X) < x\(X). Then, there exists x such that x\(X) < x < x\(X) and (I7,21|) then implies 
that g(x) > X. On the other hand, x = x\(g(x)) < x±(X) implies g(x) < X because x\ 
is increasing, which gives us a contradiction. Hence, we have x\ = x±. It follows directly 
from the definitions, since x\ + J\ = Id, that J\ = J\ and XJ\ = U\. It follows from the 
definition t|3.34e|) and (16. 2|) that V\ = V\ and V% = V% . Thus we have proved that tpi = ipl- 
In the same way, we prove that ^2 = 1P2, which concludes the proof that L o M = Idjr . 

Step 3. We prove that M o L = Id. Given (u,R,S,n,v) G T>, let V = (^1,^2) = 
L(u, i?, 5, /z, f ) and (u, i?, 5*, /2, v) = M(^). We want to prove that (u, R, S, fi, v) = (n, R, S, fi, v). 
Let g be the function defined as before by (|7.17|) . The same computation that leads to (|7.20p 
now gives 

(7.22) p,((-oo,x)) +x = g{x). 

Given I £ 1, we consider an increasing sequence Xi which converges to x\{X) and such 
that 00, x^) + Xi < X. Passing to the limit and since x 1— ► /u((— 00, x)) is lower 
semi-continuous, we obtain fi((— 00, x±(X))) + X\(X) < X. We take X = g(x) and get 

(7.23) n((-oo,x)) +x <g{x). 

From the definition of <?, there exists an increasing sequence Xi which converges to g(x) 
such that xi(Xi) < x. The definition (I3.34ah of x\ tells us that 00,2)) + x > Xi. 
Letting i tend to infinity, we obtain fi((— 00, x)) +x > g(x) which, together with (|7.23l) . 
yields 

(7.24) fi({-oo,x)) + x = g(x). 

Comparing (17.241) and (|7.22|) we get that p, = fj,. Similarly, one proves that v = v. It is clear 
from the definitions that u = u. The fact that R = R and S = S follow from (|3.34e|) and 
(|6.2p . Hence, we have proved that (u,R,S,p,,P) = (u,R,S,ii,v) and M o L = Idx>. Step 
4. We prove (|7.15p . For any ip = (^i,^) G we denote ipr = St^- Let 4> = (f,g) G G 2 
and (f>T = (/t, 5t) G G 2 be defined as in (17. lip so that II^ = ip ■ _1 and TlipT = ipT • '/'t 1 - 
By using (|7.8I) . we get 

s T o n(v>) = s T (^ • cp' 1 ) = s T {^) ■ (\r x 

and therefore St n(^) and St($) are equivalent. Then, (I7.15P follows from (|7.12l) . □ 

We now come to our main theorem. 

Theorem 7.9. Given (uo, Ro, So, fj,o, v$) G V, let us denote (u, R, S, /i, v)(t) = St(uo, Ro, So, no, Uq). 
Then u is a weak solution of the nonlinear variational wave equation (|1.1|) . that is, 

(7.25) / ((j> t - (c(u)(j)) x )Rdxdt + (fa + (c(u)<j)) x )S dxdt = 
Jr 2 Jm 2 

for all smooth functions (j) with compact support and where 



(7.26) 



R = ut + c(u)u x , S = ut — c(u)u : 
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Moreover, the measures ji{t) and v(t) satisfy the following equations in the sense of distri- 
bution 

(7.27a) + v) t - (c(fi - v)) x = 
and 

(7.27b) {- c {n-v)) t -{n + v) x = Q. 
The mapping St '■ T) — > T> is a semigroup, that is, 

St+t' = St o s t > 

for all positive t and t! . 

Proof. From Lemma 16.31 we know that (|7.26l) is fulfilled. On can check that (I7.25P is 
equivalent to 

(7.28) Rt - (c(u)R) x + S t + (c(u)S) x + c'(u) = 

lc[u) 

in the sense of distributions, which makes senses, as R, S belong to L? and c(u), c'(u) are 
bounded. After a change of variables, we have 

(R<j) t -(c(u)R)4 )x )(t,x)dtdx= I {R(fa-c{u)^ x )){t,x)(t x x Y -t Y x x )dXdY 

R 2 JR 2 



c(u)R((f)t — c(u)4> x ))(t,x)txXY dXdY 
2 / c(U)U x <t>(t,x) Y dXdY 



by (|6.14aP and because 4>(t,x) Y = <f>t{t, x)t Y + 4> x {t, x)x Y = —^r{4>t — c4> x )(t,x). We 
integrate by parts and obtain 



{R<t> t - (c(u)R)(j) x )(t, x) dtdx = 2 / (c(U)U x )Y<t>(t, x) dXdY 

(7.29) = 2 [ (c(U)Uxy + c'(U)U x U Y )<t>{t,x)dXdY. 
In the same way, one proves that 

(7.30) / (Sfo + (c(u)S)<j) x )(t,x)dtdx = 2 [ (c(U)U X y + c'(U)U x U Y )<t>(t,x) dXdY. 

JM 2 JM 2 

We have, after a change of variables, 
(7.31) 

R2 - 2RS + S2 -c'(u)6dtdx = 2 /" ( r2 ~ 2RS + s2 j 



( 2c(u) 2 c '( n )^) (*> x ) x xxy dXdY. 



2c(u) J R 2\ 2c{uf 

We introduce the set A = A\ U A2 where 

(7.32) A 1 = {(X, 7) G R 2 x x (X, Y) = 0, x Y {X, Y) > and c'(U)(X, Y) / 0} 
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and 

(7.33) A 2 = {(X,Y) G R 2 | x Y (X,Y) = 0, x x (X, Y) > and c(U)(X, Y) ^ 0}. 
We claim that 

(7.34) meas(A) = meas(^4i) = meas(j42) = 0. 

We prove this claim later. By using (|7.34h . we get 

R 2 — IRS + S 2 , , . , , , _ f / R? — IRS + S 2 , , . , \ . . ,„ 

— c!{u)cj)dtdx = 2 / — — </{u)<t>)(t,x)xxx Y dXdY 

2c[u) J a c ^ 2c(u) z / 

= / f— ZY + 2£/x£/y + — xx)c'(UU(t,x) dXdY 
J A c \x x x Y J 

(7-35) = (2^ + 2CWy + ^)c>(U)<f>(t,x)dXdY 

(7.36) = / + 2U X U Y + ^-)c'(U)cl>(t,x)dXdY. 

Note that (17.341) is necessary to get (I7.36P from (|7.35l) as the integrand in (|7.35l) does not 
vanish on A. After combining (I7.29P , (|7.30p and (|7.36j) . and using the governing equations 
(l2~T3l) . we get 

/ (R + S)4> t - (c(u)(R - S))<p x - c'(u) ( R ~ S / 4>dtdx 

= - [ (ic(U)UxY -%{Jxx Y + Jyx x ) +2c , (U)U Y Ux)<t>(t,x)dXdY 
= 0, 

which proves (|7.28p and therefore (I7.25P holds. It remains to prove the claim (|7.34h . Let 
us introduce the set 

Ax{X) = {Y G R | (X,Y) G A x }. 

Let us prove that, for almost every X G R, meas(^4i(X)) = and therefore, by Fubini's 
theorem, meas(Ai) = 0. We consider a point Yq G A\{X) and a rectangle = [JQ, X r ] x 
[y";,y r ] which contains (X, Yq). Since Z G T~l[p), there exist 5 > such that (xx + 
Jx)(A",y) > 5 for almost every X G R and all Y G R. Since x x G Wy'°°(0), the function 
xx is continuous with respect to Y for almost any given X G R. Formally the argument 
goes as follows: We consider a fixed given X G R and denote f(Y) = x x (X,Y). For any 
Yq G Ai(X), we have, by definition, xx(X,Yq) = /(Yq) = 0. By using (I2.13bp . we get 

f'(Y )=xxY(X,Y) = 

because xx = implies Ux = 0, see (I4.18cp . We do not have enough regularity to differen- 
tiate (l2.13bP ; but if nevertheless we do so, then we formally obtain 

f"(Y ) = x X yy(X,Yq) = ^{U Y x XY + Uxyxy)(X,Y ) 
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= ^(JXX 2 Y )(X,Y ) 

where we have used again the fact xx(X,Y ) = Ux(X,Y ) = 0. We have Jx(X, Y ) = 
(x x + Jx)(X,Y ) > S and x Y (X,Y ) > 0, c' 2 (U(X, Y„)) > because (X,Y ) G A r {X). 
Hence, /"(Yq) > and it implies that f(Y) > for all Y different from Yq in a neighborhood 
of Yo, so that the points in A\{X) are isolated. Let us now prove this result rigorously. 
Again, we consider Yq S A\ (X) and a rectangle ft = [Xi,X r ] x [Y],YV| which contains 
(X,Yq). Without loss of generality, we assume that Y is a Lebesgue point for the function 
Y i— > xy(X,Y) and, therefore, since xy(X,Yq) > 0, there exists S > such that 

rY 



(7.37) 



x Y {Y) dY > 5'{Y -Y ) 



Yo 



in a neighborhood of Yo- We can choose 5' > such that, in addition, 

c'(U(X,Y)) > 5 and J X (X,Y) > 5 

in a neighborhood of Yo (we recall that U is continuous). We have, after using the governing 
equations (12 . 13|) . 



xx 



Y 



Y d 

— {U Y xx + UxXY)dY 
2c 



Y d 



Yo 



Y rY 

— {U Y / x XY dY + XY I U X YdY)dY 
Yo 
Y 



2c 



Yo 



Y d 



Y d d ~\ 

x X y dY + XY I {t^{x y Jx + Jyxx ~ 7^U Y U x ))dY) dY. 



IY ^ v JYo JYo 
Since xxy and Uxy are bounded (by (|2.13p ) and xx{X,Yq) = Ux(X,Yq) = 0, we have 
that x x (X,Y) < C\Y-Y \ and U X (X,Y) < C\Y-Y \ in a neighborhood of Y for a 
constant C which depends only on Hence, 









/ x X Y dY 




f 


JYo 




JY 



2c 



< C 

JYo 

<C(Y-Y ) 2 . 



(Uyxx + Uxxy) dY 



dY 







r 


Y-Y 


'Yo 





Y d 



Thus, 
(7.38) 

In the same way, one proves that 



Y 



it- (U Y / x XY dY)dY 
Yo 2c JYo 



<C\Y-Y n 



Y d 



Yo 2c XY Jyo \^ {JYXX 



2c 



U Y U x ))dY)dY 



<C\Y-Y n 
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For Y > Yq, we have 



= —vk XYdY ) 



>^(Y-Y^ 



in a neighborhood of Yo- We can check that the same inequality holds for Y < Yq. Finally, 
we obtain that, in a neighborhood of Yo, 

(7.39) x x (X,Y) > ^(Y -Y ) 2 -C\Y -Y \ 3 > ^-(Y -Y ) 2 . 

To complete the argument, we consider the sets 

A\{X) = {Y G Ai(X) n [Y h Y r ] | xx(X,Y Q ) = 

and x x (X, Y)>0 for all Y G [Y - ~, Y + ~] \ {y }} 

for any integer A:. By (|7.39|) . we have 

A 1 (X)n[Y l ,Y r } = U k>0 A k 1 (X). 

At the same time, since A\(X) consists of points separated by a distance of at least \, we 
have m.eas(A\(X)) = 0. Hence, after taking sequences of Y[ and Y r which tend to plus and 
minus infinity, respectively, we get meas(Ai(X)) = so that meas(j4i) = and the proof 
of the claim (|7.34j) is complete. Let us prove (|7.27aj) . that is, 



(fit — c(u)(fi x ) dfidt + (4>t + c(u)(fi x ) dvdt = 
Jm? 

for all smooth function <f> with compact support. We have, after a change of variables, that 
it - c(u)<fi x ) dfj,(tfjdt 

t(t, x(t, a)) - c(u(t, x(t, s)))(fi x (t, x(t, s)))V 4 (t, X(t, s))X s (t, s) ds)dt, 



where we have added the dependence in t of the values of 0(t) = L(it, R, S, u, fi, f)(t) (which 
gives x(t, s), Vi(t, X), u(t, s) and X(t, s) in the equation above). We proceed to the change 
of variables (X, Y) ^ (t(X, Y), s = ^(X + Y)) whose Jacobian is equal to x ^ e "^ y and get 



(7.40) / ( f {<t> t -c{u)<f> x )dn(t))dt 
Jib ^ Jib / 



i t (i, x) - c(u(t, x))<fi x {t, x))J x (X, Y)X s (t, S )^2L±^0^I1 dXd Y. 

2c{U(X,Y)) 
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Since t(X(t, s), y(t, s)) = t, by definition, we get txX s + tyy s = and, since X(s)+y(s) 
2s, we have X s + y s = 2. Hence, (xx + xy)X s (t, s) = 2xy and (|7.40|) implies 



b t - c(u)(f> x ) d/J,(t))dt 

(&(t, x) - c(u{t, x))<j> x (t, x))Jx(X, Y) **\\ dXdY - 
! c(U{X,Y)) 

Since 4>(t,x) Y = ~^{<t>t ~ c(u)<f> x )(t,x), it yields 

b t - c(u)(f> x )dn(tj)dt = - / <t>(t,x) Y JxdXdY. 
Similarly, one proves that 

/ ( I {<t>t + c(u)<t> x )dv{t))dt= I <P{t,x)xJYdXdY 

so that 
(7.41) 

bt — c(u)(j) x ) dfidt+ / (4>t + c(u)cj) x ) dvdt = / (— <j>(t, x)y Jx+<fi(t, x)x Jy) dXdY = 0, 




by integration by parts, as the support of (j) is compact. Similarly one proves (|7.27bj) . 
Note that the integrand in (I7.4ip is equal to the exact form d{4>dJ) and equation (|7.27aP 
is actually equivalent to ddJ = while (|7.27bjl is equivalent to ddK = 0. The proof of the 
semigroup property follows in a straightforward manner from the results that have been 
established in this section. We have 

St ° St 1 = M o St ° L o M o St 1 ° L 

= MonoS T oLoMonoS T 'oL by fl7TT3jl 

= MoIIo5toIIo5t'oL by d7TT6]l 

= M o n o S T o S t > o L by (1715]) 

= Mo S T o S T > oL by i7l3| 

= M o St+t 1 ° L by Theorem 15.51 

= St+t 1 ■ 

□ 

The semigroup of solution we have constructed is conservative in the sense given by the 
following theorem. 

Theorem 7.10. Given (uo, Ro, So, fio, uq) G V, let us denote (u, R, S, ^, v){t) = St(uo, Ro, Sq, no, vo) 
We have 

(i) For allt£R 
(7.42) /i(t)(R) + u(t)(R) = /io(M) + M 
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(ii) For almost every i£l, the singular part of n{t) and v(t) are concentrated on the 
set where c'(u) = 0. 

This theorem corresponds to Theorem 3 in [5]. We use a different proof based on the 
coarea formula. 

Proof. Let us prove (i). We consider a given time that we denote r (to avoid any con- 
fusion with the function t(X,Y)). As in the proof of the previous theorem, we add the 
dependence in time of the values of 0(r) = ~L(u, R, S,u, [i, v)(t). In particular, the curve 
(X(t, s),y(r, s)) corresponds to the curve where time is constant and equal to r, that is, 
t((X(r, s),y(r, s))) = t. By definition (see (I6.3bl) . and Definition I5.ip . we have, for any 
Borel set B, 



(7.43) n T (B)= [ Vi{X{ T ,s))X(T,s)ds 

(7.44) = f J x (X(T,s),y(T,s))X(T,s)ds. 

J{seR | x(x(T,s),y(r,8))eB} 

Here, the notation may be confusing as x denotes two different but of course very related 
functions. In (|7.43|) . we have x(t, s) = Z 2 (r, s), which corresponds to the space variable Z 2 
parametrized by s at time r while, in (|7.44j) . x(X,Y) = Z%(X,Y), corresponds to the value 
of the space variable Z 2 , where Z(X,Y) is the solution of (|2.13l) on the whole M 2 plane. 
We have Z 2 (t,s) = Z 2 (X(t, s), X(t, s)) by (14. 15ft abd Definition 15.11 Correspondingly, we 
have 



u T (B) = / J Y (X(T,s),y(r,s))y(T,s)ds. 
J{seR | x(x( T ,s),y(T,s))eB} 

Hence, 

Hr(R) + u T (R) = [ J x (X(T,s),y(r,s))X(T,s) + J Y (X(r,s),y(r,s))y(T,s)ds 
= lim J(X(T,s),y(r,s))ds 

s— >oo 

= lim J(X(0,s), y(0, s))ds, by Lemma OH 

s— >oo 

= ^o(R) + u (R) 

Let us prove (ii). Let fi T = (// r )ac + (/^r)sing be the Radon-Nykodin decomposition of /x T . 
We want to prove that, for allmost every time r 6 K, we have 

(7.45) (/v)sing({* G M | c'(u(r, x)) + 0}) = 0. 

Let us introduce 

4 = {s£R| x x (X{r,s),y(T,s)) > 0}. 

The set A T corresponds to A in (I6.11|) in the proof of Lemma I6.2L In this same proof, we 
obtain that, for any Borel set B, 

Gu T ) ac (B)=Mr(£n(s(T, 04 T ) C )) C ) 
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so that 

(/V)sing(£) = Ht{B n (X(T, (^ T ) C ))), 

because meas(x(r, (A T ) C )) = 0. Hence, 

(7.46) (MrWOB) = / Jx^sjJtr.s))^^)*. 

j{seK | x(X{T,s),y(T,s))&Bn{A T ) c } 

We introduce the set 

E = {(X,Y) G R 2 | x x (X,Y) = and c'{U{X,Y)) / 0}. 
By using (|7.46|) . we get 

(7.47) (J, T ({x G R | c'(u(t, x)) ^ 0}) = / Jx(*(r, s),3>(r, a))*(s) 

J{sgir | (X( T ,s),y(T,s))eE} 

By the coarea formula, see [I], we get 

/ H 1 (Ent- 1 (T))dr= [ Jt\ + 1\ dXdY = [ -^-dXdY = 
Jr Je v Je c {U) 

because of (17341 . Here, H l denotes the one- dimensional Hausdorff measure. Hence, we 
have that, for allmost every time r G R, the set E n i (t) has zero one-dimensional 
Hausdorff measure. We claim that, if fj, T ({x G R | c'(u(r,x)) ^ 0}) > 0, then H X {E n 
f (t)) > 0. Indeed, let us define, for a given r, the mapping r T : s i— ► (<^(t, s), 3^(r, s)) 
from 1 to I 2 . We rewrite (|7.47l) as 

I^ t ({xGR\c'(u(t,x))^0})= [ J x (X(T,s),y(T,s))X(s)ds. 

JTr\E) 

In particular it implies that, if n T {{x G R | c'(n(r, x)) 7^ 0}) > 0, then meas(r~ 1 (i?)) > 0. 
By the area formula, we have 

H\E) > H\T T o r- 1 ^)) = / (X 2 + y 2 ,) 1 ' 2 ds > meas^; 1 ^)) 

because (X 2 + y 2 ) 1 ' 2 > \{X S + y s ) = 1. Hence, our claim is proved and it follows that 
fj- T ({x G R I c'(u(t, x)) 7^ 0}) > for at most almost every r G R and we have proved 

4231]). □ 

Theorem 7.11 (Finite speed of propagation). Given t > and x G R, for any two 
initial datas Co = {uq, Rq, So, uq, ^0) an d Co = (uo, Rq, So, £lq,vq) in T>, we consider the 
corresponding solutions (u, R, S, fi, v)(t) and (u, R, S, fi, v){t) given by Definition \6.4\ If the 
restrictions of Co an d Co are equal on [x — ret,x + ret], then the two solutions coincide at 
(t, x), that is, u(t, x) = u(t, x) . 

Proof. For a given Co = (uo, Ro, So, no, vo), we define Co equal to Co on [x — ret, x + ret] and 
zero otherwise, i.e., 



(u ,Ro,S ){x) 



(uo, Ro, So)(x) if x G [x — ret, x + ret] 
otherwise 
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Figure 6. The semigroup S t . 



and 

fio(E) = /j,o(E n [x — ret, x + ret]), vo(E) = Uq(E Pi [x — ret, x + ret]) 

for any Borel set -E. It is enough to prove that the theorem holds for this particular Co- We 
have to compute the solutions for (q and Co ■ Let us denote xi = x — ret , x r = x + ret. We 
set if} = C(Co) and ip = C(Co). 

Step 1. We want to compute V as a function of ip. We denote X; = xi, Y\ = xi, 

[Xi,X r ] x [yj,y r ]. Let us prove that 



X r = x r + [io([xi,x r ]), Y r = x r + ^o([^ 5 ^r]) and f2 



(7.48) 



and 



(7.49) 



xi(X) 



x 2 (Y) 



X 

xi(X + no(-oo,xi)) 
X - fi ([xi,x r }) 



Y 

x 2 (Y + v (-oo,xi)) 
Y - i> ([xi,x r ]) 



ifX<X t 

if Xi < X < X r 

ifX r <X 



if Y < Y { 

if Yi < Y < Y r 

if Y r < Y. 



From the definition (I3.34al) . we have 
(7.50) X1 (X) = sup{x' e 



x + /ig(— oo, x 



<X} 



First case: X < x\. For any x' such that x' + Jlo(—oo,x') < X we have x' < X. Hence, 
x' < xi and jlo(—oo,x') = 0. It follows that x\{X) = X. Second case: X\ < X < 
X r . For any x' such that x' + po(—oo,x') < X, we have x' < x r . Let us assume the 
opposite, i.e., x' > x r , then jio(—oo, x') = /io((— oo, x')f] [xi, x r \) = fJ*o([xi, x r ]) and therefore 
x' + fJLo([xi, x r ]) < X < x r + A*o([^/ ; a?r])) which gives a contradiction. We can assume 
without loss of generality that x' > x\ because for x' = xi, fio(—oo,x') = and we have 
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x' + p,o(—oo,x') = x' = xi < X . Thus we have x' G [x/,x r ] and (17.50P rewrites 

(7.51) x\(X) = sup{x' G [a?j,x r ] | a;' + /io[xi,x') < X}. 
We now want to prove that, for Xi < X < X r , 

(7.52) xi(X + /x (— oo, Xi)) = sup{x' G [xj, x r ] | x' + ,Uo[xj, x') < X}. 

For any x' such that x' + hq{— oo, x') < X + /io(— oo, x/), we have x' < x r . Let us assume the 
opposite, i.e., x' > x r , then x r + {j,q([xi, x r ]) < x' + hq(\xi, x')) < X implies a contradiction 
with the assumption that X < X r . For x' = xi, we have x' + ^o( — oo, x') < X + ^o( — oo,x;) 
so that we can assume without loss of generality that x' > x/. Hence, 

x\(X + /xo(— oo, xj)) = sup{x' G R | x' + fj,o(—oo, x ) < X + /io(— oo, x/)} 

= sup{x' G [xi,x r ] | x + /io(-oo, x') < X + /xo(— oo, x^)} 

and (|7.52l) follows. By comparing (|7.5ip and (|7.52|l , we get xi(X) = x\{X + /xq( — oo,xi)) 
for X; < X < X r . Third case: X r < X. For x' = x r , we have x' + /2q(— oo,x') < X r = 
x r + fj,o[xi,x r ] <X. Hence, xi(X) = sup{x' G [x r ,oo) | x' + fio(— oo,x') < X}. Since, for 
x' > x r , jtZo(— oo,x') = jUo([xi, x r ]), it follows that X\{X) = X — /Uo([x/, x r ]). This concludes 
the proof of proved (|7.48p . One proves in the same way (I7.49p . Let (ft = (f,g) G G 2 where 
f '. X i ^ X ~\~ fJ,o{— oo,xi) and g :Y ^ Y + Vq(— oo,x r ). We denote ip = ip ■ (ft. We have 
proved that 

(7.53) xi(X) = xi(X) for X l < X < X r 
and 

(7.54) x 2 {Y) = x 2 (Y) for X t < Y < X r . 

We denote 9 = C(^) and 9 = C(ift). 
Step 2. We prove that 

(7.55) X(s) = X(s) and 3?(s) = y(s) 

for s G [si, s r ] where si = \{X[ + Y{) and s r = \{X r + Y r ). By using the definitions of x\ 
and X2, we obtain that, for any x G R, 

(7.56) xi(x + fJ>o(— oo, x)) = xi(x + /io(— oo, x]) = x 

and that the corresponding statement for X2 holds. Let us now prove that 

(7.57) X{si)=X l and y( Sl ) = Y h 
It follows from (f7T56l) that 

(7.58) xipQ) = x 2 (y 2 ) = Xl 

as we have xi(Xj) = xi(x; + /j, (— oo, xi)) = x\ = x 2 (x; + u (— oo,x/)) = x 2 (X). For an Y 
X < X/, we have xi(X) < xi(X;) = x/. Let us prove that x\(X) < xi{X\). We assume 
the opposite, i.e., that x\(X) = x\{Xi) = x;. Then, there exists an increasing sequence 
Xi such that lim^ooXj = x/ and Xj + Ho{— oo, Xj) < X + ^o( — oo, x;). It implies that 
xi + ^o(— oo, xi) < X + jUo(— oo, x/) because of the lower semicontinuity of x i— > fj, (—oo, x). 
Hence, x\ < X/, which is a contradiction. Thus we have proved that, for any X < X;, 
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xipf) < xx{Xi) = x 2 (Yi) < x 2 (2s - X). Hence, X{s{) = X { and (17371) holds. By using 
similar arguments, one also proves that 

(7.59) X(s r ) = X r , y{s r ) = Y r and x x {X r ) = x 2 (Y r ) = x r . 
and the corresponding results for X and y, that is, 

(7.60) X{si)=X u y(si)=Y u x l {X l )=x 2 {Y l ) = x l 
and 

(7.61) X(s r )=X r , y(s r ) = Y r , x x {X r ) = x 2 (Y r ) = x r . 

In particular, we have proved (|7.55|) for s = si and s = s r . For any s S (si,s r ), either 
X\ < X{s) < X r or Yi < y(s) < Y r . We consider only the case where X{ < X(s) < X r 
as the other case can be treated similarly. By definition of X , there exists an increasing 
sequence Xi such that lim^ooXj = X(s) and x\(Xi) < x 2 (Y,j) where Yi = 2s — X, L . For i 
large enough, we have X[ < X; L < X r and, by (I7.53p . we get 

(7.62) x x {X i )=x l {X i )<x 2 (Y i ) 
If Yi < Y r then x 2 {Yi) = x 2 (Yi) and 

(7.63) x x {Xi) < x 2 (Yi). 

If Yi > Y r then (|7.63|) holds also. Indeed, let us assume the opposite. By the monotonicity 
of x\ and x 2 , we get 

(7.64) Xl (X r ) > Xl {Xi) > x 2 (Yi) > x 2 (Y r ). 

By (|7.59l) . we have x\(X r ) = xi(X(s r )) = x 2 (y(s r )) = x 2 (Y r ) and therefore (I7.64p implies 
that x 2 (Yi) = x 2 (Y r ) = x r . From the definitions of x 2 and x 2 , we know that there exists 
a decreasing sequence Xj such that limj^oo Xj = x 2 (Yi) and xj + z^o(— cxd, xj) > Yi + 
uo(—oo,xi). Letting j tend to infinity, we get x 2 (Yi) + vq{— oo, x 2 {Yi)\ > Yi + vq(— oo,xi). 
Hence, as x 2 (Yi) = x 2 {Y r ) = x r , 

Y r = X r + Vo[xi,X r ] > Yi 

which is a contradiction and we have proved that (|7.63l) holds. If Yj < y(s), we get 

(7.65) *i (*(*)) = xi(X(s)) = x 2 (y(s)) = x 2 (y{s)) 

from ([HHD and (17^91) . If Y t = y(s), x 2 (y(s)) = x x = x 2 (y(s)), by (l739]l and (TTiUl) so 
that (I7T6511 also holds. Then, it follows from (I7T63D and (I7T65D that X{s) = X{s) and the 
proof of (|7.55p is complete. 

Step 3. Let Z = S6» and Z = S6. We prove that 

(7.66) t(X, Y) = t(X, Y), x(X, Y) = x(X, Y), U{X, Y) = U{X, Y) 

for all (X, Y) G $7. Since x\ = x,\ on pQ,X r ] and x 2 = x 2 on [y;,y r ], we get, from the 
definition of L, that 

u 1 = u u Vi = v u j 1 = j 1 + j 1 {x 1 ), k 1 = R 1 + k 1 (x l ) 
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on [Xi, X r ] and 

u 2 = u 2 , V 2 = V 2 , J 2 = J 2 + J 2 (Y l ), k 2 = K 2 + K 2 {Y l ) 

on [Y/,1^.]. Since, by (|7.55j) . the two paths (X,y) and (X,y) in C(Q) are equal, one can 
check, by using the definition of the mapping C, that it implies that 

i(s)=i(s), x(s) = x(s), U(s) = U(s), J(s) = J(s) + J( Sl ), K{s) = K{s) + K{ Sl ) 

for s G [si, s r ] and 

V = V, W = W 

on [Xi,X r ] and [Y;,y r ], respectively. The elements 6 and are equal in $7 except that 
the energy potentials J and K differ up to a constant. However one can check that the 
governing equation (|2.13l) is invariant with respect to addition of a constant to the energy 
potentials. Hence, by the uniqueness result of Lemma [4.101 which holds on finite domains, 
we get (T7T661) . 

Step 4. We prove that there exists (Xq,Yq) G £1 such that 
(7.67) t{X ,Y )=t and x(X ,Y )=x. 



We have 
so that 



xi(Xi) = x 2 (Yi) = x — Kt and x\{X r ) = x 2 (Y r ) = x + nt 



x(Xi,Y[) = xi and x(X r ,Y r ) = x r . 
Let P = (X r ,Yi) denote the right-corner of O. We have 



x{P)-xi= [ T x x (X,Y l )dX= [ T c{U)tx{X,Y{)dX 
Jx x Jx l 

and ^ 

i(P)= [ r i x {X,Yi)dX. 
JXi 

Hence, using the positivity of tx and the assumption that i < c < k, we get 

(7.68) x(P) -x { < Kt(P). 

Similarly, one proves that x r — x{P) < Kt{P), which added to (j7.68J) . yields x\ — x r < 2ni(P) 
or, after plugging the definition of xi and x r , 

(7.69) t < i(P). 

The mapping (X,Y) i— > (i(X, Y), x(X, Y)) is surjective from M 2 to M 2 and there exists 
Po = (Xq,Yq) G M 2 , which may not be unique, such that (|7.67|) is fulfilled. Let us assume 
that 

(7.70) {i(X,Y),x(X,Y))^(t,x) 

for all (X,Y) G Q. By using (|7.69|) and the monotonicity of the function t and x in the 
X and Y directions, we infer that either Xo > X r and < Y or Y < Y\ and Xq < X r . 
We treat only the first the other case can be treated similarly. We have Yq < Y r 
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as, otherwise, x(Xq,Yq) > x(X r ,Y r ) = x T . We introduce the point Pi = (X r ,Yo) 6 fi. 
Let us assume x{P\) > x(Pq) = x. By the monotonicity of x, we get that x{P\) = x(Pq) 
and x x (X,Y ) = for X G [X r ,X ]. It implies that txp^Yb) = for x £ [Xi,X r ] and 
therefore t(P\) = i(Po) = t. However this contradicts the original assumption (I7.70P and 
we must have that x{P\) < x. By following the same type of computation that lead to 
(|7.68j) . we now get 

x > x(P\) = x r + xy(X r , Y) dY > x r — Kt{P\) > x r — nt > x, 

which is a contradiction. Hence, (|7.70l) cannot hold and we have proved (I7.67p . 

Step 5. We now conclude the argument. By definition, we have n(t,x) = U(Xq,Yq) for 
any (X ,Y ) such that (TOTTD holds. By (I7T66D . it follows that U(X ,Y ) = U(X ,Y ) = 
u(t,x) and t(Xo,Yo) = t and x(Xq,Yq) = x. It gives U {f {Xq) , giY^)) = u(t,x) and 
t(f(X ),g(Y ))=t a ndx(f(X ),g(Y ))=x, so that u(t, x) = u(t, x), by ([03]). 




Figure 7. We have t = t(X ,Y ) and x = x(X ,Y ). In the new set of 
coordinates (X,Y), the domain of dependence is given by rectangles. We 
define the points (Xi, Yj) and (X r , Y r ) so that they correspond to the points 
(x — Kt,0) and (x + nt,0). It then follows, from the boundedness of the 
function c(u) that (Xq,Yq) is contained in Q,. 



□ 



8. Examples 

There is a lack of explicit solutions for any choice of c except the trivial case of the linear 
wave equation for which c is constant. We here discuss two examples; first the linear case 
with general initial data, and second, a nonlinear case with very simple initial data. 

8.1. The linear wave equation. In the case of the linear wave equation, the coefficient 
c is constant and the equivalent system (|2.13h rewrites as 

(8.1) Z XY = 0. 
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We consider general initial data (uo, Ro, 5b, /1q, uq) G V, let (V'i, ^2) = L(uo, Ro, Sq, /j,o, vq) 
and 6 = (X,y,Z,V,W) = 0(^1,^2). From JO), we get that 

Z X (X,Y) = V{X) and Z Y (X,Y) = W{X). 

Given a point (X, Y) G R 2 , we denote s = y~ X (Y) and s\ = X ~ X {X) so that 

(8.2) y(s ) = Y and X( Sl ) = X. 

We have 



Z(X,Y) = Z(s ) + I" Z x (X,Y)dX = Z(s )+ [ X V(X)dX 

J X{Y) JX(Y) 



and 

i-Y pY 

Z{X,Y)=Z( S1 )+ Z Y (X,Y)dY = Z( Sl )+ W(Y)dY. 

Jy(x) Jy{x) 

By averaging these two equations, we get 

Z{X,Y) = \{Z{s Q )+Z{ Sl )) + \{ ! V(X)dX+[ W(Y)dY) 
1 1 Jx(Y) Jy(x) 



After a change of variables, it yields 



(8.3) Z(X, Y) = i(Z( So ) + Z( Sl )) + \ £ (V(X(s))X(s) - W(y(s))y(s)) ds. 

We recall that for 9 G G, we have V 2 {X{s))X{s) = W 2 (y(s))y(s), see, for example, (l3~32l . 
For the first component Z\{X,Y) that we denote t(X, Y), we have Z\(s) = t(s) = for all 
sgl and, after using (|3.8b|) . we get 



t(x,Y) = 1 [ S1 (v 2 (x( s ))x( s ) + w 2 (y(s))y( s )) ds 



= ^-(Z 2 ( Sl ) - Z 2 (s )), by ([321) 

(8.4) = ^-(x!(X) - x 2 (y)), by p34b| ). 

2c 

As far as the second component Z 2 (X,Y) = x(X, Y) is concerned, it follows directly from 
JEH) and (13.24hh that 

(8.5) x(X,Y) = ~(xi(X)+x 2 (Y)). 

For the third component Zs(X,Y) = U(X,Y), we have Z^{sq) = uo(xi(X)) and Z^{s\) = 
uq{x 2 {Y)) ■ After using (|3.34e|) and (|3.3ip . we get, after a change of variables, that 

(V 3 (X(s))X(s)= (R (x l (X{s)))x' 1 (X(s))X(s)ds = — / R (x)dx. 



s Jsq 2c J X2 (y) 
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We use the fact that x 1 (X(s 1 )) = x 2 (y(s 1 )) = x 2 (Y), which follows from (ET231) and (1Q) . 
Similarly, we obtain that f£ W(y(s))y(s) ds = -± S (x) dx. Hence, (JHHJ) yields 

1 i p'2(V) 

(8.6) C/(X,y) = -(no(x 1 (X))+n (x 2 (y))) + - / (i? + 5 ) dx. 

From dSaj) and (1831) . it follows that = - ci(X,F) and s 2 00 = + 

ct(X,Y). Therefore, after using (I6.13p . we recover d'Alembert's formula from (18.6(1 . i.e., 

u(t,x) = -(u (x - ct) + u (x + ct)) + — / (i? + 5'o)dx 

2 4c J x _ ct 

for the solution of the linear wave equation. Let us now look at the energy. We use the same 
notation as in the proof of Theorem 17.101 For a given time t, (X(t, s),y(t, s)) denotes the 
curve corresponding to a given time, that is, t(X(t, s),y(t, s)) = t (Beware of the notation, 
i(-,-) denotes a function while t, without argument, denotes a constant). For any point x, 
we have 

fx(t)(-oo,x) = I J x (X(t,s),X{t,s))X s (t,s)ds. 

Jx{X(t,s)y{t,s))<x 

From pT4l) and ([875]) . we get that x{X (t, s), y(t, s)) < x if and only ifxi(X(t,s)) <x + ct. 
Since J X (X,Y) = V 4 (A) = J[(X), we get 



H(t)(-oo,x)= / J[(X(t,s))X s (t,s)ds. 

J x\(X{t,s))<x+ct 

After a change of variables, it yields 

n(t)(—oo,x)= / J[ (X) dX = fio(— oo, x + ct). 

J x 1 (X)<x+ct 

Hence, for any Borel set B, we have 

f jL(t)(B) = f i (B + ct). 

Similarly, we get 

u{t){B) = v Q (B -ct). 

8.2. An example with singular initial data. Let 

(8.7a) u (x) = 1, R (x) = S (x) = 

for all x G M and 

(8.7b) vq = 2/i = 25 

where 5 denotes the Dirac delta function. Our intention is to consider initial data for which 
all the energy is concentrated in a set of zero measure (in this case the origin) and that is 
why we choose uq equal to a constant H Since uq does not belong to L 2 (M), the theory we 



2 If we choose wo = (the only constant in L 2 (R)) then, since c'(0) = 0, one can check from the governing 
equations (|2 . 1 3 that there is no evolution of the solution, and we have that u(t,x) = 0, 2/i(t) = u(t) = 28 
is the conservative solution. 
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have developed does not apply directly. However, we can consider the sequence of solutions 
(un, Rn, Sn, [in^n) given by the semigroup St for the following initial data 



for x G [-N,N] 
otherwise 



and Rg(x) 



yjuj — u yjbj — 0, v$ = 2/io = 25. Given a compact domain in time and space, we 
know that for N large enough the solutions will coincide on this compact domain due to the 
finite time of propagation, see Theorem 17.111 Thus we can define the solution of (jl.ip for 
the initial data (18. 7p as the limit of the solutions u N when N tends to oo. We see that, by 
using the same type of construction, we can actually construct solutions for any initial data 
such that uq,Rq,Sq belong to L^ 0C (M) and fi, v are (not necessarily finite) Radon measures 
(note that, by definition, a Radon measure is finite on compacts). 

For the initial data (uq, Rq, So,fJLo, vo) given by (18.71) . let us denote ("01, ^2) = L(iio, Ro, So, no, vo) 
as defined in Definition 13.81 and = {X , y, Z, V, W) = C^i,^) defined by Definition! 
We first find 



xi{X) 



X 



X 



if X < 0, 
if < X < 1, 
if X > 1, 



Y if Y < 0, 

if < Y < 2, 

Y -2 if Y > 2, 



which yields 



Furthermore 

Ji(X) = 

Ui(X) = 
V X {X) = 

K X {X)-- 

Next, we obtain 



(8.9a) X{s) = < 



{(X,Y) \ X1 (X) = x 2 (Y)} 

{(X, X) I X < 0} U ([0, 1] x [0, 2]) U {(X, X + l) 



X > 1}. 




if X < 0, 
if < X < 1, 
if X > 1, 



if X < 0, 
if < X < 1, 
if X > 1, 



s if s < 0, 

if0<s<l, 

2s - 2 if 1 < s < 3/2, 

s-1/2 if 3/2 < a, 



MY) 

V 2 (Y) 
K 2 (Y) 





Y 

2 



if Y < 0, 
if < Y < 2, 
if Y > 2, 







-y/ c (i) 

-2/c(l) 



s 

2,9 

2 

Ls + 1/2 



if r < 0, 

if < Y < 2, 

if y > 2. 



if 5 < 0, 
if < s < 1, 
if 1 < s < 3/2, 
if 3/2 < s, 
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and 
(8.9b) 

and U(s) 
(8.9c) 



x(s) 



3/2 



if s < 0, 

if < s < 3/2, 

if 3/2 < s, 



1 and 



J(s) 



(8.9d) 

(8.9e) 
and V3 
(8.9f) 

(8.9g) 



if s < 0, 

2s if0<s<3/2, 

3 if 3/2 < a, 



c(l)Vi(X) = V 2 (X) 



c (i)Wi(y) = -w 2 (y) 



K(s) 





__ 2s 

2(s-2) 
c(l) 

1_ 

c(l) 



w. 



and 



c(l)V 5 (X) = V 4 (X) 



c(l)W 5 (Y) = -m(x) 





1/2 


if X < 0, 


< 





if < X < 1 


1/2 


if 1 < X, 


fl/2 


if y < 0, 


< 





if < Y < 2 




1/2 


if 2 < y, 







if X < 0, 


: < 


1 


if < X < 1, 






V 


if 1 < X, 




/ 




if y < 0, 


: < 


1 


if < y < 2, 







if 2 < y 



if s < 0, 
if < s < 1, 
if 1 < s < §, 
if 3/2 < s, 



In the case of the linear wave equation, the solution is explicit. In Figure El we plot the 
curve (X,y) and the curve t(X,Y) = T, for a given T. In Figure [8], the letters A to F 
denote the regions which are delimited by the neighboring solid or dashed black lines. The 
values of Z in these different regions are given in Table CD 

If we consider the choice 

(8.10) c i u ) 2 = P cos 2 u + a sin 2 u 

where a and (3 are strictly positive constants, then no explicit solutions are available for 
the initial data (18. 7p . Due to the finite speed of propagation, we know that 



(8.11) U(X,Y) = 1 on (-oo,0] x (-oo,0] (J [0,1] x [0, 2] (J [ 1, 00) x [2, 00). 

The measures v and \x will become regular measures for t nonzero. We take a = 0.2 and 
(3 = 0.1. The solution is illustrated on Figs. [H I9HT21 Here we have used the numerical 
method described in Section [H 
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t(X,Y) = 
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Figure 8. Plot of the initial data curve and the curve for a given time T. 
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t(X,Y) 





X-Y 
2c 


Y 

2c 


X-Y-l 
2c 


X-l 
2c 


X-Y+l 
2c 


x(X,Y) 





X+Y 
2 


Y 
2 


X+Y-l 
2 


X-l 
2 


X+Y-3 
2 


U(X,Y) 


1 


1 


1 


1 


1 


1 


J(X,Y) 


X + Y 





X 


1 


l + Y 


3 


K(X,Y) 


X-Y 
c 





X 

c 


1 

c 


l-y 

c 


_i 

c 



Table 1. The values of the solution Z = (t, x,U, J, K) of the linear wave 
equation for the initial data given by (|8.7p in the different domains of the 
plane (see Figure [8j the regions A-F are delimited by the dashed and solid 
dark lines). 



9. A NUMERICAL METHOD FOR CONSERVATIVE SOLUTIONS 

Next we describe a general numerical approach to obtain conservative solutions of the 
NVW equation. Traditional finite difference methods will not yield conservative solutions, 
and we here use the full machinery of the analytical approach to derive an efficient numerical 
method for conservative solutions. 



THE NONLINEAR VARIATIONAL WAVE EQUATION 85 

We discretize the problem as follows. Given N, s m ; n and s max > we set h = (s ma _ x — s m i Q )/N 
and S{ = s m ; n -+- ih for i = 0, . . . , N. Let 

Xi = x(si), Yj = y( Sj ), p id = [Xi, Xj] 

for i = 0, . . . , N and j = 0, . . . , N, We compute the solution of (|2.13j) on the domain 
SI = [Xq, Xn] x [Yq, Y/v]- The algorithm follows the same type of iteration as in the proof of 
Lemma f4. 101 and we use the same notation here. We approximate the form Zx(X,Yj) dX 
on the interval by the constant Vjj and the form Zy(Xi,Y) dY on the interval 

[Yj, Yj + i] by the constant Wij. We denote by Z^- and ZV - the approximation of Z on the 
segments Pi-ij — P%,j and Pjj — Pij+i, respectively. The initial curve is approximated on 
the piecewise horizontal and vertical line {J^ 1 ([Pi,i, Pi,-i+i] U [P^j+i, Pj+i^+i]) and we set 

Z^ = Z{ Si ), and V it i = - * / ' V(X) dX for % = 1, . . . , N, 

M — Ji-i-i J Xj_i 

Z v ti = Z( Si ), and Wi ti = - 1 - W(X) dX for i = 0, . . . , N - 1, 

where 2, V and W are given by (18.91) . If Xj — (respectively Yj + i — Yj) is equal to zero, 
then we set V^i (respectively W^j) to zero or an arbitrary value (this value will not have 
any impact on the computed solution) . We compute the solution iteratively on vertical and 
horizontal strips: Given n £ {0, . . . , N}, we assume that the values of 

(9.1a) Zlj, V i: j for 1 < i < n, < j < n, 

(9.1b) Z^, W^j for < i < n, < j < n, 

have been computed. Then, we set iteratively, for j = n + 1, . . . , 1, 

V n+1 ,j^ = V n+lij - (Yj - Yj^)F(±(zZ +l!j + Z^MVn+ij, W ntj ^), 

^n+l.j-l = %n,j-l + (X n +l — Y n )V n+ \ : j, 

W n+1 ,j-x = Wnj-i + (X n+1 - X n )F(~(Z% +ld + Z^MVn+lj, W n)i -i) ; 
and, for i = n + 1, . . . , 2, 

Zj_i jTl+ i = ^j_ 1>n + (ln+1 - Y n )Wi-i t n, 

Vi-X,n+1 = Vi_l,n + (Y t +1 " ^n)i ? (^(^-l ) „ + ^l,n))(^-l,n, Wi_i,„), 
^i-l,n+l = ^i,n+l ~~ (-^t ~~ X-l)^i,n+l, 

Wj_i jn+ i = Wi,n+1 - (-Xi - Xj_i)F(- {Zi <n+1 + ^ n+ i))(V^,n+l, Wi,n+l), 

and 

^0,n+l = ^l,n+l _ (^1 _ X )Vx, n+ i, 
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W , n+ i = W 1>n+1 - (X 1 - X )F(-(Z£ n+1 + Zl n+1 ))(V 1>n+1 , W hn+1 ). 

We have defined the quantities in (|9.1I) for n replaced by n + 1. By induction we have 
computed the solution on the whole domain f2. To compute the solution at a given time 
T, we have to extract a curve (X,y) £ C such that t(X(s), y(s)) = T for all s£i We 
proceed by iteration and compute a set of grid points that approximates well such a curve, 
for example, by taking 

i(k) = sup{i G {0, . . . , N} | t(X ijk , Y ijk ) < T}, 

m = k, 

for k = 0, . . . , N. For a given T, the function u(T, x) can be seen as the curve (x, u(T, x)) 
in M 2 which is parametrized by x £ M, and we approximate this curve by the points 

(^(^(fe),i(jfc),^i(fc),i(fc))) ^PQ(fc),i(fe)>^(fc),j(fc))) 
for k = 0, . . . ,N. This method has been used to produce the results presented in Figured) 




Figure 9. Plot of the isotimes, that is, the curves for which t(X, Y) is a 
constant. Note the box in the middle in which t is constant and equal to 
zero. 
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